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Overview 


Nonpremixed  jet  flames  have  been  extensively  studied  to  understand  the  combustion  processes  in 
rocket  engines.  The  stabilization  and  structure  of  jet  flames  determine  the  lift-off  height  of  the  flame 
and  are  therefore  integral  to  engine  design.  In  view  of  the  nonpremixed  nature  of  the  fuel  and  oxidizer 
jets,  there  are  two  dominant  modes  of  flame  stabilization.  The  traditional  view  is  that  mixing  between 
the  two  jets  will  eventually  lead  to  a  region  of  strong  reactivity  and  hence  auto-ignition.  An  alternate 
view,  advanced  recently  and  for  largely  non-autoignitive  situations,  is  stabilization  through  the  so-called 
tribrachial  flame  (also  known  as  triple  flame),  in  which  a  lean  and  a  rich  premixed  flame  wing  with  a 
trailing  diffusion  flame  branch.  The  point  where  the  three  branches  intersect,  the  triple  point,  is 
considered  to  be  the  stabilization  point.  The  dynamic  balance  between  the  local  flame  propagation 
speed  and  the  incoming  flow  speed  is  the  stabilization  mechanism. 

Regardless  of  the  mechanisms,  assessment  of  flame  stabilization  must  necessarily  include  the 
simultaneous  and  coupled  considerations  of  the  fluid  flow  in  the  combustion  chamber,  the  mixing 
between  the  fuel  and  oxidizer  jets,  and  the  chemistry  and  dynamics  either  of  the  mixing  layer  or  of  the 
flame.  In  response  to  these  considerations,  a  comprehensive  research  program  comprised  of 
experiment,  computation  and  theory,  has  been  conducted,  yielding  both  qualitative  understanding  and 
quantitative  data  on  the  various  issues  of  interest. 


Research  Accomplishments 

Broadly,  the  research  has  been  conducted  along  two  thrusts,  namely:  (1)  the  dynamics  of  flame 
stabilization  in  environments  of  various  thermodynamic  and  dynamics  characterizations,  and  (2)  the 
chemistry  and  dynamics  of  strained  flames,  particularly  those  involving  the  cool  flame,  NTC  (negative 
temperature  coefficient)  chemistry,  noting  that  interests  in  cool  flame  chemistry  and  dynamics  have 
been  a  major  activity  in  recent  fundamental  combustion  research.  Highlights  of  the  accomplishments 
are  summarized  next.  A  total  of  six  papers  reporting  the  work  were  published  in  leading  journals,  and 
are  appended  herein. 

1.  Dynamics  of  Flame  Stabilization 

1.1  Stabilization  of  nonpremixed  coflow  flames:  The  structure  and  stabilization  mechanism  of 
laminar  nonpremixed  autoignitive  DME/air  coflow  flames  were  investigated.  Computations  were 
performed  at  30  atmospheres  with  uniform  inlet  velocities  of  3.2  m/s  for  both  streams,  and  the  coflow 
air  boundary  temperatures  were  700,  800,  900,  and  1100  K.  The  heat  release  rate  and  species  profiles 
were  examined  for  each  case.  Further  investigation  with  Chemical  Explosive  Mode  Analysis  (CEMA)  and 
Lagrangian  Flamelet  Analysis  (LFA)  were  performed  to  identify  the  controlling  chemistry  and  elucidate 
the  dominant  combustion  mode  and  stabilization  mechanism.  At  700  -  900  K,  autoignition  was  observed 
to  be  the  dominant  stabilization  mechanism,  and  NTC  chemistry  determines  the  stabilization  point  in 
mixture  fraction  space.  Conversely,  at  1100  K,  the  kinematic  balance  between  the  premixed  flame 
propagation  velocity  and  the  incoming  flow  velocity  becomes  the  dominant  stabilization  mechanism, 
and  the  classical  triple  flame  structure  was  observed.  Extended  stabilization  regimes,  in  terms  of 
increasing  boundary  temperature,  are  therefore  identified,  including  frozen  flow,  kinetically  stabilized, 
autoignition-propagation-coupled  stabilized,  kinematically  stabilized,  and  burner  stabilized  regimes. 

This  work  is  reported  in  Publication  #1:  "Autoignition-affected  stabilization  of  laminar 
nonpremixed  DME/air  coflow  flames,"  by  Sili  Deng,  Peng  Zhao,  Michael  E.  Mueller  and  Chung  K.  Law, 
Combustion  and  Flame,  Vol.  162,  pp.  3437-3445  (2015). 
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1.2  Dependence  of  Stabilization  on  Pressure  and  Co-flow  Temperature:  The  structure  and 
stabilization  mechanism  of  the  nonpremixed  autoignitive  coflow  flames  in  1.1  was  investigated  at 
elevated  temperatures  and  pressures.  Computations  with  detailed  chemistry  were  performed  at  30  atm 
with  uniform  inlet  velocities  imposed  for  both  streams.  The  heat  release  rate  profiles  were  first 
examined  for  each  case  to  demonstrate  a  multibrachial  thermal  structure.  Species  concentrations  and 
temperature  were  sampled  along  mixture  fraction  iso-contours,  and  CEMA  was  performed  to  identify 
the  controlling  chemistry  at  representative  points.  One-dimensional  LFA  was  also  performed  and 
compared  with  the  two-dimensional  computations  to  elucidate  the  relative  importance  of  diffusion 
processes  parallel  and  normal  to  the  mixture  fraction  gradient.  Various  coflow  temperatures  with 
different  inlet  velocities  are  examined  to  elucidate  their  influences  on  the  multibrachial  structure  as  well 
as  the  stabilization  mechanism.  NTC-affected  inhomogeneous  autoignition  and  the  coupled  effects  with 
premixed  flame  propagation  on  stabilization  are  further  studied.  It  is  found  that,  at  high  coflow 
boundary  temperatures  or  low  inlet  velocities,  the  classical  tribrachial  flame  structure  is  achieved,  and 
autoignition  contributes  less  to  the  stabilization  due  to  reduced  heat  and  radical  accumulation.  The 
kinematic  balance  between  the  local  flow  speed  and  flame  propagation  speed  is  found  to  be  the 
dominant  stabilization  mechanism.  On  the  contrary,  kinetic  stabilization  is  achieved  at  lower  coflow 
temperatures  or  higher  inlet  velocities  as  autoignition  becomes  dominant.  Due  to  the  transition  of  the 
dominant  chemical  pathways  during  autoignition,  the  kinetically  stabilized  structure  is  usually 
multibrachial.  The  transition  of  different  stabilization  mechanisms  can  be  made  by  changing  either  the 
boundary  velocity  or  temperature  of  the  coflow.  Based  on  these  results,  a  regime  diagram  is  constructed 
that  identifies  the  possible  stabilization  regimes:  blowout,  kinetically  stabilized,  autoignition- 
propagation-coupled  stabilized,  kinematically  stabilized,  and  burner  stabilized. 

This  work  is  reported  in  Publication  #2:  "Stabilization  of  laminar  nonpremixed  DME/air  coflow 
flames  at  elevated  temperatures  and  pressures,"  by  Sili  Deng,  Peng  Zhao,  Michael  E.  Mueller  and  Chung 
K.  Law,  Combustion  and  Flame,  Vol.  162,  pp.  4471-4478  (2015). 

1.3  Stabilization  in  Oscillating  Flows:  The  above  study  was  then  extended  to  stabilization  in  flows 
with  sinusoidally  oscillating  inlet  velocities  in  order  to  identify  effects  of  oscillation  frequency  on  the 
flame  dynamics.  To  benchmark  the  unsteady  cases,  a  normalized  displacement  velocity  was  defined  to 
differentiate  flame  propagation  from  autoignition,  and  this  definition  was  validated  against  the  steady 
cases.  In  the  oscillating  reacting  flow,  transition  between  a  multibrachial  autoignition  front  and  a 
tribrachial  flame  occurs  periodically.  However,  unlike  the  harmonic  velocity  oscillation,  the  combustion 
mode  transition  is  hysteretic.  The  oscillation  cycle  starts  with  the  largest  inlet  velocity,  with  the 
multibrachial  thermal  structure,  located  downstream,  being  governed  by  autoignition  chemistry.  As  flow 
velocity  decreases,  the  autoignition  front  moves  upstream  and  transitions  to  a  tribrachial  flame  near  the 
lower  velocity  limit,  similar  to  the  steady  flow,  as  autoignition  chemistry  becomes  weaker  with 
decreasing  upstream  residence  time.  As  the  flow  velocity  increases  again,  the  tribrachial  flame  is 
convected  downstream  and,  ultimately,  due  to  the  radical  and  heat  accumulation  in  time,  autoignition 
eventually  occurs  and  becomes  the  dominant  pathway.  The  finite  induction  time  for  autoignition  results 
in  the  hysteretic  behavior  during  the  decreasing-  and  increasing-velocity  cycles.  The  hysterisis  diminishes 
at  lower  oscillation  frequency  as  there  is  more  time  for  chemistry  to  respond  to  the  hydrodynamic 
changes  and  consequently  approach  steady  state.  At  the  relatively  low  oscillation  frequencies 
investigated,  first-stage  NTC  chemistry  is  less  affected  by  flow  dynamics  with  only  second-stage 
autoignition  and  flame  chemistry,  which  accounts  for  the  majority  heat  release,  coupled  with  the  flow 
oscillation. 

This  work  is  reported  in  Publication  #3:  "Flame  dynamics  in  oscillating  flows  under  autoignitive 
conditions,"  by  Sili  Deng,  Peng  Zhao,  Michael  E.  Mueller  and  Chung  K.  Law,  Combustion  and  Flame,  Vol. 
168,  pp.  75-82  (2016). 
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1.4  Analysis  of  Stabilization  of  Triple  Flame  in  Oscillating  Flows:  Combustion  instability  of  a  lifted 
triple  flame  downstream  of  the  injector  nozzle  in  an  acoustic  field  is  studied.  The  analysis  is  focused  on 
the  premixed  flame  segment,  and  is  based  on  the  linear  equation  for  a  perturbed  flamefront.  Assuming 
that  the  radius  of  the  flamefront  curvature  is  much  larger  than  that  of  the  internal  flame  structure,  the 
premixed  front  is  approximated  as  being  quasi-planar.  Both  acoustic  stabilization  of  the  hydrodynamic 
flame  instability  and  excitation  of  the  parametric  instability  by  sound  waves  are  considered,  and  the 
limitation  for  stable/unstable  regimes  as  a  function  of  flame  parameters  and  acoustic 
frequency/intensity  is  determined.  The  velocity-wise  nature  of  the  sound  is  scrutinized,  namely,  the 
acoustics  are  imitated  by  a  pocket  of  harmonic  flow  oscillations.  First,  the  Landau  limit  of  infinitely  thin 
flames  is  employed.  Then  the  formulation  is  extended  to  account  for  the  finite  flame  thickness,  and  it  is 
demonstrated  how  unstable  regimes  are  modified  by  the  flame  thickness. 

This  work  is  reported  in  Publication  #4.  "Coupling  of  harmonic  flow  oscillations  to  combustion 
instability  in  premixed  segments  of  triple  flames,"  by  V'yacheslav  Akkerman  and  Chung  K.  Law, 
Combustion  and  Flame,  Vol.  172,  pp.  342-348  (2016). 

2.  Cool  Flame  Chemistry  in  Flame  Stabilization 

2.1  Existence  and  Propagation  of  Premixed  Cool  Flames:  Cool  flames,  being  essential  features  of 
chemical  kinetics  of  large  hydrocarbon  fuels,  are  closely  related  to  the  negative  temperature  coefficient 
phenomenon  and  engine  knock.  In  this  work,  the  coupling  of  cool  flame  chemistry  and  convective- 
diffusive  transport  is  computationally  and  experimentally  investigated.  A  1-D  planar  premixed  cool  flame 
induced  by  a  hot  pocket  is  first  simulated  for  DME/02/N2  mixtures  with  detailed  chemistry  and 
transport,  demonstrating  the  existence  of  a  residence  time  window  for  quasi-steady  propagation.  Then 
with  residence  time  limited  by  aerodynamic  straining,  a  steady-state  premixed  cool  flame  is  simulated  in 
a  counterflow  of  heated  N2  against  a  DME/02/N2  mixture.  It  is  found  that  with  a  high  strain  rate, 
corresponding  to  short  residence  time,  low-temperature  heat  release  is  suppressed,  resulting  in  a 
stretched  low-temperature  S-curve  system  response;  and  that  with  a  sufficiently  low  strain  rate, 
corresponding  to  long  residence  time,  ignition  induced  by  low-temperature  chemistry  would  transition 
to  a  high-temperature,  intensely  burning  flame.  Consequently,  a  steady-state  premixed  cool  flame  exists 
only  for  residence  time  in  a  strain  rate  window.  A  symmetric  counterflow  configuration  is  then 
simulated  to  determine  the  cool  flame  temperature  and  flame  speed  at  a  fixed  local  strain  rate,  showing 
very  different  controlling  chemistry  and  characteristics  as  compared  to  the  normal  laminar  flames 
governed  by  high-temperature  chemistry.  In  a  companion  experimental  investigation,  premixed  cool 
flames  in  the  counterflow  were  observed  with  a  high-sensitivity  CCD  camera  in  the  UV  spectrum, 
with/without  a  bandpass  filter  corresponding  to  the  characteristic  wavelength  of  excited  HCHO.  The 
chemiluminescence  from  the  cool  flame  is  found  to  become  more  intense  with  increasing  equivalence 
ratio,  even  for  rich  mixtures,  while  the  position  of  the  cool  flame  is  insensitive  to  variation  in  the 
equivalence  ratio  at  the  same  strain  rate.  These  observations  qualitatively  agree  with  the  numerical 
simulations,  demonstrating  the  essential  features  of  premixed  cool  flames. 

This  work  is  reported  in  Publication  #5:  "Initiation  and  propagation  of  laminar  premixed  cool 
flames,"  by  Peng  Zhao,  Wenkai  Liang,  Sili  Deng  and  Chung  K.  Law,  Fuel,  Vol.  166,  pp.  477-487  (2016). 

2.2  Ignition  and  Extinction  of  Nonpremixed  Cool  Flames:  Following  2.1,  and  recognizing  that  the 
low-temperature  chemical  kinetics  is  promoted  at  elevated  pressures,  the  ignition  and  extinction  of  non¬ 
premixed  cool  flames  at  elevated  pressures  were  experimentally  and  computationally  investigated  in 
the  counterflow.  Specifically,  the  hysteretic  ignition  and  extinction  behavior  of  the  nonpremixed  cool 
flame  was  for  the  first  time  observed  and  quantified.  S-curve  analysis  was  conducted  to  demonstrate 
the  thermal  and  chemical  structure  of  the  cool  flame  and  to  elucidate  the  dominant  chemical  pathways 
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during  the  ignition  and  extinction  processes.  The  dominant  low-temperature  chemical  reactions  shift 
from  those  responsible  for  radical  runaway  to  exothermic  reactions  that  sustain  the  cool  flame. 
Increasing  the  ambient  pressure  and/or  the  oxygen  concentration  in  the  oxidizer  stream  promotes  the 
heat  release  from  the  cool  flame,  and  hence  results  in  an  extended  hysteresis  temperature  window 
between  ignition  and  extinction.  It  is  further  noted  that  while  the  observed  cool  flame  ignition 
temperatures  were  well  predicted  by  computation,  significant  discrepancies  existed  for  the  extinction 
temperatures  based  on  the  well-adopted  reaction  mechanism  used. 

This  work  is  reported  in  Publication  #6:  "Ignition  and  extinction  of  strained  nonpremixed  cool 
flames  at  elevated  pressures,"  by  Sili  Deng,  Dong  Han  and  Chung  K.  Law,  Combustion  and  Flame,  Vol. 
176,  pp.  143-150  (2017). 
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The  structure  and  stabilization  mechanism  of  laminar  nonpremixed  autoignitive  DME/air  coflow  flames 
were  investigated.  Computations  were  performed  at  30  atmospheres  with  uniform  inlet  velocities  of 
3.2  m/s  for  both  streams,  and  the  coflow  air  boundary  temperatures  were  700,  800,  900,  and  1100K. 
The  heat  release  rate  and  species  profiles  were  examined  for  each  case.  Further  investigation  with 
Chemical  Explosive  Mode  Analysis  (CEMA)  and  Lagrangian  Flamelet  Analysis  (LFA)  were  performed  to 
identify  the  controlling  chemistry  and  elucidate  the  dominant  combustion  mode  and  stabilization  mech¬ 
anism.  At  700-900  K,  autoignition  was  observed  to  be  the  dominant  stabilization  mechanism,  and  NTC 
chemistry  determines  the  stabilization  point  in  mixture  fraction  space.  Conversely,  at  1100  K,  the  kine¬ 
matic  balance  between  the  premixed  flame  propagation  velocity  and  the  incoming  flow  velocity  becomes 
the  dominant  stabilization  mechanism,  and  the  classical  triple  flame  structure  was  observed.  Extended 
stabilization  regimes,  in  terms  of  increasing  boundary  temperature,  are  therefore  identified,  including 
frozen  flow,  kinetically  stabilized,  autoignition-propagation-coupled  stabilized,  kinematically  stabilized, 
and  burner  stabilized  regimes. 

©  2015  The  Combustion  Institute.  Published  by  Elsevier  Inc.  All  rights  reserved. 


1.  Introduction 

Nonpremixed  jet  flames  have  been  extensively  studied  to 
understand  the  combustion  processes  in  rocket  and  diesel  engines. 
The  stabilization  and  structure  of  jet  flames  determine  the  lift-off 
height  of  the  flame  and  are  therefore  integral  to  engine  design. 
Due  to  the  mixing  process  of  the  fuel  and  oxidizer  streams  in  lifted 
flames  at  nonautoignitive  conditions,  the  combustion  mode  is  par¬ 
tially  premixed,  leading  to  the  observation  of  a  two-dimensional 
tribrachial  flame  (also  known  as  triple  flame)  [1];  specifically,  a 
lean  and  a  rich  premixed  flame  wing  with  a  trailing  diffusion  flame 
branch.  The  point  where  the  three  branches  intersect  is  called  the 
triple  point  and  is  generally  considered  to  be  the  stabilization  point 
for  nonautoignitive  situations.  The  dynamic  balance  between  the 
local  flame  propagation  speed  and  the  incoming  flow  speed  is  char¬ 
acterized  as  the  stabilization  mechanism.  A  recent  review  by 
Chung  [2]  discussed  the  stabilization,  propagation,  and  instability 
of  tribrachial  flames,  including  the  effects  of  concentration  gradi¬ 
ent  [3-5],  velocity  gradient  [6],  and  burned  gas  expansion  [7-10]. 
These  studies,  however,  were  limited  to  nonautoignitive  condi¬ 
tions,  but  real  engines  are  operated  at  elevated  pressures  and 
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temperatures,  where  autoignition  is  activated  and  could  interact 
with  the  tribrachial  flame. 

Chung  and  co-workers  1 1—13]  further  conducted  a  series  of 
experiments  to  investigate  the  autoignition  characteristics  of  lam¬ 
inar  Ci  to  C4  fuel  jets  in  a  heated  air  coflow  and  found  that,  above 
certain  coflow  temperatures,  lifted  flames  could  be  established 
through  autoignition.  In  these  studies,  both  the  tribrachial  struc¬ 
ture  for  most  autoignited  cases  and  a  repetitive  behavior  of  extinc¬ 
tion  and  reignition  at  the  critical  condition  near  blowout  were 
observed.  However,  the  role  that  autoignition  plays  in  the  stabi¬ 
lization  mechanism  as  well  as  its  influences  on  the  tribrachial 
flame  structure  are  still  less  understood. 

Furthermore,  practical  hydrocarbon-based  fuels  generally  have 
two-stage  ignition  processes,  in  which  the  first  stage  ignition  is 
governed  by  low  temperature  chemistry  and  the  second  stage  igni¬ 
tion  by  high  temperature  chemistry.  In  both  low  and  high  temper¬ 
ature  regimes,  the  ignition  delay  time  decreases  as  the  initial 
temperature  increases.  However,  in  the  intermediate  temperature 
regime,  the  transition  of  the  ignition  chemistry  results  in  increased 
overall  ignition  delay  time  as  the  initial  temperature  increases, 
exhibiting  the  negative  temperature  coefficient  (NTC)  phenomena, 
which  has  been  extensively  studied  in  homogeneous  systems  as  a 
major  feature  of  large  hydrocarbon  autoignition  [14].  For  engine 
applications,  however,  the  coupling  between  NTC  chemistry  and 
transport  processes  should  be  considered,  for  nonuniformities 
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invariantly  exist  in  realistic  combustion  systems.  When  the  trans¬ 
port  time  scale  becomes  comparable  to  that  of  the  NTC  chemical 
time  scale,  the  two  processes  are  expected  to  be  strongly  coupled. 
As  a  consequence,  the  global  response  of  the  inhomogeneous  sys¬ 
tem  can  also  be  affected  by  NTC  chemistry.  Recently,  a  series  of 
computational  and  experimental  studies  adopting  the  non- 
premixed  counterflow  configuration  by  Law  and  co-workers  [15- 
17]  have  demonstrated  that,  with  the  existence  of  nonuniformities 
in  the  flow,  species,  and  temperature  fields,  the  ignition  character¬ 
istics  of  nonpremixed  flames  can  be  fundamentally  affected  by  NTC 
effects,  especially  at  elevated  pressures  and/or  reduced  strain  rates. 

Therefore,  NTC-affected  stabilization  of  nonpremixed  lifted  jet 
flames  can  be  potentially  important,  yet  few  literatures  provide 
detailed  analysis.  Krisman  et  al.  [18]  recently  conducted  a  numer¬ 
ical  study  of  dimethyl  ether  (DME)/air  mixing  layer  at  40  atmo¬ 
spheres  and  air  coflow  temperatures  ranging  from  700  to  1500K 
and  observed  multibrachial  structures  in  the  heat  release  rate  pro¬ 
files.  The  mixture  fractions  corresponding  to  the  stabilization 
points  defined  based  on  the  hydroxyl  radical  (OH)  mass  fraction 
and  the  first  stage  autoignition  kernels  based  on  the 
methoxymethylperoxy  radical  (CH30CH202)  were  compared  with 
the  most  reactive  mixture  fractions  computed  from  homogeneous 
autoignition  under  the  same  initial  conditions.  A  transport  budget 
analysis  based  on  selected  species  was  performed  to  differentiate 
deflagration  from  autoignition. 

In  light  of  the  reported  multibrachial  structure,  showing  a  mod¬ 
ified  flame  shape  from  autoignition  in  the  mixing  layer,  further 
investigation  is  warranted  to  identify  the  detailed  chemical  struc¬ 
ture  and  stabilization  mechanism  of  the  multibrachial  flame.  For 
example,  tools  for  computational  diagnostics,  especially  for  identi¬ 
fying  locally  dominant  chemical  reactions,  can  be  employed  to 
understand  the  controlling  chemistry.  Moreover,  a  direct  compar¬ 
ison  to  homogeneous  autoignition  is  insufficient  to  understand 
the  transport  processes  in  the  current  configuration.  In  the 
two-dimensional  mixing  layer,  transport  processes  in  two  direc¬ 
tions  are  important:  parallel  and  normal  to  the  mixture  fraction 
gradient,  which  are  due  to  transverse  stratification  of  temperature 
and  species  and  streamwise  flow  and  (flame  back)  diffusion, 
respectively.  These  considerations  would  significantly  improve 
the  understanding  of  the  role  of  autoignition  upstream  of  the  flame 
structure  and  quantitatively  identify  the  controlling  kinetics  and 
stabilization  mechanism. 

In  the  present  study,  nonpremixed  DME/air  coflow  flames  were 
computed  at  30  atmospheres  with  the  oxidizer  stream  heated  to 
activate  autoignition.  With  fixed  inlet  velocities,  only  the  oxidizer 
stream  boundary  temperature  was  varied  to  investigate  the  corre¬ 
sponding  lifted  flame  morphology,  chemical  structure,  and  domi¬ 
nant  reaction  pathways.  In  the  following  we  shall  first  present 
the  computational  details  of  the  study.  The  thermal  and  chemical 
structures  are  then  described  with  heat  release  rate  and  selected 
species  profiles  in  Section  3.  The  evolution  of  the  controlling  chem¬ 
ical  pathways  are  subsequently  identified  with  Chemical  Explosive 
Mode  Analysis  (CEMA)  and  the  stabilization  mechanism  deter¬ 
mined  with  Lagrangian  Flamelet  Analysis  (LFA)  in  Section  4. 
Finally,  the  transition  of  the  dominant  stabilization  mechanism  is 
analyzed  in  Section  5,  with  extended  stabilization  regimes  con¬ 
structed  for  completeness. 


adiabatic  slip  wall,  and  its  diameter  is  large  enough  such  that 
increasing  the  width  of  the  domain  does  not  influence  the  compu¬ 
tation.  Uniform  inlet  velocities  of  3.2  m/s  were  specified  for  both 
fuel  and  air  streams  and  kept  the  same  for  all  the  cases  to  establish 
lifted  flames.  For  the  unsteady  establishment  of  the  flame,  a  con¬ 
vective  outflow  is  utilized  at  the  outlet  boundary,  which  simplifies 
to  a  Neumann  condition  for  the  steady  problem. 

The  flow  field  was  initialized  on  a  coarse  mesh  within  a  large 
domain.  At  first,  all  the  chemical  source  terms  were  set  to  zero 
until  the  nonreacting  flow  reached  steady-state.  Chemical  source 
terms  were  then  activated;  the  mixture  autoignited;  and  the  flow 
field  reached  steady-state.  The  domain  was  then  truncated,  and 
the  mesh  was  refined  to  fully  resolve  the  chemical  structure.  All 
the  results  presented  in  the  present  work  were  obtained  from 
the  steady-state  solutions. 

The  Navier-Stokes  equation  with  buoyancy  in  the  streamwise 
direction  and  the  conservation  equations  of  mass,  species,  and 
energy  were  solved.  The  species  diffusivities  were  determined 
from  a  constant,  nonunity  Lewis  number.  The  Lewis  numbers  for 
individual  species  (summarized  in  the  Supplementary  Data)  were 
pre-calculated  from  a  one-dimensional  flamelet  with  the  same 
boundary  conditions  and  the  mixture-averaged  transport  model 
and  evaluated  at  the  maximum  temperature  location.  The  con¬ 
served  scalar  mixture  fraction  Z  was  specified  as  unity  and  zero 
for  the  fuel  jet  and  coflow  at  the  inlet,  respectively,  and  computed 
by  solving  its  transport  equation  with  unity  Lewis  number  [19]. 
This  definition  of  mixture  fraction  is  consistent  with  the  one  used 
in  the  flamelet  calculation  in  Section  4.2. 

Dimethyl  ether  was  chosen  in  this  work,  for  it  is  a  clean  biofuel 
and  one  of  the  smallest  hydrocarbons  exhibiting  NTC  behavior. 
Detailed  reaction  models  for  low  and  high  temperature  DME  oxi¬ 
dation  [20-23]  have  been  developed  and  validated  in 
burner-stabilized  flames  [24],  nonpremixed  counterflow  ignition 
[25],  and  laminar  flame  speeds  [26].  The  present  computations 
were  conducted  using  a  skeletal  mechanism  of  39  species  27], 
including  both  low  and  high  temperature  oxidation  pathways, 
which  was  reduced  from  the  well  validated  detailed  mechanism 
of  Zhao  et  al.  [23]. 

The  low-Mach  number  formulation  of  the  governing  equations 
is  solved  using  NGA,  which  is  based  on  the  numerical  methods  of 
Desjardins  et  al.  [28].  The  momentum  and  scalar  equations  are  dis¬ 
cretized  with  a  second-order  centered  scheme  and  a  third-order 
WENO  scheme  [29],  respectively,  on  a  staggered  mesh.  The  itera¬ 
tive  second-order  semi-implicit  Crank-Nicolson  scheme  of  Pierce 
and  Moin  [30]  is  adopted  for  temporal  integration.  At  each  time 
step,  the  chemical  source  terms  for  the  species  and  energy  equa¬ 
tions  are  evaluated  independently  from  the  transport  terms  using 
the  CVODE  package  [31]. 

Uniform  grids  in  the  axial  direction  were  adopted  for  the  com¬ 
putations,  and  the  grid  spacing  was  set  as  Ax  =  2.2  pm.  A  nonuni¬ 
form  grid  was  used  in  the  radial  direction  with  a  minimum  spacing 
of  2.5  pm  to  resolve  the  mixing  layer  corresponding  to  the  separa¬ 
tion  wall  and  geometric  progression  stretch  rates  less  than  3% 
towards  both  the  centerline  and  the  outer  boundary.  The  dimen¬ 
sions  of  and  number  of  grid  points  in  the  computational  domain 
for  each  computation  are  summarized  in  Table  1. 


2.  Computational  details 


Table  1 

Computational  domain  and  number  of  grid  points. 


The  flow  configuration  is  an  axisymmetric  DME  stream  at  300  K 
in  a  heated  coflow  of  air  (700,  800,  900,  and  1100  K)  at  30  atmo¬ 
spheres.  The  fuel  nozzle  diameter  D  is  0.8  mm,  and  the  fuel  and 
air  are  initially  separated  with  an  adiabatic,  no-slip  wall  with 
thickness  D/20.  The  coflow  outer  boundary  is  specified  as  an 


Coflow  temperature  [K] 

700 

800 

900 

1100 

Lx  [mm] 

28 

7 

3.5 

3 

Lr  [mm] 

6 

3.9 

3.9 

6 

Nx 

12,290 

3072 

1536 

1282 

Nr 

192 

176 

176 

192 
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Fig.  1.  Velocity,  temperature,  and  H202  profiles  along  Zst  on  the  nominal  and  two  times  coarser  (in  each  direction)  meshes  for  an  air  temperature  of  800  K. 


A  grid  convergence  study  was  performed  for  the  air  tempera¬ 
ture  800  K  case,  for  it  has  the  most  complex  structure,  which  is  dis¬ 
cussed  in  the  following  sections.  As  shown  in  Fig.  1,  grid 
convergence  was  achieved  for  velocity,  temperature,  and  species 
profiles.  Grid  convergence  was  also  verified  for  the  air  temperature 
1100K  case,  which  shows  similar  results  and  is  therefore  not 
shown  here. 


3.  Thermal  and  chemical  structure 

To  visualize  the  flame  structures,  the  heat  release  rate  profiles 
for  the  four  cases  (700,  800,  900,  and  1100  K)  are  shown  in  Fig.  2. 
Qualitatively,  the  most  upstream  point  on  the  largest  heat  release 
contour  (the  leading  point),  colored  by  red,  will  be  referred  to  as 
the  stabilization  point. 

At  700  K,  a  tribrachial  thermal  structure  is  observed,  and  the 
stabilization  point  is  located  around  Z  =  0.15,  which  is  richer  than 
the  triple  point,  where  the  three  branches  intersect.  Moreover, 
compared  to  the  classical  triple  flame  structure,  the  middle  heat 
release  rate  branch,  corresponding  to  the  nonpremixed  flame,  is 
significantly  weaker  than  the  other  two  branches. 

At  800  K,  the  stabilization  point  is  not  located  on  the  tribrachial 
structure  any  more.  Instead,  it  is  located  near  Z  =  0.23  and  con¬ 
nects  two  trailing  heat  release  branches,  where  a  tribrachial  flame 
structure  is  attached  to  the  leaner  branch  (LB)  of  the  bibrachial 
reacting  front.  A  schematic  of  the  structure  is  shown  in  Fig.  3. 

As  the  air  boundary  temperature  increases  to  900  K,  the  stabi¬ 
lization  point  shifts  back  to  Z  =  0.14.  Moreover,  a  long  trailing 
branch  at  richer  mixture  fraction  is  attached  to  the  main  tribrachial 
structure,  resulting  in  a  tetrabrachial  structure.  Compared  with  the 
structure  shown  in  the  800  K  case,  the  main  tribrachial  structure 
stabilizes  further  upstream,  for  it  depends  less  on  the  radical  accu¬ 
mulation  ahead  of  the  flame.  Therefore,  it  catches  up  with  the 
reacting  front  at  richer  mixture  fraction,  and  they  merge  into  the 
apparent  tetrabrachial  structure. 

A  further  increase  in  the  boundary  temperature  results  in  a 
structure  that  is  very  similar  to  the  classical  triple  flame,  except 
for  the  fact  that  there  is  also  heat  release  ahead  of  the  stabilization 
point  at  Z  =  0.13.  Some  of  the  multibrachial  structures  were  also 
observed  by  Krisman  et  al.  [18],  using  different  definitions  for 
branches,  and  it  was  concluded  that  the  autoignition  chemistry 
could  affect  the  flame  structure  and  the  stabilization  mechanism. 

To  first  qualitatively  demonstrate  the  chemical  structure  of  the 
flame,  selected  species  profiles  were  examined,  shown  in  Figs.  4-7. 
The  methoxymethylperoxy  radical  (CH3OCH202)  and  hydroxyl  rad¬ 
ical  (OH)  were  chosen  as  indicators  of  low  and  high  temperature 
chemistry,  respectively.  The  hydroperoxyl  radical  (H02)  and  hydro¬ 
gen  peroxide  (H202)  were  chosen,  for  they  form  in  the  preheat 


zone  of  a  flame  or  before  autoignition  but  quickly  vanish  in  the 
post  flame  zone  or  after  ignition  [32]. 

For  all  four  cases,  similar  profiles  can  be  seen  for  some  species. 
First,  low  temperature  chemistry,  indicated  by  the  CH30CH202  rad¬ 
ical,  is  found  to  be  important  at  richer  mixture  fractions,  where  the 
temperature  is  also  lower.  Second,  the  OH  radical  peaks  at  and 
downstream  of  the  maximum  heat  release  locations  and  correlates 
well  with  the  tribrachial  structure  shown  in  the  heat  release  rate 
profiles,  indicating  the  presence  of  high  temperature  chemistry. 
Third,  the  H02  mass  fraction  peaks  in  a  thin  region.  Compared  with 
the  heat  release  rate  contours,  this  thin  region  outlines  the  flame 
front  and  the  reactive  mixture  at  the  rich  mixture  fractions  and 
indicates  the  importance  of  the  exothermic  three-body  recombina¬ 
tion  reaction  H  +  02  +  M  H02  +  M. 

However,  there  are  also  differences  in  the  chemical  structure 
among  the  different  cases.  For  example,  for  the  800  and  900  K 
cases,  another  OH  local  maxima,  which  is  two  orders  of  magni¬ 
tudes  smaller  than  the  peak  value  on  the  tribrachial  structure, 
appears  at  richer  mixture  fractions,  immediately  downstream  of 
where  the  CH30CH202  radical  and  H202  disappear,  indicating 
autoignition.  Moreover,  more  pronounced  differences  between 
the  three  lower  boundary  temperature  cases  and  the  1100  K  case 
are  shown  in  the  H202  profiles:  for  the  lower  boundary  tempera¬ 
ture  cases,  H202  accumulates  along  the  mixture  fraction 
iso-contours  until  it  decomposes  in  the  flame  region,  while,  for 
the  1100K  case,  the  H202  accumulation  is  an  order  of  magnitude 
lower,  due  to  the  reduced  residence  time  from  the  nozzle  exit  to 
the  flame  base. 

4.  Computational  diagnostics  and  analysis 

The  above  heat  release  rate  and  species  profiles  demonstrate 
the  thermal  and  chemical  structure  of  the  reacting  fronts  at  differ¬ 
ent  boundary  temperatures.  However,  more  detailed  computa¬ 
tional  diagnostics  and  analysis  are  needed  to  further  demonstrate 
the  controlling  chemistry  and  the  stabilization  mechanism. 

4A.  Chemical  explosive  mode  analysis 

In  addition  to  the  analysis  based  on  selected  species  profiles, 
Chemical  Explosive  Mode  Analysis  (CEMA)  [33,34]  was  conducted 
to  identify  the  controlling  chemistry  in  these  complex  reacting 
flows.  Briefly,  the  eigenvalues  of  the  Jacobian  matrix  of  the  chem¬ 
ical  source  term  based  on  the  local  species  concentrations  and  tem¬ 
perature  are  evaluated  and  determined  as  the  chemical  modes.  The 
largest  real  part  of  all  eigenvalues,  which  is  defined  as  the  chemical 
explosive  mode,  describes  the  rate  of  system  runaway.  The  normal¬ 
ized  contribution  of  each  reaction  to  the  chemical  explosive  mode 
is  quantified  with  a  participation  index  [34]. 

In  the  present  study,  the  dominant  reactions  at  representative 
locations,  such  as  those  upstream  and  near  the  flame  base,  are 
identified,  based  on  the  explosive  mode  and  participation  index. 


1  For  interpretation  of  color  in  Fig.  2,  the  reader  is  referred  to  the  web  version  of 
this  article. 
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Fig.  2.  Heat  release  rate  [J/m3  s]  profiles.  The  iso-contours  of  Zst,  Z  =  0.2,  and  Z  =  0.3  are  outlined  from  right  to  left  in  solid  lines,  respectively.  The  CEMA  sampling  points  at 
800  and  1100  K  are  indicated  along  the  iso-contours. 


Fig.  3.  A  schematic  of  the  thermal  structure  of  the  800  K  case.  LPF,  RPF,  and  NPF 
denotes  the  lean  premixed,  rich  premixed,  and  nonpremixed  flame  branches  on  the 
tribrachial  structure,  respectively.  LB  and  RB  denotes  the  leaner  and  richer  branches 
of  the  reacting  front,  respectively. 

For  each  case,  the  local  species  concentrations  and  temperature 
were  sampled  along  the  Zst,  Z  =  0.2,  and  Z  =  0.3  iso-contours,  as 
indicated  in  Fig.  2,  and  processed  by  CEMA  to  demonstrate  the  evo¬ 
lution  of  the  dominant  reactions. 


The  results  are  summarized  in  Fig.  8,  where  three  representa¬ 
tive  locations  along  the  Zst  iso-contour  approaching  the  flame  front 
and  two  locations  ahead  of  and  at  the  reaction  front  at  Z  =  0.2 
were  sampled.  For  the  three  lower  coflow  temperature  cases,  sim¬ 
ilar  chemical  patterns  were  found.  Upstream  of  the  flame  front,  the 
chemical  explosive  mode  is  positive,  indicating  that  the  mixtures 
have  the  potential  to  explode;  downstream  of  the  flame  front, 
the  chemical  explosive  mode  becomes  negative,  meaning  that 
the  mixtures  are  composed  of  burned  products.  Following  the  Zst 
iso-contour,  the  hydrogen  peroxide  chain  branching  reaction 
(H202  +  M  4=^  OH  +  OH  +  M)  is  the  reaction  that  has  the  largest 
contribution  to  the  explosive  mode,  showing  the  dominant  role 
of  autoignition  chain  branching  [35].  The  characteristic  DME  low 
temperature  chemistry  is  also  important  upstream  of  the  flame, 
where  methoxymethylperoxy  radical  formation 
(CH3OCH2  +  02  4=^  CH30CH202)  and  isomerization 

(CH3OCH202  CH20CH202H)  promote  the  explosion,  while 
the  /2-scission  reaction  (CH20CH202H  OH  +  CH20  +  CH20) 
retards  the  explosion.  Approaching  the  flame  front,  the  H  radical 
recombination  reaction  (H  +  02  +  M  4=>  H02  +  M)  becomes 
important  for  the  700-900  K  cases,  due  to  the  fact  that  the  H  rad¬ 
icals  generated  at  the  reaction  zone  diffuse  upstream  and  undergo 
three-body  recombination  reactions  under  the  high  pressure,  low 
temperature  condition.  Further  downstream  where  the  heat 
release  rate  peaks,  the  hydrogen  branching  reaction 
(H  +  02  4=^  0  + OH)  becomes  the  most  important  chain  branching 
reaction  as  it  is  activated  at  high  temperatures  [35]. 

CEMA  conducted  along  the  Z  =  0.2  iso-contour,  which  crosses 
the  rich  heat  release  front  in  the  800  and  900  K  cases,  shows 
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Fig.  4.  Hydroxyl  (OH)  radical  mass  fraction  profiles.  The  mixture  fraction  iso-contours  are  the  same  as  in  Fig.  2. 
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Fig.  6.  Hydroperoxyl  (H02)  radical  mass  fraction  profiles.  The  mixture  fraction  iso-contours  are  the  same  as  in  Fig.  2. 
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Fig.  7.  Hydrogen  peroxide  (H202)  mass  fraction  profiles.  The  mixture  fraction  iso-contours  are  the  same  as  in  Fig.  2. 
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different  chemical  mode  evolution.  The  H202  chain  branching  reac¬ 
tion  is  always  the  dominant  reaction  that  promotes  the  explosive 
mode,  while  the  H  radical  recombination  reaction  and  the  H 
branching  reaction  are  less  important  ahead  of  the  rich  heat  release 
front  and  at  the  front. 

On  the  contrary,  although  low  temperature  chemistry  is  still 
important  for  the  1100  K  case  upstream  of  the  reaction  zone  and 
the  hydrogen  chain  branching  reaction  promotes  explosion  in  the 
reaction  zone,  the  hydrogen  peroxide  chain  branching  reaction  is 
not  very  important  for  all  the  sampled  locations.  Since  the  hydro¬ 
gen  peroxide  reaction  is  the  crucial  chain  branching  reaction  for 

DISTRIBUTION  A:  Distribution 


the  autoignition  process,  it  is  concluded  that  the  1 100 1<  case  is  less 
affected  by  autoignition  chemistry  than  the  lower  boundary  tem¬ 
perature  cases. 

4.2.  Lagrangian  flamelet  analysis 

The  above  species  profile  analysis  and  CEMA  results  have 
demonstrated  that  autoignition  chemistry  is  crucial  to  the  complex 
flame  structure  in  the  700-900  K  cases.  However,  the  role  that 
autoignition  plays  in  the  stabilization  still  needs  further  investiga¬ 
tion.  To  elucidate  the  role  of  autoignition  for  the  current  flow 
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Fig.  8.  Normalized  participation  index  at  800  K  and  1100  K.  Sampled  locations  are  indicated  in  Fig.  2. 
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configuration,  a  direct  comparison  with  the  homogeneous  counter¬ 
part  for  a  Lagrangian  flow  particle  is  insufficient,  for  no  transport 
process  is  considered  in  homogeneous  autoignition.  In  the  current 
two-dimensional  configuration,  however,  transport  processes  are 
present  in  two  directions:  parallel  and  normal  to  the  mixture  frac¬ 
tion  gradient.  These  transport  processes  might  be  important:  first, 
the  temperature  and  species  stratification  parallel  to  the  mixture 
fraction  gradient  can  significantly  modify  the  ignition  characteris¬ 
tics,  especially  for  fuels  with  NTC  chemistry  [15,17].  Second,  flame 
propagation  normal  to  the  mixture  fraction  gradient  can  also  influ¬ 
ence  the  autoignition  front  through  thermal  and  radical  back 
diffusion. 

To  demonstrate  the  dominant  transport  direction  as  well  as  the 
stabilization  mechanism  in  the  current  work,  one-dimensional 


unsteady  flamelet  analysis  was  conducted  to  account  for  unsteadi¬ 
ness  (convection),  chemical  reactions,  and  diffusion  parallel  to  the 
mixture  fraction  gradient,  while  neglecting  the  transport  process  in 
the  normal  direction.  As  a  consequence,  the  unsteady  flamelet  is 
able  to  capture  inhomogeneous  autoignition,  with  diffusion 
allowed  only  in  one  direction.  Following  the  mixture  fraction 
iso-contour,  the  spatial  information  from  the  two-dimensional 
computation  could  be  interpreted  as  the  time  history  of  the  corre¬ 
sponding  mixture  in  the  Lagrangian  frame.  If  the  one-dimensional 
unsteady  flamelet  predicts  this  time  history,  only  the  transport 
processes  parallel  to  the  mixture  fraction  gradient  are  important, 
and,  therefore,  the  thermal  structure  is  stabilized  by  inhomoge¬ 
neous  autoignition.  Conversely,  if  the  unsteady  flamelet  solutions 
do  not  agree  with  the  two-dimensional  computations,  the 
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transport  processes  along  the  mixture  fraction  iso-contour  are  not 
negligible  compared  to  the  gradient  direction.  Therefore,  premixed 
flame  propagation  is  the  dominant  stabilization  mechanism,  or,  at 
least,  stabilization  is  strongly  affected  by  flame  back  diffusion. 

In  the  present  work,  the  unsteady  flamelet  model  developed  by 
Pitsch  et  al.  [36],  referred  to  as  Lagrangian  Flamelet  Analysis  (LFA), 
was  adopted.  Due  to  the  mixing  processes,  the  scalar  dissipation 
rate  /,  which  can  influence  the  flamelet  solution  significantly, 
decreases  in  the  streamwise  direction.  Therefore,  this  dissipation 
rate  variation  must  be  considered  when  computing  a  flamelet  as 
it  evolves  downstream. 

In  the  present  study,  the  unsteady  flamelet  was  computed  with 
FlameMaster  [37],  and  the  dissipation  rate  was  specified  as  a  func¬ 
tion  of  the  flamelet  time.  The  flamelet  time  was  computed  from  the 
two-dimensional  computational  results,  along  the  stoichiometric 
mixture  fraction  Zst  iso-contour: 

[x  1 

t  =  J0  (u  +  uz)(x')|(Z  =  Zst)dX-  (1) 


This  formulation  is  otherwise  the  same  as  that  of  Pitsch  et  al.  36], 
except  that,  in  addition  to  the  axial  component  of  fluid  convection 
velocity  u ,  the  axial  component  of  the  mixture  fraction 
iso-contour  propagation  speed  relative  to  the  fluid  convection  uz 
is  also  taken  into  account.  The  expression  for  the  constant  property 
scalar  iso-surface  velocity  relative  to  the  local  fluid  motion  was 
derived  by  Pope  [38  ,  and  the  current  work  adopts  the  formulation 
derived  by  Lignell  et  al.  [39]  for  variable  properties: 


uz  =  - 


V  •  (pDzVZ) 
PlVZI 


n, 


(2) 


where  Dz  is  the  mixture  fraction  diffusivity,  which  is  defined  in 
Section  2,  and  p  the  density.  The  normal  vector  n,  defined  as 


VZ 

|VZ|’ 


(3) 


indicates  the  direction  of  this  diffusion  induced  relative  velocity. 
The  dissipation  rate  along  the  Zst  =  0.1005  iso-contour  obtained 
from  the  two-dimensional  computation  was  then  correlated  with 
this  flamelet  time  and  provided  as  the  input  for  the  LFA  calculation. 
The  dissipation  rates  at  other  mixture  fractions  were  computed 
assuming  the  following  form  40]: 


x(Z).xfe)expe2|etfc;'(2Z||';) 

exp(-2[erfc  (2 Zst)]  ) 


X(Zst)/(Z;Zst). 


(4) 


calculated  and  compared  with  the  two-dimensional  computational 
results  for  Zst,  Z  =  0.2,  and  Z  =  0.3. 

As  shown  in  Fig.  10,  two  ignition  stages  can  be  seen  at  Zst  and 
Z  =  0.2  for  the  700  K  case,  while  at  Z  =  0.3  only  one  ignition  dom¬ 
inated  by  low  temperature  chemistry  is  observed,  due  to  the 
reduced  initial  temperature.  At  all  three  mixture  fractions  exam¬ 
ined,  the  flamelets  agree  with  the  two-dimensional  computations 
very  well.  For  the  800  K  case,  both  the  flamelet  and 
two-dimensional  computation  experience  almost  identical  time 
histories,  where  two-stage  ignition  occurs  at  all  three  mixture  frac¬ 
tions.  As  the  initial  temperature  further  increases,  corresponding  to 
the  increase  in  the  boundary  temperatures  in  the  CFD  computa¬ 
tion,  the  two-stage  ignition  phenomenon  is  less  pronounced. 
However,  the  900  K  case  still  shows  good  agreement  between  the 
flamelet  profile  with  the  time  history  of  the  two-dimensional  com¬ 
putation  at  Z  =  0.2  and  0.3,  while  LFA  slightly  lags  behind  the  CFD 
computation  at  Zst,  similar  to  the  800  K  case. 

On  the  contrary,  for  the  1100K  case,  the  ignition  delay  time 
computed  with  the  one-dimensional  flamelet  assumption  is  signif¬ 
icantly  longer  than  the  two-dimensional  counterpart,  indicating 
that  transport  processes  along  the  mixture  fraction  iso-contours 
must  be  important  and  that  autoignition  is  less  important  to  the 
stabilization  mechanism. 

5.  Stabilization  mechanism 


To  validate  this  formulation  in  the  current  configuration,  the  dissi¬ 
pation  rates  along  different  mixture  fraction  iso-contours  were 
sampled  from  the  two-dimensional  computations,  normalized 
using  Eq.  (4),  and  compared  with  the  sampling  along  the  Zst 
iso-contour.  As  shown  in  Fig.  9,  the  normalized  dissipation  rates 
at  different  mixture  fractions  all  collapse  to  the  value  at  Zst. 
Therefore,  only  the  dissipation  rate  samplings  along  the  Zst 
iso-contour  were  needed  to  perform  the  unsteady  flamelet 
calculation. 

To  account  for  the  differential  diffusion,  species  Lewis  numbers 
for  LFA  were  specified  the  same  as  in  the  two-dimensional  compu¬ 
tations.  The  governing  equations  for  species  and  temperature  fol¬ 
low  Eqs.  (24)  and  (25)  in  Pitsch  and  Peters  [19]. 

In  the  current  work,  the  time  history  of  the  dissipation  rate  /st 
was  specified  in  LFA  according  to  the  two-dimensional  computa¬ 
tion.  To  avoid  the  ill-defined  Lagrangian  time  in  the  recirculation 
zone,  time  zero  was  defined  at  a  downstream  location  ten  times 
the  thickness  of  the  wall.  Accordingly,  the  species  and  temperature 
profiles  along  the  radial  cut  at  this  location  were  specified  as  the 
initial  conditions  for  the  flamelet.  Based  on  these  initial  conditions 
and  Xst  time  history  profiles,  the  unsteady  flamelets  were 


With  the  above  analysis  based  on  species  profiles,  Chemical 
Explosive  Mode  Analysis,  and  Lagrangian  Flamelet  Analysis,  the 
transition  of  the  stabilization  mechanism  and  the  coupling 
between  autoignition  chemistry  and  flame  propagation  can  be 
clearly  identified.  In  the  current  study,  two  fundamental  stabiliza¬ 
tion  mechanisms  are  relevant:  the  kinetic  stabilization  mechanism, 
due  to  the  balance  between  the  autoignition  delay  time  and  flow 
residence  time,  and  the  kinematic  stabilization  mechanism,  due 
to  the  balance  between  the  local  premixed  flame  propagation 
velocity  and  the  local  flow  velocity. 

In  this  stratified  composition  and  temperature  field,  autoigni¬ 
tion  and  flame  propagation  are  coupled  through  thermal  and  rad¬ 
ical  interactions,  for  the  accumulation  of  the  upstream  radicals  and 
heat  release  from  autoignition  accelerate  the  flame  propagation 
velocity.  The  flame  also  transfers  heat  and  radicals  through  back 
diffusion  processes  to  the  upstream,  which  could  also  facilitate 
autoignition. 

In  the  current  study,  the  stabilization  mechanism  was  deter¬ 
mined  by  comparing  the  two-dimensional  computations  with  the 
one-dimensional  inhomogeneous  autoignition  predicted  by  LFA. 
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Fig.  11.  Extended  regimes  of  the  stabilization  mechanism  as  the  coflow  boundary  temperature  increases. 


When  these  two  time  history  profiles  agree  well,  the  case  is  char¬ 
acterized  as  kinetically  stabilized.  Specifically,  at  700  K,  the 
one-dimensional  LFA  agrees  very  well  with  the  two-dimensional 
computation  at  all  of  the  mixture  fractions  examined.  Therefore, 
the  700  K  case  is  characterized  as  kinetically  stabilized.  As  the 
boundary  temperature  increases,  the  influences  from  premixed 
flame  propagation  become  important,  as  predictions  by  LFA  lag 
behind  the  CFD  results  for  some  mixture  fractions.  At  800  K,  an 
autoignition  front  stabilizes  the  multibrachial  structure  at  rich 
mixture  fractions,  due  to  the  shorter  ignition  delay  time  resulting 
from  the  NTC  chemistry,  and  a  modified  triple  flame  structure  sta¬ 
bilizes  slightly  downstream  of  this  front  at  leaner  mixture  frac¬ 
tions,  as  shown  in  Fig.  2.  Further  increasing  the  boundary 
temperature  results  in  higher  flame  propagation  velocities. 
Therefore,  the  flame  front  at  leaner  mixture  fraction  depends  less 
on  radical  accumulation  ahead  of  the  flame  and  propagates 


upstream,  and  the  stabilization  is  influenced  by  both  inhomoge¬ 
neous  autoignition  and  premixed  flame  propagation.  The  transition 
to  a  kinematically  stabilized  flame  structure  is  achieved  for  the 
1100  K  case,  where  the  local  flame  propagation  velocity  balances 
the  incoming  flow  velocity.  Therefore,  the  flame  structure  stabi¬ 
lizes  close  to  the  nozzle  exit  and  depends  least  on  radical  accumu¬ 
lation  from  upstream.  Consequently,  the  one-dimensional  LFA 
predictions  departure  from  the  two-dimensional  computation 
significantly. 

Based  on  the  understanding  obtained  from  the  current  study, 
further  extension  of  the  stabilization  regime  can  be  made,  as 
shown  in  Fig.  11.  For  fixed  inlet  flow  velocity,  when  the  boundary 
temperature  is  sufficiently  low,  the  mixture  cannot  be  autoignited, 
and  it  is  essentially  a  frozen  flow.  Even  when  an  external  ignition 
source  is  applied,  the  flame  cannot  keep  up  with  the  excessive  high 
flow  velocity,  such  that  the  flame  blows  out.  When  the  boundary 
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temperature  is  high  enough  to  activate  autoignition,  it  occurs  far 
downstream,  but  the  flame  propagation  velocity  still  cannot  keep 
up  with  the  flow  velocity.  As  a  consequence,  a  pure  kinetically  sta¬ 
bilized  autoignition  front  can  be  achieved,  which  is  similar  to  the 
700  K  case.  Conversely,  when  the  boundary  temperature  is  suffi¬ 
ciently  high,  the  flame  stabilizes  close  to  the  inlet,  where  the 
upstream  can  be  treated  as  frozen,  due  to  reduced  residence  time, 
which  is  similar  to  the  1100  K  case.  Therefore,  a  kinematically  sta¬ 
bilized  classical  triple  flame  structure  is  achieved.  Further  increase 
in  the  boundary  temperature  results  in  an  attached  flame  with  the 
increased  flame  speed.  Although  not  included  in  the  current  paper, 
an  attached  flame  was  computed  at  1 500  K.  In  between  the  kinet¬ 
ically  and  kinematically  stabilized  regimes,  there  is  a  transitional 
regime  governed  by  both  mechanisms,  which  corresponds  to  the 
800  and  900  K  cases.  Due  to  the  NTC  behavior  of  the  autoignition 
chemistry,  the  stabilization  point,  in  terms  of  mixture  fraction 
space,  varies,  and  the  complex  multibrachial  flame  structure 
appears. 

6.  Conclusions 

In  the  present  study,  two-dimensional  nonpremixed  DME 
flames  in  heated  air  coflows  were  computed.  The  computations 
were  conducted  at  30  atmospheres  to  observe  the  influence  of 
NTC  chemistry  on  the  stabilization  mechanism.  A  uniform  and 
fixed  inlet  boundary  velocity  was  specified,  and  four  coflow  tem¬ 
perature  (700,  800,  900,  and  1100K)  cases  were  studied. 

The  heat  release  rate  profile  and  characteristic  species  profiles 
for  low  and  high  temperature  chemistry,  autoignition,  and  pre¬ 
mixed  flame  propagation  were  examined.  Further  investigation 
based  on  Chemical  Explosive  Mode  Analysis  and  Lagrangian 
Flamelet  Analysis  enabled  the  determination  of  the  evolution  of 
the  controlling  chemical  pathways  and  the  stabilization 
mechanism. 

The  700  K  case  was  characterized  as  kinetically  stabilized,  for, 
neglecting  the  diffusion  processes  along  mixture  fraction 
iso-contours,  the  one-dimensional  LFA  agrees  with  the 
two-dimensional  CFD  responses.  As  the  boundary  temperature 
increases,  the  leading  point  of  the  heat  release  profile  shifts  to 
richer  mixture  fractions  and  then  shifts  back  due  to  the  NTC  effect 
on  the  autoignition  process  and  the  coupling  between  autoignition 
and  premixed  flame  propagation  chemistry.  Stabilization  is  also 
affected  by  both  inhomogeneous  autoignition  and  premixed  flame 
propagation,  as  in  the  800  and  900  K  cases.  The  1100  K  case  was 
characterized  as  kinematically  stabilized,  for  it  exhibits  the  classical 
triple  flame  structure,  with  stabilization  achieved  due  to  the  bal¬ 
ance  between  the  premixed  flame  propagation  velocity  and  the 
local  incoming  flow  velocity. 

Based  on  the  present  study,  the  following  extended  stabilization 
regimes  are  identified.  For  sufficiently  high  inlet  velocity,  as  the 
boundary  temperature  increases  from  the  cold  case,  the  flow  is 
chemically  frozen,  and  the  mixture  is  nonautoignitive;  even  the 
flame  generated  by  an  external  ignition  source  will  blow  out. 
When  the  mixture  can  be  autoignited,  the  kinetically  stabilized 
autoignition  front  gradually  transits  to  a  kinematically  stabilized 
classical  triple  flame,  where  the  premixed  flame  front  propagation 
velocity  balances  the  local  incoming  flow  velocity.  The  triple  flame 
will  eventually  become  attached  when  the  boundary  temperature 
is  sufficiently  high,  and  the  flame  speed  is  sufficiently  fast. 

Further  study  on  the  effects  of  fuel  dilution  and  inlet  velocity  on 
the  lift-off  height,  stabilized  flame  structure,  and  stabilization 
mechanism,  is  suggested.  For  example,  the  boundary  velocity  can 
be  varied  to  change  the  flow  residence  time,  and  dilution  can  be 


added  to  the  fuel  stream  to  change  the  chemical  time  scale. 
Moreover,  if  autoignition  is  the  dominant  stabilization  mechanism 
and  the  scalar  dissipation  rate  sufficiently  low,  the  nonmonotonic 
lifted  height  variation  can  be  observed  as  the  boundary  tempera¬ 
ture  changes,  which  is  plausible  based  on  homogeneous  autoigni¬ 
tion  and  nonpremixed  counterflow  observations  [17].  Finally, 
high  pressure  experiments  are  needed  to  provide  potential  exper¬ 
imental  substantiation  of  these  predicted  multibrachial  structures. 
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The  structure  and  stabilization  mechanism  of  laminar  nonpremixed  autoignitive  DME/air  coflow  flames  was 
investigated  at  elevated  temperatures  and  pressures.  Computations  with  detailed  chemistry  were  performed 
for  DME  and  heated  coflow  air  at  30  atm  with  uniform  inlet  velocities  (2.4, 3.2,  and  8.0  m/s)  imposed  for  both 
streams.  The  heat  release  rate  profiles  are  first  examined  for  each  case  to  demonstrate  a  multibrachial  ther¬ 
mal  structure.  Species  concentrations  and  temperature  were  sampled  along  mixture  fraction  iso-contours, 
and  Chemical  Explosive  Mode  Analysis  (CEMA)  was  performed  to  identify  the  controlling  chemistry  at  rep¬ 
resentative  points.  One-dimensional  Lagrangian  Flamelet  Analysis  (LFA)  was  also  performed  and  compared 
with  the  two-dimensional  computations  to  elucidate  the  relative  importance  of  diffusion  processes  parallel 
and  normal  to  the  mixture  fraction  gradient.  Various  coflow  temperatures  with  different  inlet  velocities  are 
examined  to  elucidate  their  influences  on  the  multibrachial  structure  as  well  as  the  stabilization  mechanism. 
NTC  (negative  temperature  coefficient)-affected  inhomogeneous  autoignition  and  the  coupled  effects  with 
premixed  flame  propagation  on  stabilization  are  further  studied.  It  is  found  that,  at  high  coflow  boundary 
temperatures  or  low  inlet  velocities,  the  classical  tribrachial  flame  structure  is  achieved,  and  autoignition 
contributes  less  to  the  stabilization  due  to  reduced  heat  and  radical  accumulation.  The  kinematic  balance 
between  the  local  flow  speed  and  flame  propagation  speed  is  the  dominant  stabilization  mechanism.  On  the 
contrary,  kinetic  stabilization  is  achieved  at  lower  coflow  temperatures  or  higher  inlet  velocities  as  autoigni¬ 
tion  becomes  dominant.  Due  to  the  transition  of  the  dominant  chemical  pathways  during  autoignition,  the 
kinetically  stabilized  structure  is  usually  multibrachial.  The  transition  of  different  stabilization  mechanisms 
can  be  made  by  changing  either  the  boundary  velocity  or  temperature  of  the  coflow.  Based  on  these  results 
and  previous  work  (Deng  et  al.,  2015)  [12],  a  regime  diagram  is  constructed  that  identifies  the  possible  stabi¬ 
lization  regimes:  blow  out,  kinetically  stabilized,  autoignition-propagation-coupled  stabilized,  kinematically 
stabilized,  and  burner  stabilized. 

©  2015  The  Combustion  Institute.  Published  by  Elsevier  Inc.  All  rights  reserved. 


1.  Introduction 

Two-dimensional  tribrachial  structures  (also  known  as  triple 
flames)  [1]  are  observed  in  nonpremixed  laminar  lifted  flames 
at  nonautoignitive  conditions.  The  dynamic  balance  between  the 
local  flame  propagation  speed  and  the  incoming  flow  speed  at  the 
triple  point  is  generally  considered  as  the  stabilization  mechanism 
[2].  However,  practical  engines  operate  at  elevated  pressures  and 
temperatures.  As  a  consequence,  the  propensity  for  autoignition  is 
significantly  enhanced,  and,  therefore,  the  thermal  and  chemical 
structure  of  the  tribrachial  flame,  as  well  as  the  stabilization  mech¬ 
anism,  could  be  affected  by  the  autoignition  process.  For  example, 
experimentally,  Chung  and  co-workers  have  investigated  autoignited 
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lifted  propane/nitrogen  [3],  methane/hydrogen  [4],  and  other  neat 
Ct-C4  hydrocarbon  flames  [5],  and  compared  their  lift-off  heights 
with  homogeneous  autoignition  delay  time.  Furthermore,  the  au¬ 
toignition  process  of  most  large  hydrocarbons  under  practical  engine 
conditions  could  possibly  lie  in  the  negative  temperature  coefficient 
(NTC)  regime,  in  which  the  overall  ignition  delay  time  increases  as 
the  initial  temperature  increases.  The  NTC  phenomenon  is  relevant  to 
engine  knock  [6]  and  has  been  extensively  studied  in  homogeneous 
systems  [7].  As  Law  and  co-workers  [8-10]  recently  demonstrated, 
ignition  characteristics  in  a  nonpremixed  system  can  also  be  affected 
by  NTC  chemistry,  especially  at  elevated  pressures  and/or  extended 
residence  times.  These  computational  and  experimental  studies 
were  conducted  in  the  nonpremixed  counterflow  system  where  the 
residence  time  is  well  characterized.  When  the  flow  residence  time 
and  NTC  chemistry  timescales  are  comparable,  the  two  processes 
are  strongly  coupled,  resulting  in  modified  system  response,  such  as 
autoignition  behavior. 


DISTRIBUTION  A:  Distribution  approved  for  public  release. 


4472 


S.  Deng  et  al.  /  Combustion  and  Flame  162  (2015)  4471-4478 


To  investigate  autoignition  with  NTC  chemistry  effects  in  non- 
premixed  lifted  flame  stabilization,  Krisman  et  al.  [11]  recently  con¬ 
ducted  a  numerical  study  of  dimethyl  ether  (DME)/air  nonpremixed 
flames  at  40  atmospheres  and  elevated  air  coflow  temperatures 
(700-1500  K)  and  observed  multibrachial  thermal  structures.  The  au¬ 
toignition  response  in  the  two-dimensional  computation  was  com¬ 
pared  with  that  of  homogeneous  autoignition  under  the  same  ini¬ 
tial  conditions.  A  transport  budget  analysis  of  methoxymethylperoxy 
(CH30CH202)  and  hydroxyl  (OH)  radicals,  which  represent  the  low 
and  high  temperature  chemistry,  respectively,  was  performed  to  dif¬ 
ferentiate  deflagration  from  autoignition.  The  stabilization  points  of 
the  multibrachial  structure  was  determined  with  CH30CH202  and 
OH  radical  mass  fractions  for  low  and  high  temperature  autoignition 
chemistry,  respectively,  and  varied  as  boundary  temperature  and  ve¬ 
locity  changed. 

More  recently,  to  further  elucidate  the  chemical  structure  of  the 
multibrachial  structure  and  the  roles  of  autoignition  and  flame  chem¬ 
istry  in  the  stabilization  mechanism,  the  authors  [12]  performed  a 
numerical  study  of  nonpremixed  DME/air  coflow  flames.  Chemical 
Explosive  Mode  Analysis  (CEMA)  was  adopted  to  identify  locally 
dominant  reactions,  and  Lagrangian  Flamelet  Analysis  (LFA)  was 
adopted  to  identify  the  dominant  combustion  mode.  The  comparison 
with  the  two-dimensional  computation  was  able  to  quantify  the 
relative  importance  of  transport  processes  parallel  and  normal  to  the 
mixture  fraction  gradient  and  elucidate  the  dominant  stabilization 
mechanism.  For  increasing  coflow  boundary  temperature  at  con¬ 
stant  inlet  velocities,  the  stabilization  mechanism  transitioned  from 
kinetic  to  kinematic  stabilization. 

In  the  present  study,  nonpremixed  DME/air  coflow  flames  at  el¬ 
evated  temperatures  and  pressures  were  further  studied.  The  objec¬ 
tive  of  the  current  work  is  fourfold:  first,  to  elucidate  transport  effects 
on  stabilization,  parallel  to  the  previous  work  [12],  which  focused  on 
chemical  effects;  second,  to  demonstrate  the  effects  of  NTC  chem¬ 
istry  on  the  multibrachial  structure;  third,  to  understand  the  transi¬ 
tion  between  kinetic  and  kinematic  stabilization  mechanism  and  the 
coupling  effects;  and,  fourth,  to  structure  a  complete  regime  diagram 
that  includes  both  chemical  and  transport  effects. 

As  a  final  note,  practical  engine  conditions  are  highly  turbulent, 
and  the  autoignition  phenomenon  depends  on  both  chemistry  and 
turbulent  mixing.  For  example,  in  a  DNS  study,  Yoo  et  al.  [13]  ob¬ 
served  the  cyclic  movement  of  the  stabilization  point  of  a  turbulent 
lifted  ethylene  jet  flame  in  highly-heated  coflow,  which  is  a  conse¬ 
quence  of  consecutive  autoignition  events  in  the  high-speed  jet  and 
coflow.  When  fuels  involve  more  complicated  chemical  characteris¬ 
tics,  such  as  NTC  effects,  turbulence  plays  different  and  potentially 
competing  roles.  As  demonstrated  in  a  more  recent  computational 
work  by  Echekki  and  Ahmed  [14],  although  scalar  dissipation  rate 
tends  to  delay  ignition  due  to  heat  and  radical  loses  from  nascent 
kernels,  enhanced  mixing  ensures  much  larger  volumetric  heat  re¬ 
lease  rate  after  ignition.  However,  neither  of  these  works  has  consid¬ 
ered  high  pressure  regimes  and  analyzed  the  complicating  effects  of 
NTC  in  detail.  Therefore,  in  order  to  better  understand  flame  stabiliza¬ 
tion  and  provide  insights  for  future  studies  on  turbulent  lifted  flames 
at  elevated  temperatures  and  pressures,  the  current  work  focuses  on 
laminar  conditions. 

2.  Computational  details 

An  axisymmetric  DME  stream  at  300  K  is  surrounded  by  heated 
coflow  air  at  30  atmospheres.  The  fuel  nozzle  diameter  D  is  0.8  mm, 
and  the  fuel  and  air  are  initially  separated  by  an  adiabatic,  no-slip 
wall  with  thickness  D/20.  The  diameter  of  the  coflow  is  3.9  mm  with 
adiabatic,  slip  wall  boundary  conditions.  This  diameter  was  chosen 
to  be  wide  enough  such  that  further  widening  of  the  domain  did  not 
influence  the  computational  results.  Uniform  inlet  velocities  of  2.4, 
3.2,  and  8.0  m/s  for  both  streams  were  specified.  The  outlet  boundary 


Table  1 

Computational  domain  and  number  of  grid  points. 


Inlet  velocity  [m/s] 

2.4 

3.2 

8.0 

Lx  [mm] 

3.5 

3.5 

15 

Nx 

1536 

1536 

3072 

Nr 

176 

176 

176 

condition  is  a  convective  outflow,  which  is  a  Neumann  condition  at 
steady-state. 

The  governing  equations,  transport  model,  and  chemical  model 
were  adopted  to  be  the  same  as  in  Deng  et  al.  [12].  In  brief,  the  Navier- 
Stokes  equation  with  buoyancy  effects  in  the  streamwise  direction 
and  the  conservation  equations  of  mass,  species,  and  energy  were 
solved.  The  species  diffusivities  are  determined  assuming  a  constant, 
nonunity  Lewis  number  and  kept  the  same  as  in  the  previous  work 
[12].  The  conserved  scalar  mixture  fraction  Z  is  specified  as  unity  and 
zero  for  the  fuel  stream  and  coflow,  respectively,  and  is  computed  by 
solving  a  conserved  scalar  transport  equation  with  unity  Lewis  num¬ 
ber  [15].  A  DME  skeletal  mechanism  of  39  species  [16],  which  was 
reduced  from  the  well-validated  detailed  mechanism  of  Zhao  et  al. 
[17],  was  adopted  as  the  chemical  model. 

The  governing  equations  with  the  low-Mach  number  formulation 
are  solved  using  NGA  [18].  The  momentum  and  scalar  equations  are 
discretized  with  a  second-order  centered  scheme  and  a  third-order 
WENO  scheme  [19],  respectively,  on  a  staggered  mesh.  An  iterative 
second-order  semi-implicit  Crank-Nicolson  scheme  is  adopted  for 
temporal  integration  [20].  The  chemical  source  terms  for  the  species 
and  energy  equations  are  integrated  using  the  CVODE  package  [21  ]. 

Leveraging  previous  grid  convergence  studies  [12],  uniform  grid 
spacing  in  the  axial  direction  was  set  to  Ax  =  2. 2-4. 8  |xm,  depending 
on  the  case.  Nonuniform  grid  spacing  in  the  radial  direction  was  set  to 
minimum  A r  =  2.5  |xm  to  resolve  the  mixing  layer  near  the  thin  wall, 
and  the  grid  stretch  rate  was  less  than  3%.  Details  about  the  numerical 
discretization  are  summarized  in  Table  1. 

3.  Transport  effects 

Transport  effects  on  nonpremixed  coflow  flame  stabilization  were 
demonstrated  by  fixing  the  coflow  temperature  at  900  K  while  vary¬ 
ing  the  uniform  inlet  velocities  as  2.4,  3.2,  and  8.0  m/s.  The  3.2  m/s 
case  was  computed  in  previous  work  [12]. 

3.1  Thermal  and  chemical  structure 

The  heat  release  rate  profiles  for  the  three  cases  are  shown  in 
Fig.  1.  A  qualitative  determination  of  the  stabilization  point  is  the 
most  upstream  point  on  the  largest  heat  release  contour  (the  lead¬ 
ing  point),  colored  by  red.  The  mixture  fraction  iso-contours  of  Zst  = 
0.1005,  Z  =  0.2,  and  Z  =  0.3  are  delineated  in  solid  black  lines,  from 
right  to  left. 

When  the  inlet  velocity  is  the  lowest,  2.4  m/s,  a  tribrachial 
thermal  structure  is  observed  very  similar  to  that  of  the  classical 
triple  flame.  The  triple  point  at  Z  =  0.15,  where  the  three  large  heat 
release  branches  intersect,  is  also  the  stabilization  point.  Some  heat 
release  can  be  found  upstream  of  the  tribrachial  thermal  structure  for 
the  partially  reacting  mixture  at  elevated  temperature  but  is  much 
less  than  the  heat  release  from  the  flame  structure.  As  the  inlet  ve¬ 
locity  increases  to  3.2  m/s,  another  branch  with  large  heat  release  is 
found  attached  to  the  tribrachial  structure  around  Z  =  0.2.  The  stabi¬ 
lization  point  is,  again,  the  same  as  the  triple  point.  This  structure  has 
been  analyzed  in  our  previous  work  [12].  However,  as  the  inlet  ve¬ 
locity  further  increases  to  8.0  m/s,  the  stabilization  point  is  no  longer 
on  the  tribrachial  structure.  Instead,  it  is  found  to  be  near  Z  =  0.25 
and  is  the  intersection  point  of  two  trailing  heat  release  branches. 
Attached  to  the  leaner  branch,  there  is  a  tribrachial  structure  that 
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Fig.  1.  Heat  release  rate  [J/m3  s]  profiles.  The  iso-contours  of  Zst,  Z  =  0.2,  and  Z  =  0.3  are  outlined  from  right  to  left  in  solid  lines,  respectively.  The  CEMA  sampling  points  for  2.4 
and  8.0  m/s  cases  are  marked  along  the  iso-contours.  (For  interpretation  of  the  references  to  color  in  this  figure,  the  reader  is  referred  to  the  web  version  of  this  article.) 


appears  similar  to  the  triple  flame  structure.  This  multibrachial 
structure  is  similar  to  the  one  found  in  previous  work  at  a  higher 
oxidizer  temperature  (800  K)  but  lower  velocity  (3.2  m/s)  [12]. 

The  controlling  chemistry  of  the  three  cases  was  studied  with 
CEMA  [22,23].  Briefly,  local  species  concentrations  and  temperature 
are  sampled  from  the  two-dimensional  computation  and  input  into 
CEMA  to  evaluate  the  eigenvalues  of  the  Jacobian  matrix  (I(  +  1  by 
K+  1,  where  I<  is  the  number  of  species)  of  the  chemical  source 
terms.  The  eigenmode  associated  with  the  eigenvalue  Ze,  which  has 
the  largest  real  part  among  all  the  eigenvalues,  is  defined  as  a  chemi¬ 
cal  explosive  mode,  if 

Rc( A-e)  >  0,  Xe  =  beetle,  (1) 


where  be  and  ae  are  the  left  and  right  eigenvectors,  respectively,  asso¬ 
ciated  with  A.e,  and  Jw  is  the  chemical  Jacobian  matrix.  The  existence 
of  the  chemical  explosive  mode  indicates  the  propensity  of  local  mix¬ 
ture  autoignition  given  an  isolated,  adiabatic,  and  constant  volume 
environment.  Furthermore,  the  detailed  chemical  reactions  that  con¬ 
tribute  to  the  chemical  explosive  mode  can  be  quantified  by  the  ex¬ 
plosion  participation  index: 


PI  = 


|(be-S)®R| 
sum(\(be  -S)  0  R|)  ’ 


(2) 


where  S  is  the  stoichiometric  coefficient  matrix,  R  is  the  vector  of  the 
net  rates  for  the  reactions,  and  0  denotes  elementwise  multiplication 
of  two  vectors. 

In  the  present  study,  such  samplings  were  conducted  along  Zst, 
Z  =  0.2,  and  Z  =  0.3  iso-contours,  as  shown  in  Fig.  1.  Based  on  the 
explosive  mode  and  participation  index,  the  evolution  of  the  domi¬ 
nant  reactions  is  shown  in  Fig.  2. 

At  2.4  m/s,  the  dominant  reactions  along  Zst  and  Z  =  0.2  iso¬ 
contours  evolve  in  similar  ways:  upstream  of  the  tribrachial  struc¬ 
ture  (points  A,  B,  and  D),  low  temperature  chemistry,  characterized 
by  reactions  involving  CH30CH202  radicals,  is  active.  Due  to  the  high 
diffusivity  of  H  radicals  and  the  elevated  pressure,  the  H  radical  re¬ 
combination  reaction  (H  +  02  +  M  <<=>>  H02  +  M)  is  important.  At  the 
most  reactive  region  (points  C  and  E),  the  H  radical  branching  reac¬ 
tion  (H  +  02  O  +  OH)  becomes  the  most  important  chain  branch¬ 
ing  reaction,  indicating  the  transition  to  high  temperature  chemistry. 
On  the  contrary,  for  the  8.0  m/s  case,  while  low  temperature  chem¬ 
istry  is  still  active  upstream  of  the  multibrachial  structure,  the  dom¬ 
inant  chain  branching  reaction  is  the  hydrogen  peroxide  branching 
reaction  (H202  +  M  OH  +  OH  +  M).  Moreover,  the  dominant  reac¬ 
tions  along  the  Zst  and  Z  =  0.2  iso-contours  evolve  in  different  ways. 
Along  the  Z  =  0.2  iso-contour,  from  point  D  to  E,  the  hydrogen  per¬ 
oxide  branching  reaction  is  always  the  dominant  reaction,  indicating 
the  role  of  low-to-intermediate  temperature  autoignition  chemistry 


[24].  Although  this  is  the  case  at  point  A  on  the  Zst  iso-contour,  the 
H  radical  chain  branching  reaction  becomes  dominant  at  point  C,  the 
most  reactive  zone,  indicating  that  the  dominant  chemical  pathway 
shifts  to  high  temperature  chemistry.  CEMA  results  for  the  3.2  m/s 
case  show  similar  transitions  as  those  of  the  8.0  m/s  case,  although 
their  thermal  structures  appear  different.  Therefore,  further  compu¬ 
tational  diagnostics  is  needed  to  identify  the  dominant  combustion 
mode. 


3.2.  Stabilization  mechanism 


The  above  CEMA  results  have  demonstrated  the  importance  of  au¬ 
toignition  chemistry  in  the  multibrachial  structure.  However,  further 
analysis  is  still  needed  to  elucidate  the  role  that  autoignition  plays  in 
the  stabilization  mechanism.  As  conducted  previously  [12],  LFA  [25] 
utilizes  the  initial  conditions  provided  by  the  two-dimensional  com¬ 
putation  to  provide  the  time  history  of  the  one-dimensional  inhomo¬ 
geneous  autoignition.  As  only  the  diffusion  processes  parallel  to  the 
mixture  fraction  gradient  direction  are  accounted  for  in  this  analy¬ 
sis,  the  comparison  of  this  one-dimensional  flamelet  and  the  two- 
dimensional  result  gives  the  relative  importance  of  transport  parallel 
and  normal  to  the  mixture  fraction  gradient  and  thus  the  relative  im¬ 
portance  of  inhomogeneous  autoignition  and  premixed  flame  propa¬ 
gation  to  the  stabilization  mechanism. 

Specifically,  the  unsteady  flamelets  were  computed  with 
FlameMaster  [26].  Species  mass  fractions  and  temperature  were 
sampled  at  ten  times  the  wall  thickness  downstream  of  the  inlet 
to  avoid  the  influence  from  the  recirculation  zone  (behind  the  thin 
wall)  on  the  initial  condition  in  LFA.  The  time  history  of  the  scalar 
dissipation  rate  was  sampled  along  the  Zst  iso-contour  from  the  two- 
dimensional  computation.  The  flamelet  time  was  computed  from  the 
two-dimensional  computational  results,  along  the  Zst  iso-contour: 

t  =  L  (u  +  uz)(x0l(Z  =  Zst)dX'-  (3) 


The  axial  component  of  the  mixture  fraction  iso-contour  propaga¬ 
tion  speed  relative  to  the  fluid  convection  uz  is  taken  into  considera¬ 
tion  in  addition  to  the  axial  component  of  the  fluid  convection  veloc¬ 
ity  u.  The  expression  for  uz  is  adopted  from  Lignell  et  al.  [27]: 


uz  =  - 


V.(pDzVZ) 

P|VZ| 


(4) 


where  Dz  is  the  mixture  fraction  diffusivity,  with  unity  Lewis  number, 
and  p  is  the  density.  The  normal  vector  n,  defined  as 


VZ 

n-  ivzr 


(5) 
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Fig.  2.  Normalized  participation  index  at  2.4  and  8.0  m/s.  Sampled  locations  are  delineated  in  Fig.  1. 
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Fig.  3.  Comparison  between  CFD  and  LFA  results. 


indicates  the  direction  of  this  diffusion  induced  relative  velocity.  Only 
the  scalar  dissipation  rate  at  stoichiometric  mixture  fraction  is  sam¬ 
pled  from  the  two-dimensional  computations;  the  dissipation  rates 
at  other  mixture  fractions  were  computed  assuming  the  following 
form  [28]: 


X(Z)  =  x(Zst) 


exp(-2[erfc  1  (2Z)  ]2) 
exp  ( -  2[erfc-1  (2Zst)]2) 


X(Zst)/(Z;Zst). 


(6) 


Validation  of  this  expression  was  provided  in  previous  work  [12]. 

The  Lagrangian  time  history  profiles  of  the  two-dimensional  com¬ 
putation  and  one-dimensional  LFA  are  shown  in  Fig.  3.  For  each  in¬ 
let  velocity  case,  the  temperature  profiles  are  compared  along  Zst, 
Z  =  0.2,  and  Z  =  0.3. 

For  the  2.4  m/s  case,  LFA  fails  to  match  the  two-dimensional  re¬ 
sult  at  all  three  mixture  fractions,  indicating  that  transport  processes 
normal  to  the  mixture  fraction  gradient  are  crucial,  which  further  in¬ 
dicates  that  flame  propagation  is  the  dominant  stabilization  mech¬ 
anism.  At  3.2  m/s,  LFA  slightly  lags  behind  the  two-dimensional  re¬ 
sult  at  Zst  but  matches  well  at  Z  =  0.2  and  Z  =  0.3.  Recalling  the  heat 
release  profile  in  Fig.  1,  these  results  indicate  that  the  tetrabrachial 
structure  consists  of  a  tribrachial  structure,  at  which  flame  propaga¬ 
tion  is  not  negligible,  and  the  richer  branch  that  intersects  with  the 
tribrachial  flame  is  an  autoignition  front,  whose  response  is  well  cap¬ 
tured  with  the  one-dimensional  flamelet  model.  As  a  result,  stabiliza¬ 
tion  of  the  3.2  m/s  case  is  characterized  as  a  mixed  mode  of  inhomo¬ 
geneous  autoignition  and  premixed  flame  propagation,  depending  on 
the  local  mixture  fraction.  At  8.0  m/s,  LFA  agrees  well  with  the  two- 
dimensional  result  at  all  mixture  fractions,  indicating  that  the  trans¬ 
port  processes  normal  to  the  mixture  fraction  gradient  are  negligible. 


Therefore,  the  stabilization  mechanism  is  characterized  as  inhomo¬ 
geneous  autoignition. 

3.3.  Autoignition  and  flame  interaction 

As  shown  by  LFA,  under  some  conditions,  inhomogeneous  au¬ 
toignition  and  premixed  flame  propagation  both  contribute  to  flame 
stabilization,  resulting  in  a  multimode  stabilized  regime.  Further¬ 
more,  the  interaction  between  autoignition  and  flame  propagation 
is  tightly  coupled.  If  the  thermal  structure  is  mainly  kinetically  sta¬ 
bilized,  heat  and  radicals  generated  by  autoignition  will  modify  the 
downstream  thermal  and  chemical  environment  and  thus  the  local 
flame  speed.  On  the  contrary,  if  stabilization  is  mainly  kinematic  in 
nature  then  heat  and  radicals  generated  by  the  flame  can  back  dif¬ 
fuse  upstream,  modifying  the  reactivity  upstream. 

To  demonstrate  these  complex  interactions  and  understand  the 
transition  between  kinetic  and  kinematic  stabilization  mechanisms, 
the  LFA  results  for  the  2.4  m/s  case  were  further  analyzed.  As  shown 
in  Fig.  3,  if  there  was  a  kinetically  stabilized  inhomogeneous  autoigni¬ 
tion  front,  this  front  would  stabilize  further  downstream  than  the 
kinematically  stabilized  flame  front.  Although  not  shown,  CEMA  of 
these  LFA  solutions  show  the  same  evolution  of  the  controlling  chem¬ 
istry  as  the  8.0  m/s  case.  In  particular,  the  low-to-intermediate  tem¬ 
perature  hydrogen  peroxide  chain  branching  reaction  dominates  the 
transition  to  autoignition.  Therefore,  the  nature  and  the  qualitative 
structures  of  the  inhomogeneous  autoignition  fronts,  as  predicted  by 
LFA  in  Fig.  3  for  the  two  lower  inlet  velocity  cases,  are  essentially  the 
same  as  the  8.0  m/s  case. 

A  general  description  of  the  initiation  of  these  multibrachial  in¬ 
homogeneous  autoignition  fronts  in  all  three  cases  is  that,  due  to 
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radical  accumulation  and  heat  release,  the  controlling  chemistry 
shifts  from  low  temperature  chemistry,  represented  by  CH30CH202 
reactions,  to  hydrogen  peroxide  branching  reactions.  At  some  mixture 
fractions,  higher  temperatures  and  more  oxidizer  supply  enable  the 
dominant  chemistry  to  transition  further  to  high  temperature  chem¬ 
istry,  as  characterized  by  the  H  radical  branching  reaction.  As  shown 
in  Fig.  4,  there  are  double  or  triple  heat  release  peaks,  depending  on 
the  mixture  fraction.  However,  for  all  mixture  fractions,  these  peaks 
correlate  very  well  with  the  inflection  points  on  the  temperature  pro¬ 
files  and  peaks  on  the  OH  radical  profiles.  Comparing  them  with  the 
CH30CH202  radical  profile,  it  is  clear  that  low  temperature  chemistry 
results  in  the  first  peak  of  heat  release  profile  and  produces  OH  radi¬ 
cals.  Hydrogen  peroxide  accumulates  until  it  decomposes  at  the  sec¬ 
ond  heat  release  peak  and  produces  OH  radicals  through  H202  +  M 
OH  +  OH  +  M.  At  Zst  and  Z  =  0.2,  there  is  a  third  heat  release  and 
OH  radical  peak,  which  is  due  to  the  H  +  02  0  +  OH  reaction.  At 

Zst,  the  second  and  third  peaks  appear  to  be  much  closer  compared 
with  those  at  Z  =  0.2,  probably  due  to  higher  temperature  and  more 
abundant  oxidizer,  such  that  the  H  radical  chain  branching  reaction 
is  activated  earlier.  Due  to  the  lack  of  oxidizer,  the  third  peak  is  not 
observed  atZ  =  0.3. 

The  multibrachial  structure  of  the  inhomogeneous  autoignition 
front  is  due  to  the  variation  in  ignition  delay  time  at  different  mixture 
fractions.  As  shown  in  Fig.  1,  the  stabilization  point  is  determined 
based  on  the  threshold  value  of  heat  release  rate.  Therefore,  the 
ignition  delay  time  of  inhomogeneous  autoignition  can  be  deter¬ 
mined  accordingly,  which  corresponds  to  the  heat  release  rate  peak 
that  exceeds  1012  J/m3  s.  For  the  8.0  m/s  case  shown  in  Fig.  5,  the 
mixture  is  autoignited  first  at  Z  =  0.24,  which  corresponds  to  the 
stabilization  point  in  Fig.  1.  Although  the  first  stage  ignition  delay 
time,  defined  by  the  first  heat  release  peak  through  low  temperature 
chemistry,  increases  monotonically  as  mixture  fraction  increases 
(temperature  decreases),  the  overall  ignition  delay  time  reaches  the 
shortest  at  Z  =  0.24  due  to  the  compensation  of  larger  heat  release 
after  the  first  ignition  stage  [8],  as  shown  in  Fig.  4.  As  a  consequence, 
although  low  temperature  chemistry  is  not  the  dominant  chemical 
pathway  at  the  stabilization  point  for  the  kinetically  stabilized  case, 
it  still  influences  the  location  of  the  stabilization  point  through 
the  mixture  fraction  dependent  heat  and  radical  accumulation 
upstream. 

However,  although  the  cases  in  the  current  study  are  all  initi¬ 
ated  by  inhomogeneous  autoignition,  the  stabilization  of  the  final 
structure  depends  on  the  dominant  transport  processes,  which  are 
influenced  by  the  inlet  velocities.  At  8.0  m/s,  heat  and  radical  back 
diffusion  from  the  autoignition  front  to  upstream  is  not  able  to  keep 
up  with  convection;  therefore,  the  reacting  front  is  kinetically  stabi¬ 
lized.  At  3.2  m/s,  as  flow  convection  is  weaker,  diffusion  processes 
become  inherently  important  and  couple  with  chemical  reactions  to 
induce  a  flame  front  propagating  upstream.  However,  the  propaga¬ 
tion  speed  of  the  flame  varies  with  composition  and  temperature. 
As  a  consequence,  around  Zst,  where  higher  temperature  and  near- 
stoichiometric  mixture  composition  enable  higher  local  flame  speed, 
the  propagation  of  the  reacting  front  balances  the  incoming  flow 
velocity.  However,  such  a  balance  fails  at  richer  mixture  fractions 
where  kinetic  stabilization  dominates  due  to  enhanced  NTC-affected 
autoignition  at  richer  mixture  fractions  as  demonstrated  in  Fig.  5.  At 
2.4  m/s,  back  diffusion  is  important  at  all  mixture  fractions  such  that 
the  reacting  front  propagates  upstream  at  the  local  flame  speed,  as 
determined  by  the  local  composition  and  temperature.  Due  to  the 
increased  temperature  and  species  stratification  and  the  reduced 
thermal  and  radical  accumulation  from  autoignition,  the  propagation 
speed  of  this  reacting  front  is  less  influenced  by  inhomogeneous 
autoignition,  as  demonstrated  with  CEMA  and  LFA.  The  structure 
of  this  kinematically  stabilized  reacting  front,  which  is  generally 
tribrachial,  is  therefore  determined  by  the  variation  of  the  local  flame 
speed. 
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Fig.  4.  LFA  profiles  along  Zst,  Z  =  0.2,  and  Z  =  0.3  of  the  8.0  m/s  case. 


4.  Stabilization  regime  diagram 

The  above  sections  have  demonstrated  the  transport  effects  on 
the  thermal  and  chemical  structure  of  the  lifted  coflow  flames  as 
well  as  the  stabilization  mechanisms.  Combining  these  results  with 
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Fig.  6.  A  qualitative  regime  diagram  for  the  stabilization  mechanisms  as  the  boundary 
temperature  and  inlet  velocity  vary. 


the  chemical  effects  demonstrated  by  changing  the  coflow  boundary 
temperature  from  our  previous  study  [12],  a  two-dimensional  stabi¬ 
lization  regime  diagram  is  proposed,  as  shown  in  Fig.  6. 

Qualitatively,  when  the  boundary  temperature  is  not  high  enough 
to  activate  autoignition,  the  lifted  flame  appears  as  the  classical  triple 
flame  and  is  kinematically  stabilized.  When  the  inlet  velocity  is  below 
or  above  certain  threshold  values,  the  triple  flame  becomes  attached 
to  the  burner  or  is  blown  out,  respectively. 

As  the  boundary  temperature  is  elevated  enough  to  activate 
autoignition,  increasing  the  inlet  velocity  while  keeping  constant 
boundary  temperature,  the  flame  stabilization  mechanism  transits 
from  burner  stabilization  to  a  kinematic  balance  between  flame 
speed  and  incoming  flow  velocity,  then  to  multimode  stabilization 
influenced  by  both  flame  propagation  and  inhomogeneous  autoigni¬ 
tion,  and  finally  to  kinetic  stabilization  governed  by  inhomogeneous 
autoignition.  It  is  expected  that  the  crossover  velocities  between 
regimes  increase  with  increasing  boundary  temperature  because 
flame  speed  generally  increases  at  higher  temperature.  However,  it 
is  difficult  to  quantify  these  boundaries  as  local  composition  and 
temperature  vary  in  the  streamwise  direction,  and,  therefore,  the 
reference  flame  speed  cannot  be  calculated  based  on  the  upstream 
boundary  conditions.  Furthermore,  local  flame  front  curvature, 
cross-stream  species  stratification,  and  flow  divergence  approaching 
the  flame  front  also  modify  the  flame  speed.  As  a  consequence,  only 
a  qualitative  trend  is  demonstrated  in  Fig.  6. 

Similarly,  if  the  boundary  temperature  increases  at  fixed  inlet 
velocity,  transition  from  blow  out  to  burner  stabilized  regimes  is 


achieved  by  moving  horizontally  across  the  regime  diagram,  which 
was  discussed  in  our  previous  work  [12]. 

5.  Conclusions 

In  the  present  work,  axisymmetric  two-dimensional  laminar  non- 
premixed  DME  lifted  coflow  flames  at  elevated  temperatures  and 
pressures  were  computed.  Transport  effects  on  the  structure  and  sta¬ 
bilization  mechanism  were  demonstrated  by  changing  the  inlet  ve¬ 
locities  while  keeping  the  coflow  boundary  temperature  constant. 

The  heat  release  rate  profiles  were  examined  to  describe  the  ther¬ 
mal  structure,  and  CEMA  was  used  to  demonstrate  the  evolution  of 
the  controlling  chemistry.  Moreover,  one-dimensional  LFA  that  cap¬ 
tures  inhomogeneous  autoignition  was  adopted  to  identify  the  dom¬ 
inant  transport  directions  and  therefore  determine  the  dominant 
combustion  mode  and  stabilization  mechanism. 

At  2.4  m/s,  the  lifted  flame  appears  to  be  the  classical  triple 
flame  stabilized  by  the  balance  between  the  local  flame  speed  and 
incoming  flow  velocity,  and  it  is  therefore  characterized  as  kine¬ 
matically  stabilized.  As  the  inlet  velocity  increases,  such  a  balance 
cannot  be  achieved  at  certain  mixture  fractions.  Instead,  inhomo¬ 
geneous  autoignition  becomes  the  dominant  combustion  mode. 
As  a  consequence,  the  multibrachial  structure  is  stabilized  by  both 
premixed  flame  propagation  and  inhomogeneous  autoignition  and 
is  characterized  as  multimode  stabilized.  At  8.0  m/s,  the  kinematic 
balance  cannot  be  achieved  anywhere  in  the  flow  field  due  to  lack  of 
back  diffusion.  A  kinetically  stabilized  inhomogeneous  autoignition 
front  is  formed  where  diffusion  processes  along  the  mixture  fraction 
iso-contours  are  negligible  compared  to  the  gradient  direction.  In 
the  kinetically  stabilized  autoignition  front,  NTC  chemistry  plays 
an  important  role,  dictating  the  stabilization  point.  This  point  oc¬ 
curs  at  the  mixture  fraction  with  the  shortest  ignition  delay  time, 
which  is  a  compromise  of  the  first  stage  autoignition  delay  time 
and  heat  release  from  low  temperature  chemistry.  Conversely,  the 
kinematically  stabilized  flames  are  less  affected  by  NTC  chemistry, 
with  only  a  minor  effect  on  the  local  flame  speed  resulting  from 
the  upstream  accumulation  of  heat  and  radicals.  Combined  with  the 
extended  stabilization  regimes  demonstrated  in  our  previous  study 
[12],  an  extended  two-dimensional  stabilization  regime  diagram  was 
constructed,  considering  both  transport  (inlet  velocity)  and  chemical 
(coflow  boundary  temperature)  effects. 

This  work  represents  a  significant  contribution  to  the  understand¬ 
ing  of  nonpremixed  flames  at  autoignitive  conditions.  However,  prac¬ 
tical  engines  conditions  are  turbulent.  In  such  flows,  the  unsteady 
motion  of  eddies  and  enhanced  temperature  and  species  dissipation 
can  strongly  influence  stabilization.  Therefore,  further  investigation 
is  required  to  understand  the  role  of  turbulence  in  the  stabilization  of 
such  flames  and  determine  whether  the  stabilization  mechanism  is 
modified  or  not. 
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The  structure  and  dynamics  of  laminar  nonpremixed  dimethyl  ether  (DME)/air  coflow  flames  were  inves¬ 
tigated  at  elevated  temperatures  and  pressures.  Computations  with  detailed  chemistry  were  performed 
for  DME  and  heated  coflow  air  at  30  atm  with  uniform  but  sinusoidally  oscillating  inlet  velocities.  These 
unsteady  cases  were  compared  with  the  steady  results  from  Deng  et  al.  (2015)[14]  to  elucidate  the  effect 
of  oscillation  frequency  on  the  flame  dynamics.  To  benchmark  the  unsteady  cases,  a  normalized  displace¬ 
ment  velocity  was  defined  to  differentiate  flame  propagation  from  autoignition,  and  this  definition  was 
validated  against  the  steady  cases.  In  the  oscillating  reacting  flow,  transition  between  a  multibrachial 
autoignition  front  and  a  tribrachial  flame  occurs  periodically.  However,  unlike  the  harmonic  velocity  os¬ 
cillation,  the  combustion  mode  transition  is  hysteretic.  The  oscillation  cycle  starts  with  the  largest  inlet 
velocity,  with  the  multibrachial  thermal  structure,  located  downstream,  being  governed  by  autoignition 
chemistry.  As  flow  velocity  decreases,  the  autoignition  front  moves  upstream  and  transitions  to  a  tri¬ 
brachial  flame  near  the  lower  velocity  limit,  similar  to  the  steady  flow,  as  autoignition  chemistry  becomes 
weaker  with  decreasing  upstream  residence  time.  As  the  flow  velocity  increases  again,  the  tribrachial 
flame  is  convected  downstream,  and,  ultimately,  due  to  the  radical  and  heat  accumulation  in  time,  au¬ 
toignition  eventually  occurs  and  becomes  the  dominant  pathway.  The  finite  induction  time  for  autoigni¬ 
tion  results  in  the  hysteretic  behavior  during  the  decreasing-  and  increasing-velocity  cycles,  which  dimin¬ 
ishes  at  lower  oscillation  frequency  as  there  is  more  time  for  chemistry  to  respond  to  the  hydrodynamic 
changes  and  consequently  approach  steady  state.  At  the  relatively  low  oscillation  frequencies  investigated 
in  the  current  study,  first-stage  NTC  chemistry  is  less  affected  by  flow  dynamics  with  only  second-stage 
autoignition  and  flame  chemistry,  which  accounts  for  the  majority  heat  release,  coupled  with  the  flow 
oscillation. 
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1.  Introduction 

Fuel  injection  and  its  subsequent  mixing  and  reaction  with 
either  a  coflowing  air  stream  or  a  highly  turbulent  oxidizing 
environment  is  an  integral  process  in  the  operation  of  many 
combustors.  Due  to  experimental  and  computational  limitations, 
simplifications  are  usually  made  to  obtain  fundamental  under¬ 
standing  of  these  reacting  flows  that  are  then  extrapolated  to 
more  complex  conditions.  For  example,  laminar  nonpremixed 
coflow  flames  at  normal  ambient  temperatures  have  been  studied 
to  elucidate  the  coupling  between  fluid  dynamics  and  chemistry, 
leading  to  the  observation  that,  in  the  fuel  and  oxidizer  mixing 
layer,  a  two-dimensional  tribrachial  structure  (also  known  as  triple 
flame)  [1]  is  obtained.  Specifically,  both  lean  and  rich  premixed 
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flame  branches  and  a  trailing  diffusion  flame  tail  intersect  at  a 
triple  point.  Based  on  the  observation  of  such  tribrachial  laminar 
flames,  the  partially  premixed  flamelet  model  [2]  was  proposed  to 
explain  lifted  flames  in  nonpremixed  turbulent  jets  [3]. 

Previous  studies  have  indicated  that  nonpremixed  coflow  flame 
structure  can  be  modified  in  unsteady  flows.  For  example,  in  the 
experimental  investigation  of  Strawa  and  Cantwell  [4],  flow  in¬ 
stability  and  flame  breakup  was  achieved  by  imposing  a  small- 
amplitude,  periodic  velocity  fluctuation  to  nonpremixed  jet  flames 
at  elevated  pressures  and  low  Reynolds  numbers.  Later,  in  the 
computational  study  of  Sanchez-Sanz  et  al.  [5],  perturbation  fre¬ 
quency  effects  on  the  thermal  and  chemical  properties  of  the  flame 
in  such  periodically  time-varying  flows  were  evaluated.  Three 
regimes  were  found  depending  on  the  flame’s  Strouhal  number, 
S  =  D//2U,  with  D  and  /  denoting  the  fuel  jet  diameter  and  per¬ 
turbation  frequency,  respectively.  For  small  Strouhal  numbers  (S  = 
0.1),  perturbations  can  travel  far  downstream,  resulting  in  an  oscil¬ 
lating  flame.  Flame  surface  flickering  was  observed  when  S  «  0.2, 
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and  vigorous  flame  pinch-off  was  observed  at  5  =  0.5.  Larger  val¬ 
ues  of  S  confine  the  oscillation  to  the  jet’s  near-exit  region  with 
the  pulsation  having  minimal  effects  on  temperature  and  concen¬ 
tration  values.  The  unsteadiness  in  flickering  flames  also  increases 
pollutant  formation,  such  as  soot  [6]  and  carbon  monoxide  [7]. 
Mohammed  et  al.  [8]  followed  by  Dworkin  et  al.  [9]  conducted 
computational  and  experimental  studies  of  20  Hz  periodically- 
forced  methane/air  coflow  diffusion  flames.  Acetylene  production 
increased  [8]  and  the  oxidation  of  CO  to  C02  was  inhibited  [9]  in 
the  downstream  region  of  the  flame  at  certain  times  during  the 
flame’s  cyclic  history. 

Although  demonstrating  unsteady  effects  on  flow-chemistry 
coupling,  these  experimental  and  computational  investigations 
mainly  focused  on  simple  fuels,  such  as  methane,  and  were  lim¬ 
ited  to  nonautoignitive  conditions.  However,  as  demonstrated  by 
Chung  and  co-workers  through  a  series  of  experiments  on  lami¬ 
nar  lifted  C0  to  C8  flames  in  heated  coflow,  autoignition  can  be 
activated  even  at  atmospheric  pressure  [10,11].  Although  flames 
were  initiated  by  autoignition,  the  lifted  height  did  not  always 
correlate  well  with  the  ignition  delay  time.  Computationally,  Kris- 
man  et  al.  [12]  demonstrated  that,  under  more  realistic  engine  con¬ 
ditions  of  elevated  temperature  and  pressure  with  practical  fuels 
that  have  more  complex  chemical  kinetics,  the  transport-chemistry 
coupling  becomes  more  complicated.  Depending  on  the  condi¬ 
tions,  either  the  traditional  tribrachial  flame  or  autoignition  can 
be  dominant.  Their  findings  were  confirmed  and  further  discussed 
by  Deng  et  al.  [13,14]  through  a  series  of  computational  studies 
of  nonpremixed  dimethyl  ether  (DME)/air  coflow  flames  at  30  atm 
with  varying  inlet  velocities  and  coflow  temperatures,  recognizing 
that  DME  possesses  the  low-temperature  chemistry  (LTC),  negative 
temperature  coefficient  (NTC)  behavior.  A  regime  diagram  was  pro¬ 
posed  demonstrating  that  the  tribrachial  flame  is  favored  at  lower 
inlet  velocity  and  higher  coflow  temperature,  while  autoignition  is 
dominant  at  higher  inlet  velocity  and  relatively  lower  coflow  tem¬ 
perature. 

In  the  present  study,  unsteady  nonpremixed  DME/air  coflow 
flames  under  autoignitive  conditions  are  computationally  studied 
to  elucidate  the  coupling  between  unsteady  fluid  dynamics  and 
chemical  kinetics.  Various  oscillation  frequencies  were  imposed 
on  the  inlet  velocity,  with  the  maximum  and  minimum  velocities 
maintained  the  same  as  those  in  the  previous  steady  study  [14], 
which  correspond  to  an  autoignition  front  and  a  tribrachial  flame, 
respectively.  The  current  study  focuses  on  low  frequency  oscillation 
ranging  from  25  to  100  Hz,  which  covers  buoyancy-driven  insta¬ 
bility  frequencies  [8,9]  and  acoustic-driven  oscillation  frequencies 
in  gas  turbines  [15].  The  objective  of  the  current  study  is  three¬ 
fold.  The  first  objective  is  to  capture  the  transition  in  combus¬ 
tion  mode.  As  the  steady  cases  correspond  to  different  combustion 
modes,  it  is  expected  that,  at  certain  frequencies  of  velocity  os¬ 
cillation,  the  dominant  combustion  process  will  shift  between  the 
nonpremixed  tribrachial  flame  mode  and  the  autoignition  mode. 
The  second  objective  is  to  assess  the  thermal  and  chemical  differ¬ 
ences  during  such  transition  and  to  elucidate  the  transition  mech¬ 
anism.  The  third  objective  is  to  demonstrate  the  effects  of  oscil¬ 
lation  frequency  on  the  coupling  of  fluid  dynamics  and  chemical 
kinetics. 

2.  Computational  details 

The  geometry  in  this  work  is  the  same  as  that  in  Deng 
et  al.  [14].  Briefly,  axisymmetric  coflow  flames  at  30  atmospheres 
were  computed,  in  which  a  300  K  DME  stream  is  surrounded  by  a 
900  K  air  stream.  The  diameter  (D)  of  the  fuel  nozzle  is  0.8  mm, 
which  is  20  times  the  thickness  of  the  adiabatic,  no-slip  wall,  sepa¬ 
rating  fuel  and  coflow.  The  outer  diameter  of  the  coflow  is  3.9  mm. 
Adiabatic,  slip  wall  conditions  were  specified  at  the  outer  radial 


boundary.  The  same  inlet  velocities  were  imposed  for  both  streams 
and  are  uniform  in  space  and  sinusoidally  oscillating  in  time.  The 
maximum  (8.0  m/s)  and  minimum  (2.4  m/s)  velocities  were  set  to 
match  the  fastest  and  slowest  steady  cases  as  in  Deng  et  al.  [14]. 
Three  oscillation  frequencies  (25,  50,  and  100  Hz)  were  investi¬ 
gated,  with  the  maximum  Strouhal  number,  based  on  the  veloc¬ 
ity  and  jet  radius,  estimated  to  be  less  than  0.02  to  avoid  flame 
pinch-off.  The  domain  length  is  15  mm,  with  a  convective  outflow 
boundary  condition.  The  computational  results  were  not  affected 
by  further  widening  or  lengthening  of  the  domain.  Discretization 
of  the  domain  was  guided  by  previous  convergence  studies  [13], 
with  a  3072  (x)  by  176  (r)  grid.  Uniform  spacing  in  the  axial  direc¬ 
tion  was  set  to  Ax  =  4.8  |xm,  and  nonuniform  spacing  in  the  radial 
direction  was  set  to  minimum  Ar  =  2.5  |xm  to  resolve  the  mixing 
layer  near  the  thin  wall.  The  grid  stretch  rate  is  less  than  3%. 

The  Navier- Stokes  equations  with  buoyancy  effects  in  the 
streamwise  direction  and  the  conservation  equations  of  mass, 
species,  and  temperature  were  solved.  The  species  diffusivities 
were  determined  assuming  constant,  nonunity  Lewis  numbers  and 
kept  the  same  as  in  Deng  et  al.  [13].  The  conserved  scalar  mix¬ 
ture  fraction  (Z)  was  specified  as  unity  and  zero  for  the  fuel  stream 
and  coflow,  respectively,  and  was  computed  by  solving  a  conserved 
scalar  transport  equation  with  unity  Lewis  number  [16].  DME  was 
chosen  as  the  fuel,  for  it  is  one  of  the  simplest  fuels  that  has 
NTC  chemistry  [17],  as  noted  earlier.  A  skeletal  mechanism  of  39 
species  [18],  which  was  reduced  from  the  well-validated  detailed 
mechanism  of  Zhao  et  al.  [19],  was  adopted  as  the  chemical  model. 

A  low-Mach  number  code  NGA  [20]  was  adopted  to  solve 
the  discretized  governing  equations  on  a  staggered  mesh.  A 
second-order  centered  scheme  was  used  for  the  momentum  equa¬ 
tions,  while  a  third-order  WENO  scheme  [21]  was  used  for  the 
scalar  equations.  An  iterative  second-order  semi-implicit  midpoint 
scheme  was  adopted  for  temporal  integration  [22]  utilizing  Strang 
splitting  between  transport  and  chemistry  in  the  scalar  equa¬ 
tions  [23].  The  chemical  source  terms  for  the  species  and  temper¬ 
ature  equations  were  integrated  using  the  CVODE  package  [24]. 

3.  Results  and  discussion 

In  this  section,  the  thermal  structure  of  the  unsteady  non¬ 
premixed  DME/air  coflow  flame  is  first  described  and  compared 
qualitatively  with  the  previous  steady  study  by  Deng  et  al.  [14]. 
The  definition  of  normalized  displacement  velocity  is  then  intro¬ 
duced  and  validated  against  the  steady  cases  to  differentiate  the 
combustion  modes.  Finally,  oscillation  frequency  effects  on  the 
evolution  of  the  combustion  modes  are  analyzed. 

3.1.  Thermal  structure 

As  the  largest  and  smallest  inlet  velocity  cases  were  designed 
to  match  the  two  extreme  cases  in  Deng  et  al.  [14],  which  are 
of  different  thermal  structures,  it  is  expected  that  similar  thermal 
structures  will  be  obtained.  Furthermore,  these  thermal  structures 
might  transition  back  and  forth  in  response  to  the  oscillating  flow 
field.  Indeed,  such  transitions  were  observed  for  all  three  frequen¬ 
cies.  For  example,  the  evolution  of  the  thermal  structure  of  the 
100  Hz  case,  in  terms  of  the  heat  release  rate  profile,  is  demon¬ 
strated  in  Fig.  1.  The  oscillation  cycle  starts  with  the  largest  inlet 
velocity  of  8.0  m/s,  and  the  minimum  inlet  velocity  (2.4  m/s)  is 
achieved  at  half  cycle. 

At  8.0  m/s,  the  multibrachial  thermal  structure  is  located  fur¬ 
thest  downstream.  The  leading  point,  which  is  defined  as  the  most 
upstream  point  that  has  the  heat  release  rate  value  of  1012  W/m3, 
is  located  at  mixture  fraction  Z  =  0.24.  As  the  inlet  velocity  de¬ 
creases,  the  multibrachial  structure  moves  upstream,  without  ob¬ 
vious  change  of  the  leading  point  location,  in  terms  of  mixture 
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Fig.  1.  Heat  release  rate  [W/m3]  profile  evolution  during  one  oscillation  cycle  at  100  Hz.  The  iso-contours  of  Zst  =  0.1,  Z  =  0.2,  and  Z  =  0.3  are  outlined  from  right  to  left  in 
solid  lines,  respectively. 


fraction.  When  the  inlet  velocity  reaches  its  minimum,  the  multi- 
brachial  structure  transitions  to  a  tribrachial  structure,  and  the 
leading  point  switches  to  Z  =  0.14.  As  the  flow  velocity  increases, 
the  tribrachial  structure  is  pushed  downstream,  and  both  its  tri¬ 
brachial  shape  and  its  leading  point  mixture  fraction  remain  un¬ 
changed.  The  thermal  structure  returns  to  that  of  multibrachial 
when  the  flow  velocity  further  increases.  Such  transitions  in  struc¬ 
ture  repeat  once  a  new  oscillation  cycle  starts. 

3.2.  Differentiation  of  combustion  mode 


As  mentioned  in  the  Introduction,  Deng  et  al.  [14]  were  able  to 
relate  the  morphology  of  the  thermal  structures  to  two  different 
combustion  modes:  tribrachial  flame  and  autoignition.  Specifically, 
at  steady  state,  the  multibrachial  structure  in  the  8.0  m/s  case  is 
an  autoignition  front,  while  the  2.4  m/s  case  is  a  tribrachial  flame. 
In  this  previous  study,  species  mass  fraction  profiles  at  the  inlet 
of  the  two-dimensional  computation  were  treated  as  the  initial 
conditions  for  one-dimensional  Lagrangian  Flamelet  Analysis  [25], 
which  only  considers  diffusion  processes  parallel  to  the  mixture 
fraction  gradient  and  neglects  those  in  the  normal  direction.  When 
the  LFA  prediction  agreed  with  the  CFD  result,  transport  in  the  nor¬ 
mal  direction  of  the  mixture  fraction  gradient  was  negligible,  and 
autoignition  was  the  dominant  combustion  process.  However,  due 
to  the  unsteadiness  in  the  current  study,  such  comparison  between 
LFA  and  CFD  is  no  longer  applicable,  and  a  new  criterion  to  differ¬ 
entiate  the  modes  of  tribrachial  flame  and  autoignition  needs  to  be 
identified  and  validated  against  the  steady  cases. 

A  density-weighted  displacement  speed,  Sd,  is  often  used  to  dis¬ 
tinguish  between  deflagrations  and  spontaneous  ignition  fronts  in 
HCCI  combustion  [26],  which  is  defined  from  an  iso-line  of  species 
k  as  [27,28]: 


Sd 


1 

p7W\ 


dpYkVj,k\ 

dxj  )■ 


(1) 


where  Yk,  Vj  k,  and  cok  denote  species  mass  fraction,  diffusion  ve¬ 
locity  in  the  j-direction,  and  net  production  rate,  respectively,  and 
pu  is  the  density  of  the  unburnt  mixture.  The  choice  of  species 
k  and  its  iso-line  value  can  be  ambiguous.  Therefore,  major  prod¬ 
ucts,  such  as  C02,  H20,  H2,  CO,  and  combinations  of  these  products 
have  been  tested,  and  the  sampling  location  is  chosen  as  the  lead¬ 
ing  point,  as  defined  above,  to  enable  further  comparison  with  the 
steady  cases.  Sd  at  the  leading  point  is  insensitive  to  the  choice  of 
species  in  the  current  study,  for  less  than  5%  difference  was  ob¬ 
served  across  all  the  combinations.  Consequently,  H20  was  chosen 
for  simplicity.  Both  the  laminar  flame  speed  SL  and  the  unburnt 


Fig.  2.  Normalized  displacement  velocity  (red)  and  leading  point  location  (blue) 
time  history  profiles  at  100  Hz.  (For  interpretation  of  the  references  to  color  in  this 
figure  legend,  the  reader  is  referred  to  the  web  version  of  this  article). 


mixture  density  pu  were  obtained  from  laminar  flame  speed  cal¬ 
culations  using  the  FlameMaster  code  [29].  The  composition  and 
temperature  boundary  conditions  for  the  laminar  flame  speed  cal¬ 
culations  were  based  on  the  sampled  mixture  fraction  at  the  lead¬ 
ing  point  and  linearly  interpolated,  in  the  mixture  fraction  space, 
between  the  corresponding  inlet  values  of  the  fuel  and  coflow 
streams. 

Following  the  above  procedure,  displacement  velocities  were 
calculated  for  all  three  oscillation  frequency  cases,  with  20  points 
per  cycle,  to  demonstrate  their  evolution.  Furthermore,  as  shown  in 
Fig.  2,  SdISL  for  the  two  steady  cases  (2.4  and  8.0  m/s)  were  sim¬ 
ilarly  calculated  to  validate  this  definition  of  normalized  displace¬ 
ment  velocity  and  differentiate  between  tribrachial  flame  and  au¬ 
toignition.  For  clarity,  only  the  100  Hz  case  is  included  in  Fig.  2  to 
elucidate  its  evolution,  and  the  effect  of  the  oscillation  frequency 
will  be  discussed  in  Section  3.3. 

The  normalized  displacement  velocities  for  the  steady  autoigni¬ 
tion  front  and  tribrachial  flame  are  shown  in  Fig.  2  as  the  top  and 
bottom  horizontal  lines,  respectively.  The  Sd/SL  for  the  steady  tri¬ 
brachial  flame  is  around  unity,  while  this  value  is  around  eight 
for  the  autoignition  front.  These  values  are  similar  to  those  in 
HCCI  combustion  studies  [26]  and  therefore  can  be  used  to  bench¬ 
mark  the  unsteady  cases.  The  periodic  time  history  profile  of  Sd/SL 
is  bounded  by  but  does  not  fully  reach  the  two  steady  values, 
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Fig.  3.  Normalized  displacement  velocities  at  various  inlet  velocities  for  two  steady 
cases  and  the  100  Hz  oscillating  unsteady  cases. 


indicating  that,  while  the  chemical  structure  responds  to  the  flow 
dynamics,  such  response  is  not  fast  enough  to  reach  steady-state. 

When  SdISL  approaches  the  tribrachial  flame  limit,  its  value  is 
almost  constant,  while  the  change  near  the  autoignition  limit  is 
more  sinusoidal.  Moreover,  SdISL  changes  more  abruptly  when  the 
combustion  mode  switches  from  tribrachial  flame  to  autoignition. 
Compared  to  the  profile  of  the  normalized  leading  point  location 
(xd/D),  which  is  almost  sinusoidal,  the  profile  of  the  normalized 
displacement  velocity  is  asymmetric,  indicating  that  the  transition 
from  tribrachial  flame  to  autoignition  as  the  inlet  velocity  increases 
is  not  an  exact  reverse  process  of  the  transition  from  autoignition 
to  tribrachial  flame.  Indeed,  as  shown  in  Fig.  1,  although  the  in¬ 
let  velocities  at  0.25  and  0.75  cycle  are  the  same,  the  structures 
demonstrate  different  morphologies  during  the  cycle  of  decreasing- 
and  increasing-velocity:  there  is  hysteresis  during  the  transition. 

Such  hysteresis  is  demonstrated  more  clearly  in  Fig.  3,  where 
SdISL  is  plotted  against  the  inlet  velocity.  Given  the  same  inlet 
velocity,  the  reacting  fronts  have  different  displacement  velocities 
during  the  cycle  of  decreasing-  and  increasing-velocity.  Additional 
evidence  of  hysteresis,  shown  in  Fig.  1,  is  the  shift  in  the  location 
of  the  leading  point  in  the  mixture  fraction  space:  Z  =  0.14  when 
the  tribrachial  flame  dominates  and  Z  =  0.24  when  autoignition 
dominates.  The  shift  in  the  leading  point  mixture  fraction  as  well 
as  the  displacement  velocity  indicates  different  dominant  chemi¬ 
cal  reactions,  and  analysis  of  the  dominant  chemical  reactions  will 
reveal  the  mechanism  of  the  hysteresis. 

From  the  steady  case  analysis  [14],  the  dominant  chemical  path¬ 
ways  are  found  to  be  different  at  the  leading  point  of  the  tri¬ 
brachial  flame  and  autoignition  front.  Specifically,  the  hydrogen 
peroxide  branching  reaction  (H202  +  M  OH  +  OH  +  M)  is  the 
dominant  chain  branching  reaction  at  the  leading  point  of  the  au¬ 
toignition  front,  while  the  H  radical  branching  reaction  (H  +  02 
O  +  OH)  is  the  most  important  chain  branching  reaction  at  the 
tribrachial  flame  leading  point.  Due  to  the  longer  residence  time, 
hydrogen  peroxide  accumulation  is  much  higher  upstream  of  the 
autoignition  front  compared  to  the  tribrachial  flame  front. 

As  hydrogen  peroxide  plays  different  roles  in  the  tribrachial 
flame  and  autoignition  front,  its  spatial  profiles  along  the  Z  =  0.14 
and  Z  =  0.24  iso-contours  are  compared  in  Fig.  4,  with  the  left 
and  right  subfigures  corresponding  to  the  decreasing-velocity  and 
increasing-velocity  half  cycles,  respectively.  Qualitatively,  the  evo¬ 
lution  of  the  hydrogen  peroxide  profiles  shows  similar  trends  along 
both  mixture  fraction  iso-contours.  The  left  figures  show  that  hy¬ 


drogen  peroxide  accumulates  until  either  autoignition  occurs  or 
it  is  consumed  at  the  flame  front,  resulting  in  a  sharp  drop  in 
its  mass  fraction.  However,  depending  on  the  mixture  fraction, 
the  peak  value  of  the  hydrogen  peroxide  mass  fraction  differs  by 
three  to  five  times  between  a  steady  autoignition  front  (the  steady 
8.0  m/s  case)  and  a  tribrachial  flame  (the  steady  2.4  m/s  case), 
which  implies  its  different  significance  in  these  two  combustion 
modes  and  sets  the  benchmark  for  the  unsteady  evolution.  As  the 
inlet  velocity  decreases  from  8.0  m/s,  the  peak  Vh2o2  almost  re¬ 
mains  constant,  indicating  that  the  chemical  structure  is  very  close 
to  the  steady  autoignition  case.  As  a  consequence,  the  dominant 
chemical  pathway  remains  H202  +  M  <<=>>  OH  +  OH  +  M,  and  au¬ 
toignition  is  the  dominant  combustion  process,  resulting  in  larger 
SdISL.  However,  as  the  flow  velocity  decreases,  a  larger  gradient  is 
achieved,  resulting  in  steeper  profiles  and  smaller  Sd/SL,  according 
to  Eq.  (1).  The  large  gradient  results  in  enhanced  back  diffusion 
from  the  reacting  front  to  the  unburnt  upstream  mixture  and  ulti¬ 
mately  drives  the  transition  into  a  tribrachial  flame. 

Even  when  the  inlet  velocity  reaches  the  minimum  2.4  m/s, 
which  is  the  same  as  the  steady  case,  the  reacting  front  contin¬ 
ues  to  move  upstream.  Inlet  velocity  changes  slowest  around  the 
half  cycle,  allowing  the  chemical  structure  to  respond  to  the  hy¬ 
drodynamic  changes.  As  shown  on  the  right  of  Fig.  4,  the  peak 
Yh2o2  decreases  from  0.5  to  0.65  of  the  cycle.  At  this  stage,  au¬ 
toignition  is  not  fully  activated,  since  the  peak  Yh2o2  *s  lower  than 
the  steady  autoignition  case.  However,  the  peak  Vh2o2  is  still  larger 
than  the  steady  tribrachial  flame.  Therefore,  the  Sd/SL  of  the  tri¬ 
brachial  structure  is  close  to  but  slightly  larger  than  a  steady  flame, 
for  it  is  propagating  into  a  partially  reacted  mixture. 

As  the  inlet  velocity  further  increases  from  0.65  to  0.75  of  the 
cycle,  the  propagation  speed  of  the  tribrachial  flame  cannot  keep 
up  with  the  flow  incoming  velocity,  and  the  flame  structure  is 
therefore  convected  downstream.  During  this  flame  blow-off  pro¬ 
cess,  SdISL  remains  essentially  constant,  which  is  demonstrated  as 
the  flattened  bottom  in  Fig.  3.  Meanwhile,  the  unburnt  mixture 
upstream  of  the  flame  accumulates  radicals  and  heat  as  it  moves 
downstream  and  eventually  triggers  autoignition,  indicated  by  a 
sudden  jump  of  Sd/SL.  The  time  difference  between  the  half  cycle 
(where  the  inlet  velocity  is  2.4  m/s)  and  the  last  sample  point  be¬ 
fore  the  sudden  jump  of  SdISL  is  defined  as  the  induction  period, 
as  indicated  in  Fig.  3. 

3.3.  Effects  of  oscillation  frequency 

To  better  understand  the  coupling  between  hydrodynamics  and 
chemistry  and  the  hysteretic  behavior  in  Fig.  3,  the  effects  of  oscil¬ 
lation  frequency  are  analyzed.  Compared  to  the  100  Hz  case,  the 
other  two  cases  of  lower  oscillation  frequency  are  similar  qual¬ 
itatively  but  with  some  quantitative  differences.  Figure  5  shows 
the  evolution  of  the  normalized  leading  point  location  (xd/D)  for 
the  three  oscillating  cases,  compared  with  the  steady  benchmark 
cases.  The  peak-to-peak  variation  represents  the  oscillation  ampli¬ 
tude  during  a  cycle.  Two  distinct  trends  are  evident.  First,  as  the 
oscillation  frequency  increases,  the  oscillation  amplitude  decreases. 
Extrapolating,  if  the  oscillation  frequency  were  much  faster  than 
any  chemical  or  transport  times  scale,  the  oscillation  amplitude 
would  be  zero,  since  the  thermal  structure  could  not  respond  to 
the  velocity  changes. 

Second  and  more  interestingly,  for  all  three  frequencies,  the 
most  downstream  points  during  the  oscillation  are  almost  identi¬ 
cal  to  the  steady  case  at  8.0  m/s.  However,  as  the  frequency  in¬ 
creases,  the  most  upstream  point  during  the  oscillation  deviates 
more  from  the  steady  case  at  2.4  m/s.  The  different  effects  of  os¬ 
cillation  frequency  on  the  most  downstream  and  upstream  points 
are  attributed  to  the  difference  in  the  combustion  mode  at  these 
two  locations.  As  demonstrated  in  Fig.  3,  the  100  Hz  oscillating 
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Fig.  4.  Comparison  of  hydrogen  peroxide  mass  fraction  profiles  along  the  Z  =  0.14  and  Z  =  0.24  iso-contours  at  steady  state  and  at  100  Hz  during  the  decreasing-velocity 
cycle  (left)  and  increasing-velocity  cycle  (right). 


Cycle 

Fig.  5.  Normalized  leading  point  location  time  history  profiles  at  25,  50,  and 
100  Hz. 


case  has  already  established  quasi-steady  state  when  the  bound¬ 
ary  velocity  is  at  8.0  m/s  when  the  combustion  mode  is  kinetically 
controlled  by  autoignition.  In  a  Lagrangian  sense,  the  location  of 
an  autoignition  front  is  directly  related  to  the  ignition  delay  time, 


determined  by  chemical  kinetics.  Therefore,  the  location  of  the  au¬ 
toignition  front  at  this  quasi-steady  state  is  almost  the  same  as  the 
steady  case  with  the  same  boundary  velocity,  which  can  be  crudely 
estimated  as  the  product  of  the  ignition  delay  time  and  the  bound¬ 
ary  velocity.  For  even  lower  frequency  cases,  such  a  quasi-steady 
state  is  even  easier  to  achieve.  Consequently,  for  all  three  oscilla¬ 
tion  frequencies,  the  most  downstream  points  remain  very  close 
to  the  steady  autoignition  governed  case  with  the  boundary  ve¬ 
locity  of  8.0  m/s.  Conversely,  the  stabilization  mechanism  for  the 
steady  case  at  2.4  m/s  is  kinematically  controlled  by  the  balance 
between  the  tribrachial  flame  propagation  speed  and  the  incoming 
flow  speed.  For  the  100  Hz  case,  however,  such  kinematic  balance 
never  achieves  quasi-steady  state,  as  discussed  in  the  previous  sec¬ 
tion.  The  location  of  the  tribrachial  flame  is  then  determined  kine¬ 
matically  both  by  the  velocity  difference  between  the  unsteady  tri¬ 
brachial  flame  propagation  speed  and  the  flow  speed  and  the  time 
allowed  for  such  displacement  to  occur.  At  lower  oscillation  fre¬ 
quencies,  longer  time  is  allowed  for  the  tribrachial  flame  to  prop¬ 
agate  upstream  into  the  partially  reacted  mixture,  while  the  dis¬ 
placement  velocity  decreases  gradually  due  to  reduced  reactivity 
upstream.  The  most  upstream  point  location  will  then  asymptot¬ 
ically  approaches  the  steady  state  case  when  the  oscillation  fre¬ 
quency  is  sufficiently  low.  This  explains  why  the  minimum  xd/D 
approaches  the  steady  state  case  at  2.4  m/s  more  closely  at  lower 
oscillation  frequency  in  Fig.  5. 

Besides  the  oscillation  amplitude,  the  hysteretic  behavior  is 
also  affected  by  the  oscillation  frequency.  As  shown  in  Fig.  6,  the 
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Fig.  6.  Normalized  displacement  velocities  at  various  inlet  velocities  for  four  steady 
cases  and  three  oscillating  unsteady  cases  at  different  frequencies.  Induction  periods 
for  the  unsteady  cases  can  also  been  defined  similar  to  Fig.  3  but  are  not  shown 
here  for  clarity. 


Table  1 

Induction  time  at  different  oscillation  frequencies. 


Frequency  [Hz] 

25 

50 

100 

Induction  time  [ms] 

6 

4 

2.5 

hysteresis  of  decreasing  and  increasing  velocity  is  diminished  as 
the  oscillation  frequency  decreases,  denoted  by  the  shrinking  of 
the  enclosed  area.  Hysteresis  remains  at  slower  inlet  velocities 
since  the  decreasing-velocity  branch  is  still  autoignition  domi¬ 
nated,  but  it  takes  finite  induction  time  for  the  increasing-velocity 
unburnt  mixture  to  achieve  autoignition.  At  higher  inlet  veloci¬ 
ties,  both  branches  are  autoignition  dominant  and  therefore  col¬ 
lapse  to  a  single  path  at  sufficiently  low  frequency,  approaching  the 
quasi-steady  limit.  Ideally,  the  quasi-steady  limit  could  be  achieved 
by  investigating  even  lower  frequencies,  such  as  0.1  Hz;  however, 
such  calculation  would  be  extraordinarily  intensive,  requiring  more 
than  20  million  core-hours.  Instead,  three  steady  state  cases  re¬ 
ported  previously  [14]  and  a  new  steady  case  at  the  boundary 
velocity  of  6.4  m/s  are  also  included  in  Fig.  6  to  illustrate  the 
steady  state  limit.  The  3.2  m/s  case  is  kinematically  stabilized,  sim¬ 
ilar  to  the  2.4  m/s  case,  and  therefore  is  closer  to  the  increasing- 
velocity  branches  of  the  oscillating  cases  but  is  different  from  the 
decreasing-velocity  branches  for  the  same  reason  at  the  lower  ve¬ 
locity  condition,  as  discussed  in  the  previous  section.  Conversely, 
the  6.4  m/s  case  is  kinetically  stabilized  by  autoignition,  which 
can  be  approached  by  both  the  decreasing-velocity  and  increasing- 
velocity  branches  at  lower  oscillation  frequency. 

In  terms  of  the  relative  portion  of  the  oscillation  cycle, 
the  100  Hz  case  demonstrates  a  more  pronounced  hysteresis 
and  longer  induction  period.  However,  in  term  of  the  absolute 
time,  shown  in  Table  1,  the  25  Hz  case  takes  a  longer  time  to 
achieve  the  transition  from  small  Sd/SL  to  larger  values  during 
the  increasing-velocity  portion  of  the  cycle  ,  noting  that  the  time 
intervals  between  sampling  points  in  Fig.  6  are  not  the  same  for 
the  three  frequencies,  being  0.5,  1.0,  and  2.0  ms  for  100,  50,  and 
25  Hz,  respectively. 

This  counter-intuitive  finding  is  explained  with  Fig.  7.  Here,  the 
temperature,  hydrogen  peroxide,  and  methoxymethylperoxy  radical 
(CH30CH202)  profiles  are  compared  for  the  three  frequency  cases 
for  both  decreasing-  and  increasing-velocity  branches  at  the  same 
inlet  velocity  of  3.2  m/s  and  benchmarked  with  the  correspond¬ 


ing  steady  computation.  According  to  previous  studies  [13,14],  this 
steady  case  is  autoignition  dominated  and  stabilized  at  Z  =  0.24. 

According  to  the  temperature  profiles,  the  steady  and  unsteady 
cases  all  demonstrate  two  heat  release  stages,  with  the  first  and 
second  stage  correlated  well  with  the  depletion  of  the  CH30CH202 
radical  and  hydrogen  peroxide,  respectively.  The  CH30CH202  radi¬ 
cal  is  often  chosen  to  represent  the  NTC  chemistry  [12].  As  demon¬ 
strated  in  previous  studies  [13,14],  NTC  chemistry  is  important  in 
the  upstream  of  both  tribrachial  flame  and  autoignition  front,  and 
NTC  chemistry  is  still  important  for  the  unsteady  cases.  It  is  seen 
from  Fig.  7  that,  irrespective  of  the  oscillation  frequency  and  trav¬ 
eling  direction,  at  the  same  inlet  velocity,  Ych3och202  matches  with 
its  steady  state  profile,  indicating  that  NTC  chemistry  responds  to 
the  flow  oscillation  relatively  fast,  and  therefore  is  decoupled  from 
flow  dynamics. 

Conversely,  the  coupling  between  fluid  dynamics  and  second- 
stage  autoignition/flame  chemistry  is  important,  for  the  hydro¬ 
gen  peroxide  profiles  of  the  unsteady  cases  fail  to  collapse  onto 
that  of  the  steady  case.  Although  three  unsteady  profiles  on  the 
decreasing-velocity  branch  show  similar  peak  values,  the  lower  the 
frequency,  the  closer  the  profile  matches  the  steady  case,  which 
is  expected  from  Fig.  6.  However,  the  unsteady  profiles  on  the 
increasing-velocity  branch  show  lower  peak  values  than  the  steady 
structure,  with  the  100  Hz  case  profile  being  closer  to  the  steady 
counterpart.  Noting  from  Fig.  6  that  the  minimum  Sd/SL  in  all  three 
cases  are  larger  than  the  steady  2.4  m/s  case  and  at  lower  fre¬ 
quency  the  unburnt  mixture  has  longer  time  to  relax  to  steady 
state,  the  25  Hz  case  should  match  the  2.4  m/s  steady  case  closer, 
that  is,  a  tribrachial  flame  with  low  H202  accumulation.  Therefore, 
a  longer  time  is  required  to  accumulate  this  lost  H202  and  activate 
autoignition.  Conversely,  the  100  Hz  case  allows  less  time  to  relax 
to  a  steady  flame  structure  and  therefore  has  larger  Yh2o2  t0  start 
with  during  the  increasing-velocity  cycle,  resulting  in  a  shorter  in¬ 
duction  time  for  autoignition. 

The  different  effects  of  oscillation  frequency  on  the  first-stage 
NTC  chemistry  and  second-stage  autoignition/flame  chemistry 
can  be  understood  by  comparing  the  characteristic  time  scales  of 
these  processes.  The  ignition  delay  time  for  the  first-stage  ignition 
facilitated  by  NTC  chemistry  is  relatively  short  (~0.3  ms  [14]) 
compared  to  the  major  autoignition  process  induced  by  the 
hydrogen  peroxide  branching  reaction  (~1  ms  [14])  and  the  char¬ 
acteristic  hydrodynamic  oscillation  time  (10-40  ms  in  the  current 
work).  Therefore,  for  the  low  frequency  oscillations  investigated 
in  the  current  study,  the  fluid  dynamics  and  second-stage  au¬ 
toignition/flame  chemistry  coupling  is  mainly  responsible  for  the 
hysteretic  behavior  and  deviation  from  the  steady  state. 

At  even  higher  oscillation  frequencies,  the  NTC  chemistry  could 
interact  with  the  hydrodynamics.  Such  conditions  were  attempted 
computationally,  but  the  elevated  frequency  resulted  in  vortex 
shedding  and  local  flame  extinction,  which  add  further  complexity. 

4.  Conclusions 

Axisymmetric  laminar  nonpremixed  DME  coflow  flames  at  ele¬ 
vated  temperatures  and  pressures  with  sinusoidally  oscillating  in¬ 
let  velocities  were  computationally  investigated.  The  inlet  velocity 
oscillates  between  2.4  and  8.0  m/s  at  25,  50,  and  100  Hz.  Flame 
dynamics  in  such  oscillating  flows  and  frequency  effects  on  the 
hydrodynamics-chemistry  coupling  were  analyzed. 

The  heat  release  rate  profiles  were  examined  to  describe  the 
thermal  structure.  The  morphology  of  the  thermal  structure  tran¬ 
sitions  between  tribrachial  and  multibrachial.  The  multibrachial 
structure  is  favored  when  the  inlet  velocity  is  higher,  although 
there  is  hysteresis  during  the  transition.  Such  structures  agree  well 
with  the  steady  cases  in  Deng  et  al.  [14],  which  correspond  to 
different  combustion  modes:  tribrachial  flame  and  autoignition. 
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Fig.  7.  Temperature,  methyoxymethyiperoxy  radical,  and  hydrogen  peroxide  mass  fraction  profiles  along  the  Z  =  0.24  iso-contour  at  3.2  m/s.  DV:  decreasing-velocity,  and  IV: 
increasing-velocity. 


Normalized  displacement  velocity  was  defined  to  differentiate 
these  two  modes  in  the  current  study  and  compare  with  the  steady 
cases. 

According  to  the  steady  results,  the  normalized  displacement 
velocity  for  a  tribrachial  flame  is  around  unity  and  is  larger  for  au¬ 
toignition.  The  same  criterion  was  applied  to  the  unsteady  cases 
to  elucidate  the  evolution  of  combustion  mode.  As  the  inlet  veloc¬ 


ity  decreases,  autoignition  is  the  dominant  combustion  process  un¬ 
til  flame  chemistry  takes  over  around  the  most  upstream  location 
and  slowest  inlet  velocity.  The  tribrachial  flame  is  convected  down¬ 
stream  as  the  flow  velocity  increases.  The  radical  and  heat  accumu¬ 
lation  upstream  of  the  tribrachial  flame  finally  results  in  autoigni¬ 
tion,  showing  a  sudden  increase  in  the  normalized  displacement 
velocity. 
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Oscillation  frequency  effects  on  the  hydrodynamics-chemistry 
coupling  were  analyzed  by  examining  the  profiles  of  temperature, 
methyoxymethylperoxy  radical,  and  hydrogen  peroxide  during  the 
oscillation  process.  It  is  found  that,  at  the  three  frequencies  in¬ 
vestigated,  the  tribrachial  structure  does  not  have  sufficient  time 
to  reach  steady  state,  and  the  transition  from  tribrachial  flame  to 
a  multibrachial  autoignition  front  occurs  over  a  finite  induction 
time  as  velocity  increases.  Consequently,  the  decreasing-velocity 
and  increasing-velocity  cycles  have  different  normalized  displace¬ 
ment  velocities  and  hence  demonstrate  hysteresis.  At  lower  fre¬ 
quencies,  such  hysteresis  is  less  pronounced,  for  longer  relaxation 
time  is  allowed  to  approach  the  quasi-steady  state  condition. 

NTC  chemistry  represented  by  methyoxymethylperoxy  radical 
accumulation  and  depletion  has  shorter  time  scales,  and  there¬ 
fore  is  able  to  respond  to  the  hydrodynamic  changes.  However,  au¬ 
toignition  and  flame  establishment  have  comparable  time  scales  to 
the  oscillation  period  and  are  therefore  coupled  with  the  flow  dy¬ 
namics.  At  lower  oscillation  frequencies,  chemical  kinetics  is  closer 
to  reaching  the  two  steady-state  conditions  with  the  maximum 
and  minimum  boundary  velocities. 
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Combustion  instability  of  a  lifted  triple  flame  downstream  of  the  injector  nozzle  in  an  acoustic  field  is 
studied,  with  implication  on  the  combustion  instability  in  liquid  rocket  engines.  The  analysis  is  focused 
on  the  premixed  flame  segment,  and  is  based  on  the  linear  equation  for  a  perturbed  flamefront.  Assuming 
that  the  radius  of  the  flamefront  curvature  is  much  larger  than  that  of  the  internal  flame  structure,  the 
premixed  front  is  approximated  as  being  quasi-planar.  Both  acoustical  stabilization  of  the  hydrodynamic 
flame  instability  and  excitation  of  the  parametric  instability  by  sound  waves  are  considered,  and  the 
limitation  for  stable/unstable  regimes  as  a  function  of  flame  parameters  and  acoustic  frequency/intensity 
is  determined.  The  velocity-wise  nature  of  the  sound  is  scrutinized,  namely,  the  acoustics  are  imitated 
by  a  pocket  of  harmonic  flow  oscillations.  First,  the  Landau  limit  of  infinitely  thin  flames  is  employed. 
Then  the  formulation  is  extended  to  account  for  the  finite  flame  thickness,  and  it  is  demonstrated  how 
unstable  regimes  are  modified  by  the  flame  thickness. 
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1.  Introduction 

Combustion  instabilities  in  liquid  rocket  engines  (LRE)  are 
driven  by  various  fundamental  processes  including  mixing,  flame 
dynamics,  chemical  kinetics  and  spray  formation  [1].  Besides,  the 
dynamics  of  nozzle  injection  as  well  as  the  subsequent  fuel- 
oxidizer  mixing  and  combustion  can  also  contribute  to  the  insta¬ 
bility  onset.  In  particular,  a  resonance  can  develop  through  acous¬ 
tic  coupling  between  the  chamber  dynamics  and  the  injection  and 
combustion  processes  [2].  It  is  generally  believed  that  combustion 
within  LREs  likely  proceeds  in  the  diffusion-flame  mode  and,  due 
to  the  high  density  of  the  medium,  the  associated  reactions  occur 
infinitely  fast  as  compared  to  the  processes  of  heat  and  mass  trans¬ 
port.  Consequently,  the  majority  of  studies  on  the  LRE  combustion 
instabilities  are  focused  on  diffusion  flame  processes  in  the  limit  of 
flame-sheet  combustion. 

Generally  speaking,  four  possible  flame  stabilization  mecha¬ 
nisms  have  been  proposed  for  the  bulk  diffusion  flame,  result¬ 
ing  from  the  separately  injected  fuel  and  oxidizer,  namely,  (i)  that 
through  the  presence  and  extinction  of  diffusion  flamelets  gen¬ 
erated  by  mixing  [3];  (n)  stabilization  by  the  recirculation  zone 
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around  the  injection  nozzle  [4,5];  (in')  that  through  auto-ignition 
in  a  turbulent  flow  field,  occurring  when  one  of  the  reactants  is 
introduced  at  a  high  temperature  [6];  and,  eventually,  (iv)  stabi¬ 
lization  of  the  bulk  flame  through  a  triple  flame  structure  at  the 
base  of  the  flame  [7],  which  constitutes  the  focus  of  the  present 
paper. 

To  be  more  specific  with  the  latter,  4th  stabilization  mecha¬ 
nism,  of  Chung  [7],  we  consider  a  flamefront  propagating  through 
a  fuel/air  mixing  layer,  where  the  fuel  concentration  is  stratified 
from  lean  to  rich.  The  resulting  flame  exhibits  a  tribrachial  (triple) 
structure,  consisting  of  a  lean  and  a  rich  premixed  flame  wings 
together  with  a  trailing  diffusion  flame,  all  extending  from  a  sin¬ 
gle  point  [8].  The  propagation  speed  and  dynamics  of  these  flame 
segments  are  dominantly  influenced  by  the  fuel  concentration  gra¬ 
dient,  or  the  mixture  fraction  gradient  upstream  of  the  intersec¬ 
tion  point,  since  the  local  laminar  burning  velocities  along  the  pre¬ 
mixed  flame  wings  depend  sensitively  on  the  fuel  concentration, 
and  the  concentration  gradient  determines  the  effective  thickness 
of  a  flammable  region  and,  thereby,  the  curvature  of  the  premixed 
wings  [7].  When  the  flow  rate  from  the  nozzle  is  small,  the  flame 
is  attached  to  the  rim  of  the  nozzle.  In  contrast,  when  the  flow  rate 
is  sufficiently  large,  the  flame  can  be  lifted  off  as  shown  schemati¬ 
cally  in  Fig.  1. 

Stabilization  and  the  structure  of  triple  flames  were  addressed 
in  Refs.  [9]  and  [10],  which  were  followed  by  the  classical  pa¬ 
pers  [8,11],  see  also  Refs.  [1]  and  [12]  for  details.  In  particular,  the 
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stabilization  criterion  have  been  determined  in  Ref.  [11]  recogniz¬ 
ing  that  the  point  of  stabilization  for  a  triple  flame  is  a  locus 
where  the  lean-premixed,  rich-premixed  and  diffusive  flame  seg¬ 
ments  meet,  and  where  the  fuel/oxidizer  equivalence  ratio  is  close 
to  stoichiometric.  Based  on  the  fact  that  the  local  flame  speed 
should  balance  the  flow  velocity,  Ref.  [11]  calculated  the  stoichio¬ 
metric  laminar  flame  speed  SL  as  a  function  of  the  flame  parame¬ 
ters  such  as  the  Schmidt  number  Sc  and  the  fuel  mass  fraction  Y, 
as  well  as  the  flow  parameters,  such  as  the  nozzle  diameter  d,  the 
inlet  flow  mean  velocity  u0t  and  the  flow  velocity  profile. 

The  formulation  [11]  has  been  extended  in  a  series  of  studies.  In 
particular,  the  substantial  influence  of  the  thermal  expansion  and 
mixture  fraction  on  the  triple  flame  speed  has  been  found  [13],  and 
the  stretch  effect  [14]  as  well  as  the  heat  release  [15]  were  incor¬ 
porated  into  the  investigation.  However,  an  analysis  of  the  stabil¬ 
ity  of  the  triple  flame  through  that  of  the  leading,  premixed  flame 
segment,  when  subjected  to  imposed  oscillations  like  external  or 
flame-generated  acoustics,  has  not  been  performed.  In  this  regard, 
we  note  that  analyses  of  the  hydrodynamic  instability  of  premixed 
flames  were  performed  in  Refs.  [16-19]  and  subsequently  extended 
to  flame-sound  interactions  in  Refs.  [20-23].  It  was  shown  that 
high-frequency  acoustics  may  either  stabilize  or  amplify  the  hy¬ 
drodynamic  instability,  depending  on  various  flame  and  sound  pa¬ 
rameters.  Furthermore,  a  flame  accelerated  normal  to  itself,  such  as 
due  to  flow  oscillations  normal  to  itself,  may  be  subjected  to  the 
parametric  instability.  In  addition,  it  is  noted  that  a  flame-acoustic 
resonance  may  lead  to  violent  folding  of  a  premixed  flamefront 
[23,24]  and,  presumably,  even  flame  turbulization  [25,26]. 

In  this  regard,  in  the  present  work  we  investigate  the  effect  of 
harmonic  flow  oscillations,  such  as  might  be  induced  by  acoustic 
waves,  on  the  stability  of  tribrachial  lifted  flames.  For  this  purpose, 
the  structure  of  a  tribrachial  flame  [8,11],  the  dispersion  relation 
for  the  intrinsic  flame  instability  [17-19],  and  the  flame-acoustic 
interactions  [2,20-23]  are  combined.  It  is  noted,  in  this  respect, 
that  complementary  to  a  tangential  component  to  the  flamefront, 


such  as  it  would  in  the  premixed  flame  segments,  the  flow  may 
also  involve  wrinkles  propagating  along  the  flame  surface. 

The  analysis  is  presented  in  the  following  four  sections.  The 
generalized  formulation  on  the  coupling  of  an  imposed  har¬ 
monic  flow  oscillation  to  an  intrinsic  instability  mode  is  presented 
in  Section  2.  The  local  flame  speed  is  described  in  Section  3. 
The  development  of  the  hydrodynamic  instability  in  the  pres¬ 
ence  of  sound  waves,  either  external  or  flame-generated,  as  well 
the  acoustically-excited  parametric  instability,  are  analyzed  in 
Section  4  within  the  Landau  limit  of  zero  flame  thickness.  The  ef¬ 
fect  of  finite  flame  thickness  is  considered  in  Section  5.  The  paper 
ends  with  concluding  remarks. 

2.  Generalized  formulation 

We  consider  a  discontinuity  surface  of  zero  or  finite  thickness 
separating  two  fluids  of  different  properties.  The  surface  is  corru¬ 
gated  by  weak,  linear,  two-dimensional  perturbations  described  by 
the  function  z  =  F(x,  t )  =  J2fi(t)  cos (/qx).  Hydrodynamic  equations 
ahead  and  behind  the  front,  the  matching  conditions  across  the 
front,  and  the  reaction  evolution  equation  can  be  incorporated  into 
a  linear  dispersion  relation  for  the  perturbation  amplitude  /  in  a 
general  form  If  =  0,  where  the  instability  operator  /  describes  the 
surface  response  to  the  perturbation,  and  it  couples  the  instabil¬ 
ity  growth  rate  o  to  the  perturbation  wavenumber  k,  I  =  I  (a,  /<), 
which  can  also  include  a  multitude  of  other  hydrodynamic  and 
thermal-chemical  parameters  inherent  to  the  associated  instabil¬ 
ity  mode.  If  the  system  is  subjected  to  an  external  oscillating  field, 
such  as  that  of  acoustics,  the  dispersion  relation  reads  If  =  i>a/. 
It  is  assumed  that  the  forcing  operator  i>a  is  characterized  mainly 
by  a  dominant  acoustic  mode  of  frequency  Q,  and  amplitude  Ua, 
i>a  =  (£2,  Ua),  which  are  determined  by  the  combustion  chamber 
dynamics  for  implications  to  the  amplitude  and  frequency  of  the 
LRE  combustion  instability. 

Within  linear  approximation,  there  are  two  interacting  fre¬ 
quency  modes  between  the  intrinsic  instability  and  the  acoustic 
oscillations,  namely,  (i)  a  resonant  mode,  Q  =  co,  which  modifies 
and  weakens  the  intrinsic  instability,  and  (n)  a  doubling  mode, 
£2  =  2 co,  which  can  excite  the  parametric  instability  through  the 
intrinsic-acoustic  coupling  if  the  sound  is  strong  enough.  For  the 
resonant  mode,  we  seek  the  solution  to  the  dispersion  relation 
If  =  i>a/  in  the  form  /  =  /i  (t)  +  f2(t)  cos  (£2t),  with  fx  and  f2  be¬ 
ing  “slow  functios”  that  vary  considerably  during  a  time  interval 
much  larger  than  the  oscillation  period,  2ttI£2.  For  the  paramet¬ 
ric  mode,  a  partial  solution  to  the  dispersion  relation  acquires  the 
form  /  =  {/3  cos(£2t/2)  +  /4  sin(£2t/2)}  exp(crt).  In  both  branches 
we  then  arrive  at  the  stability  criterion  given  by  Re{a}  <  0. 

This  formulation  will  be  systematically  employed  to  quantify 
the  coupling  between  various  intrinsic  flamefront  instability  modes 
and  the  chamber  acoustics.  The  instability  modes  of  interest  in¬ 
clude,  though  not  limited  to:  (i)  the  Rayleigh-Taylor  (RT;  body- 
force)  instability  for  both  premixed  and  diffusion  flame  segments; 
(ii)  the  Kelvin-Helmholtz  (KH;  shear-layer)  instability  which  can 
be  mostly  important  for  diffusion  flames  because  of  the  inher¬ 
ently  different  properties  of  the  fuel  and  oxidizer  jets;  (Hi)  the 
diffusional-thermal  (DT;  nonequidiffusive)  instability,  which  not 
only  is  pervasive  for  premixed  flames,  with  its  intensity  diminish¬ 
ing  with  increasing  pressure  as  the  flame  thickness  reduces,  but  is 
also  relevant  for  near-limit  diffusion  flames;  and  (zv)  the  Darrieus- 
Landau  (DL;  hydrodynamic)  instability,  inherent  to  all  practical 
premixed  flames.  Furthermore,  unlike  the  DT  mode,  the  DL  in¬ 
stability  intensifies  with  increasing  pressure  and  thereby  is  espe¬ 
cially  relevant  to  LRE  premixed  flame  segments.  All  these  insta¬ 
bility  modes  and  their  relevance  to  the  diffusive  and  premixed 
flame  segments  are  summarized  in  Table  1.  This  approach  has 
recently  been  employed  to  the  KH  instability  in  a  non-reacting 
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Table  1 

Intrinsic  and  acoustic  flamefront  instability  modes. 


Instability  mode 

Premixed  flame 

Diffusion  flame 

Rayleigh-Taylor  (RT) 

Yes 

Yes 

Kelvin-Helmholtz  (KH) 

? 

Yes 

Diffusional-thermal  (DT) 

Pervasive,  |  as  p  f 

Near-limit  situations 

Darrieus-Landau  (DL) 

t  as  p  t  and 

No 

//////////////////////// 

//////////////////////// 

//////////////////////// 

Acoustical/parametric 

Yes 

Yes 

share-flow  [27].  In  the  rest  of  this  paper,  we  shall  focus  on  the 
4th,  DL  instability  and  its  coupling  to  a  resonant  harmonic  flow 
oscillation  imitating  an  imposed  or  flame-generated  sound  wave. 

3.  Local  flame  speed  at  premixed  flame  segments 


In  this  section,  we  brifely  recall  Chung’s  theory  of  bulk  flame 
stabilization  through  a  triple  flame  structure  at  the  base  of  the 
flame  ensemble  [7].  Upon  mixing  between  the  fuel  and  oxidizer 
downstream  of  a  nozzle,  a  concentration  gradient  is  established 
from  the  fuel  side  to  the  oxidizer  side,  leading  to  the  formation  of 
a  premixed  flame  segment  that  varies  from  fuel-lean  to  fuel-rich 
over  its  surface.  Downstream  of  this  premixed  flame  segment,  the 
unburned  oxidizer  from  the  fuel-lean  segment  further  reacts  with 
the  unburned  fuel  from  the  fuel-rich  segment,  forming  a  diffusion 
flame,  which  eventually  merges  with  the  bulk  diffusion  flame.  The 
flame  structure  therefore  consists  of  a  lean-premixed  flame  (LPF) 
segment,  a  rich-premixed  flame  (RPF)  segment,  and  a  downstream 
diffusion  flame  (DF)  segment,  as  shown  in  Fig.  1.  It  is  believed 
that  stabilization  of  this  structure  is  achieved  through  that  of  the 
concentration-stratified  premixed  flame  base,  at  the  stoichiometric 
location  on  its  surface.  It  is  emphasized  that  the  flame  dynamics 
near  the  attachment  point  has  an  essential  influence  on  the  over¬ 
all  flame  dynamics  for  both  premixed  and  diffusion  flames.  Specif¬ 
ically,  these  dynamics  essentially  appear  as  “boundary  conditions” 
for  the  rest  of  the  premixed  or  diffusion  flamelet  dynamics;  see, 
for  instance  Refs.  [28]  and  [29]  as  well  as  references  therein. 

According  to  Ref.  [8],  the  triple  flame  is  described  by  the  simi¬ 
larity  variable 


1 1  = 


4 Rer- , 
32  z 


(1) 


where  the  Reynolds  number  Re  =  u0d/v  is  related  to  the  inflow, 
and  the  coordinate  i)  =  i)S  along  the  stoichiometric  contour  Y  =  Ys 
may  be  expressed  as 


r  2Sc  + 1  I  ' 

,1/2 Sc 

"2Sc+l  Red' 

L  8  JZVZs  pYs. 

1  = 

_  32  Ys  zs  _ 

(2) 


with  the  mass  flux  of  the  fuel  being  I  =  7tpu0d2/4.  As  a  result,  we 
find 


_  16  zs  /  /  2Sc  +  1  Red  \  '/2Sc 

rs~VfRe)l  ) 


(3) 


If  the  first  term  under  the  square  root  in  Eq.  (3)  is  much  greater 
than  unity,  then  the  shape  of  the  stoichiometric  contour  Y  =  Ys  is 
described  by  the  explicit  self-similar  formula 


Zs  _ 

/  32  Ys  \<*-> 

/V3rs\ 

dRe 

V  2Sc  + 1  / 

l  16  d  j 

The  z-velocity  component  of  the  flow  and  the  fuel  mass  fraction 
along  the  similarity  variable  are  then  given  by 


3  J  1  3  Red  1 

8^s^(1+??2/8)2  “  32U°~zF  (i+??2/8)2’ 


2Sc + 1  J  1  2Sc + 1  Red  1 

s  “  8jFvzs  ~p  ^  +  ^2/8)2Sc  "  32  IT  (t  +  ’ 

with  the  momentum  flux  J  =  u0I  =  tt  pu^d2 /4,  and  the  lift-off 
height  of  the  flame  given  by  [8] 


Consequently,  we  have  HL  oc  d2i^2Sc-1^/(Sc-1),  where  the  propor¬ 
tionality  factor  is  determined  from  the  fuel  characteristics.  The 
lift-off  height  is  proportional  to  the  nozzle  cross  section,  and  the 
Schmidt  number,  Sc,  is  the  key  parameter  for  the  stability  of 
triple/lifted  flames,  in  that  a  stabilized  lifted  flame  can  exist  only 
for  Sc  <  0.5  or  Sc  >  1.  Eliminating  r]s  by  combining  Eqs.  (5)  and 
(6),  and  recognizing  that  the  stoichiometric  local  flame  speed  is 
balanced  by  the  axial  flow  velocity  at  the  stabilization  point,  SL  = 
u(Ys ),  we  obtain 


Together  with  Eq.  (4),  Eq.  (8)  yields 

Eq.  (9)  shows  SL  as  a  function  of  the  radial  coordinate  rs.  This 
is  a  given  mixture  parameter,  which  determines  the  flame  stabi¬ 
lization  condition.  It  is  emphasized  that  unlike  the  standard  one¬ 
dimensional  unstretched  laminar  flame  speed,  which  is  a  thermal- 
chemical  parameter,  the  local  quantity  SL  is  a  result  of  both  the 
flame  and  the  flow.  It  depends  on  the  fuel  properties  by  means  of 
Sc  and  Ys,  and  on  the  flow  parameters  via  d  and  u0.  Another  pa¬ 
rameter  characterizing  the  fuel  type  is  the  density  expansion  factor 
©  =  pu/ph.  The  latter  plays  the  key  role  in  the  present  work,  as 
will  be  shown  subsequently.  Obviously,  varying  the  fuel  type  we 
would  change  SL  and  ©,  and  thereby  modify  the  results. 

It  is  also  recognized  that  the  analysis  above,  based  on  the 
Chung  theory  [7,8]  does  not  consider  transition  to  turbulence, 
thereby  providing  the  limitation  on  the  Reynolds  number  that  can 
be  employed.  It  is  noted,  in  this  respect,  that  according  to  the  clas¬ 
sical  Hottel-Hawthorne  work  [30],  when  the  jet  velocity  is  low  and 
the  flow  is  laminar,  the  flame  height  increases  with  the  jet  velocity 
but  only  until  a  certain  maximum  at  which  the  instability  is  ob¬ 
served  at  the  flame  tip.  In  contrast,  with  a  further  increase  in  the 
jet  velocity,  the  flame  becomes  turbulent,  and  its  height  sharply 
decreases  and,  when  scaled  by  the  nozzle  diameter,  remains  al¬ 
most  constant  thereafter  [1,30].  Such  a  turbulent  flame  is  beyond 
the  present  work  and  will  be  considered  elsewhere. 


4.  Coupling  of  harmonic  flow  oscillations  to  premixed  flames 


Eq.  (3)  describes  a  “quasi-steady”  flame  shape  in  the  form  z5  = 
F0(rs).  Flowever,  the  steady  formulation  of  Section  3  does  not  ac¬ 
count  for  hydrodynamic  flame  instability  as  well  as  external  and/or 
flame-generated  acoustics  and  vortices,  which,  in  turn,  can  lead  to 
additional  destabilization  of  the  flame.  Accounting  for  the  fact  that 
the  radius  of  the  curvature  of  the  LPF/RPF  sheet  is  much  larger 
than  the  internal  flame  structure,  we  consider  the  entire  premixed 
front  as  quasi-planar;  and  the  local  mass  fraction  of  the  fuel  mix¬ 
ture  Y,  the  local  density  expansion  factor  0  =  pu/pb,  and  the  local 
planar  flame  speed  SL  are  assumed  to  be  constant.  For  simplicity, 
we  first  illustrate  our  formulation  in  the  Landau  limit  of  zero  flame 
thickness,  and  then  extend  it  to  that  of  finite  flame  thickness. 

We  consider  weak  perturbations  of  any  planar  segment  in  the 
form  z  =  F0(r)  +  F(r,  t),  with  F0(r)  given  by  Eq.  (4),  and  \F\  <<c  F0.  It 
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is  nevertheless  recognized  that  neglecting  the  radius  of  the  curva¬ 
ture  of  the  leading  edge  of  the  flamefront  could  potentially  miss 
some  physical  effects  such  as  the  triple  flame  may  actually  ap¬ 
pear  at  a  point  in  the  flow,  where  the  velocity  is  noticeably  larger 
than  that  evaluated  as  the  laminar  flame  speed  [29].  Furthermore, 
the  oscillatory  flow  disturbances  may  cause  an  oscillation  in  cur¬ 
vature  of  the  leading  edge  of  the  flame  and,  thereby,  an  oscilla¬ 
tory  modification  of  the  approach  flow  velocity.  However,  a  de¬ 
tailed  quantification  of  these  effects  requires  a  separate  analysis 
and  will  be  presented  elsewhere.  For  simplicity,  hereafter  we  deal 
with  a  two-dimensional  configuration,  though  the  formulation  can 
be  extended  to  a  three-dimensional  flow.  In  the  Landau  limit  of 
zero  flame  thickness,  the  basic  equations  in  a  two-dimensional  in¬ 
compressible  flow  (u,  w)  are  the  continuity  equation 


du  d  w 

15 - 1 — a —  = 

dz  dr 

the  Navier- Stokes  equations 


(10) 


du  du  du 

-FT-  +  U  -77 - h  W  ~p— 

dt  dz  dr 


1  dP 
p  dz 


+  vV2u  +gz, 


(11) 


9w  9w  9w  1  dP 
-or-  +  u^~  +w^-  =  —  q-  +  vV2w  +  gr, 
dt  dz  dr  p  dr 

the  matching  conditions 


(12) 


u+-U--(w+-w_)|^  =SL(0-l)yi  +  ( 

(  dF\2 
ydr)  ’ 

(13) 

dF 

w+-w_  +  (u+-u_)g^  =0, 

(14) 

P+-P-  =  -(@-1)Pu5l, 

and  the  flame  evolution  equation 

(15) 

9  F  dF  ~  1  ( 9F  + 

_+w___u_+Siyi  +  ^— j  . 

(16) 

Eqs.  (10)-(16)  can  be  incorporated  into  the  linear  dispersion  re- 

lation  for  the  perturbation  amplitude  F  =  f(t)  exp (ikr)  in  the  form 
[17-20] 


W  +  ~A0^k2f  =  Akgfcosa, 


(17) 


where  A  =  (0  -  l)/(©  +  1)  is  the  local  Atwood  number,  g  the  gen¬ 
eral  acceleration  field,  and  a  the  angle  between  the  acceleration 
vector  and  the  normal  vector  to  the  flame  sheet.  Eq.  (17)  describes 
the  linear  stage  of  the  Darrieus-Landau  and  Rayleigh-Taylor  insta¬ 
bility  modes  in  the  acceleration  field.  Without  external  accelera¬ 
tion,  g=  0,  Eq.  (17)  reduces  to  the  Darrieus-Landau  dispersion  re¬ 
lation  [16].  In  the  present  study,  the  acceleration  field  is  that  of 
acoustics.  The  sound  wave  is  characterized  by  a  dominant  mode 
uz  =  U  sin(<wt),  with  the  amplitude  U  and  frequency  co.  Then  Eq. 
(17)  becomes 


d2f 

dt2 


20 

0TT 


-  AQSfk2/  =  AkUco  cos  (cot)  f  cos  a. 


(18) 


We  consider  two  situations,  namely:  (a)  that  of  external  acous¬ 
tics,  when  U  and  co  are  free  parameters  of  the  problem,  and 
(j b)  of  flame-generated  acoustics,  when  U  and  co  depend  on  the 
flame/flow  parameters. 

It  is  noted  that  there  are  two  modes  of  the  flame-acoustic  in¬ 
teraction:  (i)  the  resonant  mode,  Q  =  co,  describing  the  develop¬ 
ment  of  the  Darrieus-Landau  (DL)  instability  in  the  presence  of 
acoustics,  and  thereby  possible  flame  stabilization  by  sound;  and 
(n)  the  mode  which  frequency  is  twice  that  of  the  oscillating  field, 
£2  =  2 co,  that  can  excite  the  parametric  instability  of  a  corrugated 
flame  segment. 


4.1.  Stabilization  of  flame  instability  by  a  resonant  harmonic 
oscillation 


Recognizing  that  the  stabilization  of  the  intrinsic  flamefront  in¬ 
stabilities  by  sound  waves  is  analogous  to  the  stabilization  of  the 
Kapitsa  pendulum  [16],  following  the  method  of  Bychkov  [21],  we 
seek  the  solution  to  Eq.  (18)  in  the  form 

/  =  /i(0  +  f2(t) cos  (cot),  (19) 


where  fa  and  f2  are  “slow  functions”  which  do  not  change  con¬ 
siderably  during  the  acoustic  time  27 r/co.  Then  d2f2/dt2  oc  co2f2  » 
SLk(df2/dt)  »  Sfk2f2,  and  substituting  Eq.  (19)  into  Eq.  (18),  we 
find 


d2h 

dt2 


co2  f2  cos  (cot) 


20  Cfcd/l 
0-lSi  dt 


A  0  Sfk2f i  -AcoUk  /i  cos  a  cos  (cot) 


Aco  U  k  f2  cos  a 


1  +  cos  (2 cot) 


=  0. 


(20) 


Separating  the  terms  of  Eq.  (20)  oscillating  with  frequency  co, 
we  obtain 


f2  =  -AUkfi  cos  a /co, 


(21) 


+  (lA2U2cos2a-A®Sl)k2h  =  0.  (22) 

Seeking  the  solution  to  Eq.  (22)  in  the  form  /1ocexp(<j1t),  we 
obtain  the  dispersion  relation 


<*\  = 


0 


0  +  1 


-1 


I  ©2  +  0  —  1  \  ( (0  -  1)  D  cos  a 

©  2 1  0  §r~ 


SLk. 


(23) 


Eq.  (23)  formally  does  not  include  the  acoustic  frequency  co  be¬ 
cause  it  was  derived  with  the  assumption  of  high-frequency  oscil¬ 
lations  as  compared  to  the  characteristic  DL  growth  rate.  It  is  seen 
from  Eq.  (23)  that  the  increase  in  the  acoustic  amplitude  reduces 
the  instability  growth  rate  0\.  The  instability  is  completely  sup¬ 
pressed  by  the  sound  waves  if  Re{cr <  0,  i.e.  when 


fUc  oscA2  ,^©  +  1  2© 

>29r'I' 


U  cos  a  >  3 


(-§^),,i(s+)*(^)i 


U  0, 


(24) 


(25) 


accounting  for  the  relation  (9).  Here  we  have  considered  the  sound 
waves  as  an  external  force  of  the  problem.  In  the  event  of  flame¬ 
generated  sound,  we  take  a  ~  0,  and  the  local  acoustic  amplitude 
can  be  estimated  by  means  of  the  local  corrugated  flame  speed  as 


u  ~  (0  -  1  )SLf  1  +  (dFo/dr)2.  (26) 

With  Eqs.  (4)  and  (26),  the  inequality  (24)  becomes 


(  4Sc 

32  Ys  | 

CV3\\~ 

2 

4Sc-l 

20(0  +  1) 

\4Sc  -  1 K ) 

2Sc  +  1  * 

V16/ 

(©-i)3 

The  result  (27)  shows  the  criterion  governing  the  suppression  of 
the  onset  of  the  hydrodynamic  instability  at  the  premixed  segment 
of  a  lifted  triple  flame  by  flame-generated  acoustics.  According  to 
Eq.  (27),  the  suppression  occurs  for  arbitrary  Re,  Sc  and  Ys  if  ©  > 
©c  «  4.84,  as  the  RHS  of  Eq.  (27)  is  negative  in  this  case. 
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4.2.  Parametric  instability 


Another  effect  of  flame-acoustic  interaction  is  the  development 
of  the  parametric  instability  of  a  premixed  flame  segment,  which 
is  excited  by  the  mode  of  twice  the  smaller  frequency  than  that 
of  sound.  The  parametric  instability  of  a  premixed  flamefront  in 
an  oscillating  field  resembles  the  classical  parametric  instability  of 
an  oscillator  [16,20].  Following  Bychkov  [21],  we  look  for  a  partial 
solution  to  Eq.  (18)  in  the  form 


/  =  [/3  cos  (^)  +  f4  sin  (^)J  exp  ( cr2t ), 


(28) 


where  the  coefficients  f3  and  f4  are  the  free  parameters,  and  o2  the 
growth  rate  of  the  parametric  instability.  Substituting  Eq.  (28)  into 
Eq.  (18),  collecting  the  terms  containing  cos(cut/2)  and  sin(&>t/2), 
and  omitting  terms  of  higher  frequencies,  we  find 
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(30) 


and  subsequently  arrive  at  the  dispersion  relation  for  the  paramet¬ 
ric  instability 
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©TT 
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The  parametric  instability  is  excited  if  Re{cr }  >  0,  i.e.  when 

2  /  _  _  \  2  /  o  . 


(31) 
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2© 
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In  the  event  of  flame-generated  sound,  the  acoustic  amplitude 
can  be  estimated  by  Eqs.  (4)  and  (26),  while  its  frequency  can  be 
estimated  by  means  of  the  speed  of  sound  cs  and  the  combustor 
width  Lcomb  or  the  lift-off  height  HL,  Eq.  (7), 


co  ~  27TCs/Lcomb ,  or  co  ~  2ttcs/Hl, 


(33) 


0  12  3  4 

0) /2(A0)l/2SLk 

Fig.  2.  Stability  diagrams  for  flames  with  small,  0  =  1.2  (a),  and  large,  0  =  8  (b), 
thermal  expansion  in  the  Landau  limit  of  zero  flame  thickness,  L  =  0. 


numbers.  According  to  Eqs.  (24)  and  (32),  there  is  no  “absolute  sta¬ 
bility  band”  if  0  >  \  (\  +  V5)  «  1.6,  as  demonstrated  in  Fig.  2b.  It 
is  noted,  in  this  respect,  that  this  threshold  thermal  expansion  as¬ 
sociated  with  the  absolute  stability  is  much  smaller  than  the  typi¬ 
cal  values  in  most  flames.  Therefore  the  absolute  stability  of  Fig.  2 
is  rarely  expected  to  occur  in  practical  combustors.  However,  the 
“conditional  stability”  with  respect  to  a  limited  range  of  pertur¬ 
bations  (limited  wavelengths/wavenumbers)  is  devoted  to  higher 
thermal  expansions,  including  those  occurring  in  the  practical  re¬ 
ality,  say  0  =  5-  8. 


5.  Effect  of  finite  flame  thickness 


respectively,  depending  on  whether  radial  or  axial  sound  waves  are 
considered.  Consequently,  Eqs.  (4),  (7),  (26),  (32)  and  (33)  consti¬ 
tute  the  criterion  for  the  parametric  instability  excitation  by  flame- 
induced  acoustics. 

Obviously,  the  premixed  flame  segment  can  be  stable  if  the  re¬ 
sult  (32)  exceeds  the  result  (24)  for  any  k,  co  and  r.  Figure  2  shows 
Eq.  (24),  dashed,  and  Eq.  (32),  solid,  with  0  =  1.2  and  ©  =  8  in 
Figs.  2(a)  and  (b),  respectively.  We  note  that  the  solid  line,  Eq.  (32), 
achieves  its  minimal  value  when  the  oscillation  frequency  is  in  res¬ 
onance  with  the  DL  parameter  of  frequency  dimension 

a)0  =  2(A@)V2SLk,  (34) 

which  is  related  to  the  characteristic  time  of  development  of  the 
DL  instability  excited  by  a  perturbation  of  the  wavenumber  k.  A 
“stability  band”  observed  for  small  thermal  expansion  in  Fig.  2a 
designates  absolute  stability  against  perturbations  of  any  wave 


Analysis  of  flame-acoustic  interaction  developed  in  Section  3  is 
limited  to  the  Landau  limit  of  zero  flame  thickness,  L  =  0.  It  is 
noted,  however,  that  accounting  for  the  internal  flame  structure 
leads  to  qualitatively  and  quantitatively  new  effects.  First  of  all, 
the  finite  flame  thickness  can  stabilize  the  DL  effects  at  sufficiently 
small  scales  [17,18].  Second,  for  non-unity  Lewis  number  Le,  the 
diffusional-thermal  instability  becomes  effective,  leading  to  cellu¬ 
lar  or  pulsating  behavior  of  Le  <  1  or  Le  >  1  premixtures,  respec¬ 
tively  [1].  Finally,  stretch  effects  are  activated  for  Le  ^  1  flames 
of  finite  thickness.  A  positive/negative  stretch  suppresses/promotes 
the  cellular  instability,  while  promoting/suppressing  the  pulsating 
instability  [1]. 

Within  the  approach  of  small,  but  finite  flame  thickness,  0  <  kl 
<1,  Ref.  [17]  extended  the  Darrieus-Landau  dispersion  relation  to 
a  =(x0(l  -k/kc),  where  cr0  is  the  DL  growth  rate  in  the  Landau 
limit,  determined  by  Eq.  (23)  without  acoustics,  U  =  0,  and  kc  the 
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DL  cut  off  wavenumber  given  by 
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(35) 


where  the  function  h($)  describes  the  temperature  dependence  of 
the  transport  coefficients,  with  h(l)  =  1,  h(©)  =  hb,  Pr  the  Prandtl 
number,  and  Mk  the  Markstein  number.  The  latter  describes  the 
flame  response  to  curvature  and  stretch,  and  it  is  related  to  the 
Zel’dovich  and  Lewis  numbers,  Ze  and  Le,  as  [17-19] 


Mk  = 


=  —  f 

©-1  J 
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It  is  noted  that  the  characteristic  size  of  the  burning  zone  can 
vary  with  the  locus,  L  =  L(r),  because  ©,  Y  and  SL  do.  Nevertheless, 
similar  to  all  these  flame  parameters,  we  assume  L  to  be  constant 
within  any  separate  quasi-planar  flame  segment. 

As  a  result,  we  readily  identify  a  preliminary  stability  criterion, 
which  also  specifies  the  limits  of  the  present  formulation.  This  cri¬ 
terion  simply  states  that  any  small  quasi-planar  premixed  segment 
of  constant  flame  parameters  should  nevertheless  exceed  the  DL 
cut-off  wavelength  Xc  =  2n  /kc.  Mathematically,  such  an  estimate 
takes  the  form 


J_d©  l_dY 

©  dr  <  <c’  Y  dr 


1  d§L  <  k 

SL7F<kc 


Ze  , 

J—  <  kc, 

^ mix 


(37) 


where  Lmix  is  a  characteristic  size  of  the  fuel-oxidizer  mixing  layer 
ahead  of  the  nozzle.  Of  course,  the  result  (37)  is  an  estimate,  while 
for  a  rigorous  stabilization  criterion  we  have  to  reproduce  the  anal¬ 
ysis  of  Section  3  accounting  for  the  flame  thickness. 

Within  the  approach  of  0  <  kL  <1,  extension  of  Eq.  (18)  takes 
the  form 


C,d2f/dt2+C2SLk(df/dt) 

where 


Ci  =  1  +AkL\Mk- 


C3S2/<2/  =  C41<U co  cos  (cot)f,  (38) 
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Following  Section  3,  we  derive  the  counterparts  of  inequalities 
(24)  and  (32), 
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Frequency:  co  /  2(A0)I/ZSL  k 

Fig.  3.  Stability  diagrams  for  zero,  L  =  0,  and  non-zero,  kL  =  0.1,  flame  thickness. 
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Fig.  4.  Stability  diagrams  for  a  flame  of  finite  thickness  with  thermal  expansion 
0  =  8  and  various  diffusivities:  Le  =  1,  Mk  =  4.2  (a),  and  Le  ±  1,  Mk  =  7  (b). 


The  results  (41)  and  (42)  are  shown  in  Fig.  3  for  kL  =  0.1  as 
compared  to  the  Landau  limit,  L  =  0,  shown  in  Fig.  2b,  with  ©  =  8 
in  both  events.  It  is  seen  in  Fig.  3  that,  unlike  the  Landau  limit, 
accounting  for  the  finite  flame  thickness  predicts  steady  LPF  and 
RPH  sheets  for  certain  conditions. 

Figure  4  presents  the  scaled  acoustic  amplitude  Ucosa/SL  ver¬ 
sus  the  scaled  wavenumber  kL  for  an  equidiffusive  premixure,  Le  = 


1,  with  Mk  =  4.2  (Fig.  4a),  and  a  non-equidiffusive  premixture,  Le  ± 
1,  with  Mk  =  1  (Fig.  4b),  with  the  expansion  factor  ©  =  8  and  sev¬ 
eral  fixed  scaled  frequencies  cdL/Sl  =  1;  2;  3  in  both  plots.  Again, 
we  observe  a  “stability  band”  for  a  non-equidiffusive  flame  in 
Fig.  4b,  but  there  is  no  such  a  band  in  Fig.  4a.  Consequently,  it 
is  seen  that  the  increase  in  the  flame  thickness  and  the  Markstein 
number  modifies  the  effect  of  the  combustion  instability. 
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Based  on  Figs.  3  and  4(a,  b),  we  therefore  conclude  that  a  fi¬ 
nite  flame  thickness  moderates  the  combustion  instability,  due  to 
the  thermal  stabilization  of  the  DL  instability  mode.  Furthermore, 
diffusional-thermal  effects  within  a  non-equidiffusive  burning  zone 
can  also  stabilize  the  flame. 

6.  Concluding  remarks 

In  summary,  we  anticipate  the  following  scenario  of  flame  sta¬ 
bilization/destabilization  of  a  triple  flame.  First,  hydrodynamic  in¬ 
stability  increases  the  total  burning  rate,  hence  generating  sound 
waves,  which  in  turn  can  suppress  the  instability.  Second,  the 
sound  wave  can  trigger  parametric  instability.  Therefore,  in  order 
to  retain  a  lifted  flame  stable  we  have  to  keep  the  flame,  setup 
and  acoustical  characteristics  at  any  point  of  the  front  within  the 
stability  range  shown  in  Figs.  2-4. 

Finally,  it  is  noted  that  there  are  two  “natures”  of  acoustics  and 
flame-acoustic  coupling:  the  pressure-wise  and  the  velocity-wise 
ones.  The  pressure-wise  nature  is  associated  with  the  acoustic  en¬ 
ergy,  with  the  stability  determined  by  the  classical  Rayleigh  cri¬ 
terion.  The  acoustic  energy  and  Rayleigh  criterion  are  beyond  the 
primary  interest  of  this  paper.  However,  it  may  oftentimes  happen 
that  while  the  Rayleigh  criterion  predicts  the  stability,  is  does  not 
occur  in  the  practical  reality  because  of  the  instability  related  to 
the  velocity-wise  nature  of  the  flame-acoustic  coupling.  It  is  this 
situation  that  constitutes  the  focus  of  the  present  work. 
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Cool  flames,  being  essential  features  of  chemical  kinetics  of  large  hydrocarbon  fuels,  are  closely  related  to 
the  negative  temperature  coefficient  (NTC)  phenomenon  and  engine  knock.  In  this  work,  the  coupling  of 
cool  flame  chemistry  and  convective-diffusive  transport  is  computationally  and  experimentally  investi¬ 
gated.  A  1-D  planar  premixed  cool  flame  induced  by  a  hot  pocket  is  first  simulated  for  DME/02/N2  mix¬ 
tures  with  detailed  chemistry  and  transport,  demonstrating  the  existence  of  a  residence  time  window  for 
quasi-steady  propagation.  Then  with  residence  time  limited  by  aerodynamic  straining,  a  steady-state 
premixed  cool  flame  is  simulated  in  a  counterflow  of  heated  N2  against  a  DME/02/N2  mixture.  It  is  found 
that  with  a  high  strain  rate,  corresponding  to  short  residence  time,  low-temperature  heat  release  is  sup¬ 
pressed,  resulting  in  a  stretched  low-temperature  S-curve  system  response;  and  that  with  a  sufficiently 
low  strain  rate,  corresponding  to  long  residence  time,  ignition  induced  by  low-temperature  chemistry 
would  transition  to  a  high-temperature,  intensely  burning  flame.  Consequently,  a  steady-state  premixed 
cool  flame  exists  only  for  residence  time  in  a  strain  rate  window.  A  symmetric  counterflow  configuration 
is  then  simulated  to  determine  the  cool  flame  temperature  and  flame  speed  at  a  fixed  local  strain  rate, 
showing  very  different  controlling  chemistry  and  characteristics  as  compared  to  the  normal  laminar 
flames  governed  by  high-temperature  chemistry.  In  a  companion  experimental  investigation,  premixed 
cool  flames  in  the  counterflow  were  observed  with  a  high-sensitivity  CCD  camera  in  the  UV  spectrum, 
with/without  a  bandpass  filter  corresponding  to  the  characteristic  wavelength  of  excited  HCHO.  The 
chemiluminescence  from  the  cool  flame  is  found  to  become  more  intense  with  increasing  equivalence 
ratio,  even  for  rich  mixtures,  while  the  position  of  the  cool  flame  is  insensitive  to  variation  in  the  equiv¬ 
alence  ratio  at  the  same  strain  rate.  These  observations  qualitatively  agree  with  the  numerical  simula¬ 
tions,  demonstrating  the  essential  features  of  premixed  cool  flames. 

©  2015  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 

Cool  flames,  discovered  in  1817  [1],  have  been  studied  exten¬ 
sively  [2].  Its  controlling  kinetics  of  the  low-temperature  chemistry 
[3]  has  been  shown  to  be  closely  related  to  the  phenomena  of  two- 
stage  ignition  and  the  negative  temperature  coefficient  (NTC)  [4], 
associated  with  the  autoignition  [5]  and  engine  knock  of  large 
hydrocarbon  fuels  [6].  They  are  normally  characterized  by  faint 
bluish  luminescence  from  excited  formaldehyde  [7],  relatively 
small  heat  release  and  very  low  C02  production.  Furthermore,  con¬ 
trary  to  the  normal  hot  flames  in  which  chemical  reactions  and 
transport  are  strongly  coupled,  cool  flames  have  been  mostly 
observed  in  homogeneous  systems  such  as  heated  closed  reaction 
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vessels  [8]  and  jet-stirred  reactors  [9],  exhibiting  oscillatory  volu¬ 
metric  heat  release  and  chemiluminescence. 

Noting  that  gas-phase  nonuniformities  invariably  exist  in  prac¬ 
tical  systems,  recently  there  has  been  considerable  interest  in 
understanding  the  coupled  effects  of  cool  flame  chemistry  and 
transport,  leading  to  the  postulation  and  identification  of  this  class 
of  laminar  flames  whose  flame  temperature  can  be  substantially 
lower  than  the  adiabatic  flame  temperature  corresponding  to  the 
(near-)complete  consumption  of  the  controlling  reactant(s).  In  par¬ 
ticular,  Law  and  Zhao  [10]  and  Zhao  and  Law  11]  demonstrated 
computationally  that  with  low  strain  rates  and/or  high  pressures, 
low-temperature  chemistry  can  strongly  couple  with  transport  to 
induce  a  secondary  S-curve  on  the  lower  branch  of  the  conven¬ 
tional  (primary)  S-curve,  with  its  own  distinct  ignition  and  extinc¬ 
tion  states.  Subsequently,  using  infrared  detection  Deng  et  al.  [12] 
experimentally  observed  a  diffusion  cool  flame  in  the  nonpremixed 
counterflow  at  atmospheric  pressure.  Furthermore,  diffusion  cool 
flames  have  been  suggested  to  explain  some  experimental  results 
on  microgravity  droplet  combustion  [13,14],  and  have  also  been 
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directly  observed  in  experiments  on  plasma-assisted  combustion 

[15]. 

Recognizing  that  recent  investigations  [10-15]  on  the  coupled 
behavior  of  low-temperature  chemistry  and  convective-diffusive 
transport  have  been  mainly  on  diffusion  flames,  it  behooves  us  to 
explore  if  steady,  premixed  cool  flames  can  also  be  initiated  and 
established,  and  if  so,  what  are  the  associated  flame  structure 
and  propagation  rates,  including  the  role  of  the  two-stage  ignition 
delay  in  establishing  such  flames.  The  relevant  issue  here  is  that 
while  a  diffusion  flame  is  primarily  controlled  by  the  counter¬ 
diffusion  of  the  fuel  and  oxidizer  reactants,  with  chemistry  playing 
a  secondary  role  in  the  global  flame  characteristics,  the  structure 
and  propagation  of  premixed  flames  inherently  depend  on  the  cou¬ 
pled  influence  of  chemical  kinetics  and  transport  [16].  Conse¬ 
quently  an  assessment  of  the  essential  role  of  low-temperature 
chemistry  on  flames  must  necessarily  involve  studies  on  premixed 
flames.  In  fact,  because  of  the  very  weak  exothermicity  of  the  low- 
temperature  chemistry,  and  the  inevitable  heat  loss  in  a  flame  sys¬ 
tem,  it  is  not  evident  a  priori  that  a  flame  embodying  this  chemistry 
can  be  readily  initiated  and  support  steady  propagation.  It  is  also 
noted  that  while  premixed  cool  flames  have  been  studied  previ¬ 
ously,  for  example  as  early  as  the  1960s  [17,18],  the  increased 
understanding  of  flame  structures  since  then  is  expected  to  yield 
further  insight  and  interpretation  of  this  important  phenomenon. 

The  present  study  consists  of  three  tasks.  We  first  extend,  in  the 
next  section,  the  classical  NTC  phenomenon  describing  the  tran¬ 
sient  ignition  of  a  homogeneous  mixture  to  one  in  which  ignition 
is  accomplished  by  a  localized  heat  source  in  the  1-D  planar 
domain,  which  is  inherently  nonhomogeneous  in  nature  and  is 
affected  by  both  chemistry  and  diffusive  transport.  Consequently 
the  interest  here  is  the  dynamics  and  chemistry  of  the  ignition  pro¬ 
cess  leading  to  low-temperature  flame  propagation  and  its  subse¬ 
quent  transition  to  the  high-temperature  flame.  We  shall  also 
show  that  the  transient  ignition  event  depends  sensitively  on  the 
temperature  of  the  ignition  source  when  low-temperature  chem¬ 
istry  is  involved. 

In  order  to  capture  and  stabilize  the  low-temperature  flame 
identified  in  Section  2,  in  Section  3  we  computationally  explore 
the  possible  existence  of  a  premixed  flame  of  low-temperature 
chemistry  in  the  counterflow,  with  its  characteristic  residence  time 
controlled  by  the  strain  rate  of  the  flow.  We  shall  show  that  such  a 
flame  is  predicted  to  exist,  with  its  chemical  structure  including 
the  states  of  ignition  and  extinction  depending  sensitively  on  the 
strain  rate  of  the  flow.  Our  third  task,  presented  in  Section  4,  pro¬ 
vides  the  first  experimental  observation  that  such  stretch- 
stabilized  premixed  cool  flames  indeed  exist,  thereby  brings  more 
insights  on  the  existence  and  characteristics  of  laminar  cool  flames, 
both  nonpremixed  [10-15]  and  premixed  [17-20]. 

The  fuel  selected  for  the  study  is  dimethyl  ether  (DME),  which  is 
one  of  the  smallest  fuel  molecules  exhibiting  the  NTC  behavior.  The 
reaction  mechanism  adopted  in  the  computation  is  a  validated 
reduced  model  [21]  consisting  of  both  low-  and  high- 
temperature  kinetics;  the  associated  validation  and  applications 
are  reported  in  [12,22]. 


2.  Initiation  and  propagation  of  the  1-D  planar  cool  flame 

We  first  examine  the  adiabatic  and  isobaric  autoignition  delays 
of  homogeneous  stoichiometric  DME/02/N2  mixtures,  in  order  to 
select  the  proper  conditions  to  study  the  ignition  and  propagation 
of  1-D  planar  flames.  Fig.  1(a)  shows  the  calculated  first-stage, 
second-stage,  and  total  ignition  delay  times  for  initial  temperature 
T,  denoted  as  Tlf  t2  and  t  respectively.  It  is  seen  that  t  varies  non- 
monotonically  with  increasing  T,  and  hence  demonstrates  the  NTC 
behavior.  Consequently,  for  an  appropriate  total  ignition  delay 
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Fig.  1.  (a)  First-stage,  second-stage,  and  total  ignition  delays  of  DME/02/N2  mixture 
with  15.8%  N2  in  volume,  under  1  atm.  The  inset  shows  the  definitions  of  various 
ignition  delays  in  a  two-stage  ignition  process,  (b)  Heat  release  rate  history  for 
ignition  processes  at  643  K,  750  K  and  967  K. 


time,  say  t  =  22.9  ms,  there  are  three  distinct  ignition  situations: 
at  the  low  temperature  state  of  T  =  643  K,  Ti  is  much  longer  than 
t2;  at  the  intermediate  temperature  state  of  T  =  750  K,  t2  is  much 
longer  than  ti;  and  finally  at  the  high  temperature  state  of 
T  =  967  K,  the  two-stage  ignition  behavior  disappears  and  combus¬ 
tion  is  totally  taken  over  by  the  intermediate-to-high  temperature 
chemistry,  exhibiting  a  single-stage  ignition  event.  These  three  dis¬ 
tinct  ignition  situations  are  further  characterized  by  their  corre¬ 
spondingly  distinct  heat  release  rates,  Q,  as  shown  in  Fig.  1(b). 
Specifically,  for  the  low  and  intermediate  temperatures  of  643  K 
and  750  K,  heat  release  occurs  in  two  bursts,  corresponding  to 
the  first  and  second  stage  ignition.  Furthermore,  the  first  ignition 
delay  for  the  750  K  case  occurs  more  readily  than  that  of  the 
643  K  case,  as  expected.  For  the  high-temperature  case  of  967  K, 
there  is  no  first-stage  ignition  and  as  such  ignition  occurs  in  a  sin¬ 
gle  burst,  through  the  second  ignition  delay  chemistry.  It  is  also 
clear  that  Q  is  much  lower  both  before  and  after  the  respective 
ignition  bursts. 

Based  on  these  observations,  we  explore  how  NTC  would  affect 
the  initiation  and  propagation  of  a  1-D  planar  flame  induced  by 
an  ignition  kernel  with  these  three  different  temperatures.  The 
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problem  studied  is  as  follows.  At  the  initial  state,  a  quiescent, 
homogeneous,  stoichiometric  DME/02/N2  mixture  with  15.8%  N2 
in  volume,  and  at  300  K  and  1  atm,  fills  a  5  cm  computational 
domain.  Then,  a  hot  pocket  of  initial  temperature  643  K,  750  K  or 
967  K,  identified  above,  is  respectively  imposed  from  0  to  1  cm  in 
the  domain  to  initialize  flame  propagation.  Due  to  the  heat  deposi¬ 
tion,  the  kernel  gas  not  only  tends  to  auto-ignite,  it  also  conducts 
heat  into  the  unreacted  mixture  causing  it  to  ignite.  It  is  therefore 
expected  that  if  the  thermal  diffusion  time  scale  is  longer  than 
the  first-stage  autoignition  delay  time  Ti  of  the  gas  in  the  ignition 
kernel,  low-temperature  chemistry  will  be  rapid  enough  to  heat 
up  the  kernel  gas  to  the  cool  flame  temperature.  The  time  evolu¬ 
tions  starting  from  the  initial  condition  will  be  shown  in  the  follow¬ 
ing,  under  the  ideal,  adiabatic  condition. 

All  simulations  were  carried  out  using  the  Adaptive  Simulation 
of  Unsteady  Reacting  Flows  (A-SURF)  code  [23,24],  which  has  been 
successfully  used  in  various  investigations  of  unsteady  flame  initi¬ 
ation  and  propagation  [23-26].  It  simulates  the  1-D  compressible 
multicomponent  reacting  flow  in  the  planar  coordinate  using  the 
finite  volume  method.  Details  of  the  governing  equations,  numer¬ 
ical  method  and  code  validation  are  given  in  [23,24].  A  multi¬ 
level  adaptive  mesh  with  a  minimum  mesh  size  of  10  pm  is  used 
to  capture  the  flame  structure  and  propagation.  Zero-gradient 
boundary  condition  is  imposed  at  the  left  boundary  (x  =  0)  and 
transmissive  condition  at  the  right  boundary  (x  =  5  cm).  To  avoid 
discontinuity  and  the  associated  numerical  difficulties,  we  used 
the  hyperbolic  tangent  function,  tank,  to  denote  the  step  function 
in  the  initial  temperature  profile.  In  order  to  show  evolution  of 
the  reaction  front  and  the  detailed  thermal  structure,  the  T  and  Q. 
profiles  at  the  representative  state  for  each  case  are  plotted  in 
Figs.  2-4.  For  clarity  in  discussion,  we  designate  the 
high-temperature  reaction  front  with  its  temperature  around  the 
adiabatic  flame  temperature  as  the  hot  flame  front  (HFF),  and  a 
low-temperature  reaction  front  with  its  temperature  below  or 
around  1000  K  as  the  cool  flame  front  (CFF);  the  reaction  fronts 
are  indicated  by  dashed  arrows. 

Fig.  2  shows  evolution  of  the  temperature  and  heat  release  rate 
profiles  for  the  low-temperature  case  of  T  =  643  K.  It  is  seen  that 
first-stage  autoignition  is  initiated  around  the  center  region  of 
the  kernel  (curve  2),  which  drives  up  the  center  T  to  around 
900  K,  with  the  corresponding  Q.  occurring  in  the  bulk  of  the  kernel 
region.  Shortly  thereafter,  the  Q  in  the  kernel  region  decreases  as 
the  first-stage  reactions  are  completed  (left  curve  3),  while  at  the 
same  time  a  low-temperature  reaction  front  with  a  distinct  peak 


in  Q  starts  to  propagate  into  the  unburnt  mixture  at  x  =  1 .4  cm 
(right  curve  3),  which  is  the  CFF  initiated  by  the  low-temperature 
kernel.  With  this  CFF  further  propagating  to  x  =  1.5  cm  (right  curve 
4),  another  high-temperature  front,  an  HFF  segment  with  much 
higher  Q.  (left  curve  4)  is  induced  around  x  =  0.9  cm  by  the 
second-stage  autoignition,  forming  an  extended  reaction  structure 
consisting  of  distinct  HFF  and  CFF  segments.  Subsequently,  as 
autoignition  ends,  diffusion  becomes  important  and  gradually 
smooths  out  this  combined  front  to  approach  that  of  the  normal 
high-temperature  premixed  flame,  shown  as  curve  5.  It  is  noted 
that  the  instantaneous  Q  for  the  HFF  in  the  (left)  curve  4  is  much 
higher  than  that  of  curve  5,  because  autoignition  occurs  much  fas¬ 
ter  than  diffusion,  and  therefore  heat  and  species  generated  from 
autoignition  cannot  readily  diffuse  away  as  autoignition  takes 
place. 

For  the  high-temperature  ignition  of  T  =  967  K,  which  is  beyond 
the  NTC  regime,  Fig.  3  shows  that  while  the  two-stage  autoignition 
does  not  occur  for  this  initial  temperature  as  shown  in  Fig.  1,  the 
initial  temperature  stratification,  spanning  from  300  I<  to  967  K 


Fig.  3.  Temperature  and  heat  release  profiles  at  various  time  steps  for  a  1-D  cool 
flame  initiated  by  a  hot  pocket  of  967  K,  under  1  atm.  [1]:  0  ms,  [2]:  20  ms,  [3]: 
27  ms,  [4]:  29.625  ms,  [5]:  29.75  ms,  [6]:  29.875  ms. 
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Fig.  2.  Temperature  and  heat  release  profiles  at  various  time  steps  for  a  1-D  cool 
flame  initiated  by  a  hot  pocket  of  643  K,  under  1  atm.  [1]:  0  ms,  [2]:  16.75  ms,  [3]: 

22  ms,  [4]:  22.236  ms,  [5]:  22.286  ms.  The  dash  arrows  indicate  the  propagation  of 
the  cool  and  hot  flame  front  segments,  as  CFF  and  HFF. 
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Fig.  4.  Temperature  and  heat  release  profiles  at  various  time  steps  for  a  1-D  cool 
flame  initiated  by  a  hot  pocket  of  750  K,  under  1  atm.  [1]:  0  ms,  [2]:  7  ms,  [3]: 
13  ms,  [4]:  23.736  ms,  [5]:  23.8  ms.  The  zoom-in  plot  shows  the  thermal  structure 
of  propagating  cool  flames  1,  2  and  3. 
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around  x  =  1  cm,  shown  as  curve  1,  implies  that  there  still  exists  a 
low-temperature  region  favoring  NTC-affected  ignition,  and  conse¬ 
quently  a  CFF  can  still  be  initiated,  as  shown  by  curves  2  and  3. 
When  time  approaches  29  ms,  while  the  CFF  still  propagates 
towards  the  unburned  mixture,  the  main  ignition  event  takes  place 
at  the  mixing  layer  around  x  =  1  cm,  leading  to  the  formation  of 
two  HFFs  propagating  left  and  right  respectively,  as  shown  in 
curves  4  and  5.  The  right-propagating  reacting  front  with  the  com¬ 
bined  CFF  and  HFF  segments  (right  branches  of  curves  4  and  5) 
again  is  smoothed  out  by  diffusion  and  eventually  approaches 
the  normal  hot  premixed  flame  (right  curve  6). 

Fig.  4  shows  evolution  of  the  T  and  Q.  profiles  for  the  intermedi¬ 
ate  temperature  case  (T=750I<).  Similar  to  the  low-  and  high- 
temperature  cases,  it  is  seen  that  subsequent  to  the  first-stage  igni¬ 
tion,  a  CFF  forms  and  propagates  into  the  unburned  mixture,  as 
seen  in  curves  2  and  3  in  Fig.  4.  Although  this  front  is  sustained 
much  longer  due  to  the  prolonged  second-stage  ignition  delay  t2, 
it  still  lasts  for  only  a  limited  residence  time.  After  23  ms,  an  HFF 
occurs  around  x  =  1 .4  cm,  as  curve  4,  and  merges  with  the  CFF  part. 
Such  a  combined  reacting  front  then  gradually  transitions  to  a 
high-temperature  premixed  flame,  as  curve  5. 

These  evolution  processes  show  that  a  premixed  cool  flame  is  a 
thin  thermal-chemical  layer  bridging  the  upstream  ambient  and 
downstream  low-temperature  quasi-equilibrium  induced  by  low- 
temperature  ignition,  as  the  CFFs.  Such  a  flame  front  with  the 
inherent  coupling  of  low-temperature  chemistry  and  transport 
could  generally  be  initiated  by  a  hot  pocket  with  heat  deposition, 
regardless  of  the  extent  of  the  thermal  stratification  induced  by 
different  kernel  temperatures.  We  have  therefore  demonstrated 
the  general  relevance  of  cool  flames  sustained  by  temperature  gra¬ 
dients,  even  for  hot  pockets  with  temperatures  beyond  the  NTC 
regime,  as  also  noted  in  [27].  Although  these  cases  have  different 
durations  for  the  self-sustained  CFF  propagation,  eventually  they 
all  transition  to  the  regular  hot  flame  with  the  corresponding  adi¬ 
abatic  flame  temperature.  In  fact,  by  comparing  Figs.  2  and  4,  it  is 
clear  that  kernels  within  the  NTC  regime  with  longer  t2  can  induce 
a  longer  self-sustained  cool  flame  propagation  period. 

To  further  investigate  the  characteristics  of  cool  flame  propaga¬ 
tion  for  the  intermediate  temperature  case,  evolutions  of  the  max¬ 
imum  temperature  and  flame  position  are  shown  in  Fig.  5.  It  is  seen 
that  the  maximum  temperature  first  stays  around  750  K  before  the 
first-stage  autoignition  takes  place.  At  about  3  ms,  a  cool  flame 
appears  and  starts  to  propagate  until  23  ms,  at  which  it  transitions 
to  a  hot  flame  with  a  much  higher  temperature,  hence  further 
demonstrating  the  existence  of  a  much  wider  residence  time  win¬ 
dow  of  20  ms  for  the  presence  of  the  cool  flame.  A  linear  fitting 
between  the  flame  position  and  the  time  during  this  stage  yields 
the  flame  speed  of  the  burnt  mixture,  Sbf  as  23.2  cm/s.  Using  the 
thermal  expansion  ratio  a  =  3.5  evaluated  from  the  unburned  and 
burned  density  ratio,  the  unburned  flame  speed  is  then  Su  =  Sbl 
a  =  6.6  cm/s;  the  magnitude  is  consistent  with  that  determined 
from  the  counterflow  cool  flame  simulation  shown  in  Fig.  14(a) 
in  the  next  section. 

Comparison  with  the  regular  flame  has  been  carried  out  by  con¬ 
ducting  a  simple  laminar  flame  calculation  for  stoichiometric  DME/ 
02/N2  with  15.8%  N2  under  300  K  and  1  atm,  the  same  condition  as 
the  cool  flame  shown  in  Fig.  5.  It  is  seen  that  the  corresponding 
unburnt  flame  speed  Su  for  a  regular  hot  flame  is  about 
246.9  cm/s,  with  an  adiabatic  flame  temperature  2976  K  corre¬ 
sponding  to  a  thermal  expansion  ratio  of  10.0,  and  a  thin  flame 
thickness  of  about  7.9E-03  cm.  However,  the  cool  flame  thickness 
similarly  obtained  by  taking  the  ratio  of  the  temperature  difference 
on  the  burned  and  unburned  side  with  the  maximum  temperature 
gradient  in  the  flame  structure,  is  about  0.137  cm,  together  with  its 
much  lower  flame  temperature  of  901  K  and  Su  of  6.6  cm/s.  From 
the  classical  thermal  theory  of  flames,  the  product  of  the  unburnt 
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Fig.  5.  (a)  maximum  temperature,  flame  position  and  (b)  flame  speed  evolution  for 
a  1-D  cool  flame  initiated  by  a  hot  pocket  of  750  K  (solid)  and  643  K  (dash),  under 
1  atm. 


flame  speed  Su  with  the  flame  thickness  S  should  scale  with  the 
mass  diffusivity  D,  i.e.,  SUS  ~  D.  The  calculated  product  SUS  for  the 
cool  flame  is  0.9  cm2/s,  while  this  product  for  the  hot  flame  is  about 
1.9cm2/s.  Considering  their  largely  different  chemical  reactivity 
and  speciation,  a  difference  by  a  factor  of  2  from  the  flame  thermal 
theory  implies  similar  transport  properties  of  the  cool  and  regular 
flames,  and  the  physical  consistency  of  the  cool  flame  speed. 

Having  determined  the  propagation  speed  for  the  flame  identi¬ 
fied  above,  it  is  important  to  recognize  that  a  “cool  flame  speed”  is 
fundamentally  not  a  well-defined  concept  as  compared  to  the  well- 
defined  laminar  flame  speed,  even  in  the  1-D  planar  configuration. 
Specifically,  the  downstream  of  a  cool  flame  is  only  one  of  quasi 
chemical  equilibrium  in  a  state  of  transition  such  that,  with  suffi¬ 
ciently  long  residence  time,  it  will  eventually  evolve  to  the  final 
equilibrium  state  characterized  by  the  adiabatic  flame  temperature 
of  the  mixture.  As  such,  a  steady-state  1-D  premixed  cool  flame 
with  an  attendant  propagation  speed  in  the  doubly-infinite 
domain,  where  infinite  residence  time  is  allowed  for  chemical 
equilibrium  downstream  and  chemical  frozen  condition  upstream, 
is  not  conceptually  viable  and  as  such  does  not  exist. 

Having  recognized  the  failure  of  a  1-D  cool  flame  to  assume  an 
eigen-mode  of  propagation  in  the  doubly  infinite  regime,  there  are 
at  least  two  scenarios  in  which  a  steadily  propagating,  globally- 
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planar  cool  flame  can  be  envisioned.  The  first  is  the  freezing  of  the 
downstream  chemistry  through  radiation  loss,  rendering  it  to  be 
asymptotically  non-reactive.  The  problem  is  not  unlike  that  of 
the  cold-boundary  difficulty  in  the  conventional  premixed  flame 
analysis,  in  that  given  infinite  time  any  reaction  with  a  finite  tem¬ 
perature  will  be  able  to  completely  react;  in  the  far  upstream  for 
the  conventional  flame  and  far  downstream  for  the  cool  flame, 
even  with  radiative  loss.  Thus  the  eigenvalue  nature  of  the  problem 
can  only  be  defined  in  some  asymptotic  sense. 

The  second  scenario  to  establish  a  steady  cool  flame  is  to  trun¬ 
cate  the  downstream  flow  time  so  as  to  prevent  transition  of  the 
cool  flame  chemistry  to  reach  the  state  of  chemical  equilibrium. 
The  downstream  state  of  this  flame  is  then  in  quasi-equilibrium 
with  the  cool  flame  chemistry.  This  can  be  readily  accomplished 
by  aerodynamically  straining  the  flow,  as  in  the  case  of  a  stagna¬ 
tion  flow  or  counterflow.  We  shall  therefore  study  in  the  following 
stretched  cool  flames  in  the  counterflow,  both  computationally  and 
experimentally.  Indeed,  since  the  laminar  flamelets  within  a  turbu¬ 
lent  flame  structure  are  mostly  positively  stretched,  the  results 
bear  relevance  to  practical  situations. 

It  is  also  worth  noting  that  the  artificial  air  we  adopted  with  sig¬ 
nificant  N2  reduction  is  not  a  necessary  condition  for  propagating 
cool  flames,  which  instead  largely  depends  on  the  thermodynamic 
conditions.  A  propagating  cool  flame  can  also  be  observed  in  regu¬ 
lar  air  if  the  pressure  is  high  enough,  as  in  the  environments  of 
rapid  compression  machine  and  engine  experiments.  As  shown  in 
the  next  section,  ignition  induced  by  low  temperature  chemistry 
becomes  much  more  pronounced  if  the  dilution  effect  from  nitro¬ 
gen  is  reduced.  The  reason  is  that  with  the  same  low  temperature 
heat  release,  less  inert  content  leads  to  higher  increment  of  the 
cool  flame  temperature. 


3.  Stretched  cool  flame  in  counterflow  -  computational  results 

In  order  to  explore  if  a  steady  cool  flame  with  truncated  down¬ 
stream  flow  time  can  be  physically  realized,  a  premixed  counter¬ 
flow  configuration  consisting  of  a  hot  N2  stream  flowing  against 
a  premixed  DME/02/N2  stream  at  room  temperature  was  studied, 
based  on  the  following  considerations.  First,  the  flow  along  the 
center  line  is  1-D  and  the  residence  time  is  well  characterized  by 
the  aerodynamic  strain.  Furthermore,  by  fixing  the  strain  rate 
and  varying  the  inert  boundary  temperature,  the  system  could 
show  multiple  solutions  corresponding  to  those  defined  by  the 
characteristic  S-shaped  response  curves,  with  the  lower  and  upper 
turning  points  respectively  designate  the  ignition  and  extinction 
states  of  the  system.  It  is  therefore  of  interest  to  explore  if  a  pre¬ 
mixed  cool  flame  can  be  stabilized  and  as  such  also  exhibit  distinct 
ignition  and  extinction  states,  in  manners  similar  to  those  of  flames 
in  practical  situations  with  finite  local  residence  times. 

The  governing  equations  for  the  counterflow  flame  simulation 
are  given  in  [28].  The  simulation  employs  the  damped  Newton’s 
method  and  time  integration  solution  scheme  to  solve  the  ordinary 
differential  equations  with  boundary  conditions  specified  on  both 
sides  of  the  potential  flow.  S-curve  marching  is  performed  using 
the  flame-controlling  method  of  Nishioka  et  al.  [29],  with  detailed 
chemistry  and  transport  database.  The  critical  states  of  ignition 
and  extinction  can  then  be  assessed  through  the  multiple  turning 
points  of  the  S-curve.  The  simulation  conditions  are  as  follows: 
the  upper  boundary  consists  of  a  mixture  of  DME,  oxygen  and 
nitrogen  at  a  fixed  temperature  of  300  K,  while  the  lower  boundary 
is  nitrogen  at  a  higher  temperature.  The  separation  distance 
between  the  nozzle  exits  is  2  cm.  To  generate  an  S-curve,  marching 
is  performed  by  changing  the  nitrogen-side  boundary  temperature 
with  fixed  strain  rate. 


Fig.  6  shows  the  S-curve  for  a  stoichiometric  DME/02  mixture 
with  15.8%  N2  at  different  strain  rates,  using  the  maximum  CH3 
concentration  as  a  representative  global  response  to  the  nitrogen 
boundary  temperature.  Due  to  the  fact  that  the  temperature  varia¬ 
tion  during  low-temperature  ignition  shows  an  S-curve  response, 
all  radicals  show  corresponding  S-curve  responses  with  the  same 
turning  point  behavior,  as  has  also  been  shown  in  Refs.  [10,11]. 
Here  CH3  is  chosen  as  the  system  indicator  because  it  can  also  indi¬ 
cate  high-temperature  chemistry  while  the  low-temperature  spe¬ 
cies  (e.g.,  CH30CH202)  cannot.  It  is  seen  that,  for  each  of  the  two 
lower  strain  rates  of  k  =  20/s  and  50/s,  two  distinctive  ignition 
turning  points,  II  and  12,  appear,  with  II  controlled  by  the  low- 
temperature  chemistry,  and  12  by  the  intermediate-  to  high- 
temperature  chemistry.  Furthermore,  when  the  strain  rate  is 
increased  to  100/s,  the  S-curve  induced  by  the  low-temperature 
chemistry  loses  its  non-monotonicity  and  becomes  stretched,  due 
to  the  suppression  of  low-temperature  heat  release  and  thereby 
thermal  feedback  at  reduced  residence  time  [11].  It  is  also  noted 
that,  for  k  =  20/s,  the  system  directly  transitions  to  the  hot  flame 
upon  attaining  the  state  of  the  first-stage  ignition,  in  that  12  is 
not  accessed  with  increasing  boundary  temperature,  and  as  such 
demonstrates  the  significance  of  low-temperature  chemistry  at 
low  strain  rates.  This  is  to  be  contrasted  with  the  case  of  k  =  50/s, 
in  that  the  system  can  actually  maintain  a  state  of  steady  burning 
characterized  by  the  low-temperature  chemistry  after  attaining 
the  first  stage  ignition.  Intense  burning  is  attained  only  when  the 
state  of  12  is  reached  with  further  increase  in  temperature. 

We  next  investigate  the  transition  to  ignition  induced  by  the 
low-temperature  chemistry.  Based  on  the  steady-state  S-curve, 
the  transition  behavior  is  studied  by  including  the  unsteady  term 
and  introducing  a  1%  perturbation  to  the  local  temperature  profile 
of  the  steady-state  solution  at  the  first  ignition  turning  point.  Fig.  7 
shows  the  transition  process  of  the  low-temperature  chemistry 
induced  ignition  II,  for  k  =  50/s.  It  is  seen  that  after  introducing 
the  perturbation  at  t  =  0,  the  system  maximum  temperature  first 
drops  due  to  thermal  diffusion.  Then  both  the  maximum  tempera¬ 
ture  and  maximum  formaldehyde  mole  fraction  suddenly  increase 
around  t  =  0.2  s  and  eventually  approach  the  steady-state  cool 
flame  solution  on  the  upper  branch  of  the  low-temperature  S- 
curve,  with  the  maximum  temperature  in  the  flow  field  being 
about  100K  higher  than  the  boundary  temperature,  suggesting 
weak  combustion  and  heat  release  from  the  low-temperature 


Fig.  6.  S-curve  response  of  stoichiometric  DME/02  mixture  with  15.8%  N2  vs.  hot  N2 
at  different  strain  rates  20/s,  50/s  and  100/s,  showing  the  different  behavior  of  low 
temperature  ignition. 
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chemistry.  This  therefore  demonstrates  that  a  steady-state  cool 
flame  can  be  established  with  a  residence  time  corresponding  to 
an  appropriate  aerodynamic  strain  rate.  However,  for  the  even 
lower  strain  rate  of  k  =  20/s,  shown  in  Fig.  8,  the  system  evolves 
towards  a  hot  flame  around  t  =  0.4  s.  Such  a  transition  process 
shows  that,  with  long  enough  flow  residence  time  at  low  strain 
rate,  ignition  induced  by  low-temperature  chemistry  can  cause 
autoignition  and  eventually  achieve  intense  burning,  consistent 
with  the  steady-state  S-curve  interpretation  of  Fig.  6. 

Figs.  7  and  8  together  then  show  the  different  transition  pro¬ 
cesses  of  low-temperature  chemistry  induced  ignition  under  dif¬ 
ferent  strain  rates,  and  further  support  the  strain  rate  window 
identified  for  the  appearance  of  a  steady-state  cool  flame  shown 
in  Fig.  6,  with  longer  residence  time  promoting  autoignition  and 
shorter  residence  time  suppressing  the  thermal  feedback.  As  such, 
we  have  demonstrated  that  the  occurrence  of  steady  premixed  cool 
flames  strongly  depends  on  the  flow  residence  time. 

We  next  conduct  a  parametric  study  to  further  show  the 
response  of  flames  induced  by  the  low-temperature  chemistry  by 
varying  the  amount  of  dilution  and  equivalence  ratio.  For  the  case 
of  50/s  strain  rate  within  the  cool  flame  strain  rate  window,  effects 
of  dilution  are  observed  by  fixing  the  equivalence  ratio  and  varying 
the  amount  of  the  inert,  as  shown  in  Fig.  9.  It  is  seen  that  with 
decreasing  dilution,  the  S-curve  becomes  more  pronounced.  This 
is  due  to  the  fact  that,  with  equivalence  ratio  and  therefore  heat 
release  rate  of  the  low-temperature  chemistry  fixed,  the  tempera¬ 
ture  increases  with  reduced  dilution.  Consequently,  the  nonlinear 
thermal  feedback  to  the  S-curve  phenomena  becomes  stronger. 
On  the  other  hand,  with  dilution  increased  to  a  certain  amount, 
the  S-curve  induced  by  the  low-temperature  chemistry  becomes 
stretched  and  the  system  response  becomes  monotonic  without 
ignition  and  extinction. 

To  identify  the  dependence  of  low-temperature  ignition  on  the 
equivalence  ratio,  mixtures  with  different  equivalence  ratios  but 
the  same  amount  of  N2  dilution  are  simulated  for  50/s  and  1  atm, 
shown  in  Fig.  10.  It  is  seen  that  the  ignition  turning  point  is  not 
sensitive  to  the  equivalence  ratio,  yielding  fairly  close  ignition  tem¬ 
peratures.  However,  with  increasing  equivalence  ratio,  the  maxi¬ 
mum  CH3  mole  fraction  gradually  increases,  indicating  stronger 
reactivity  for  the  low-temperature  chemistry  under  richer 
conditions. 

Having  identified  the  ignition  characteristics  in  the  asymmetric 
counterflow,  we  next  study  the  thermal  and  chemical  structures  of 
the  cool  flame  in  the  symmetric  twin-flame  counterflow,  using 
OPPDIF  [30]  with  the  same  kinetic  model  and  transport  properties. 
The  nozzles  are  again  separated  by  2  cm.  A  plug  flow  boundary 


Time  (s) 

Fig.  7.  Transition  from  low  temperature  chemistry  induced  ignition  to  a  steady  cool 
flame,  for  stoichiometric  DME/02  with  15.8%  N2,  1  atm,  l<  =  50/s. 


Fig.  8.  Transition  from  low  temperature  chemistry  induced  ignition  to  a  hot  flame, 
for  stoichiometric  DME/02  with  15.8%  N2,  1  atm,  k  =  20/s. 


Fig.  9.  Low  temperature  chemistry  induced  ignition  for  stoichiometric  DME/02 
with  various  amount  of  N2,  1  atm,  k  =  50/s. 


Fig.  10.  Low  temperature  chemistry  induced  ignition  for  DME/02  with  various 
equivalence  ratios  and  15.8%  N2,  1  atm,  k  =  50/s. 

condition  is  set  for  the  boundary  velocity  while  the  boundary  tem¬ 
perature  is  set  at  300  K.  Fig.  11(a)  shows  two  cool  flames  with  dif¬ 
ferent  boundary  velocities  of  1 5  cm/s  (case  A)  and  30  cm/s  (case  B), 
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Fig.  11.  (a)  Velocity  and  temperature  profiles  for  stoichiometric  DME/02/N2  cool 
flames  with  15.8%  N2  in  a  symmetric  counterflow  with  plug  flow  boundary 
condition  under  1  atm,  showing  local  minimum  (A)  and  inflection  point  (B) 
behavior  for  the  reference  cool  flame  speed  corresponding  to  strong  and  weak 
thermal  expansion,  respectively,  (b)  Reference  cool  flame  speed  for  the  mixture  of 
Fig.  11(a)  at  different  strain  rates,  showing  a  linear  behavior  as  a  function  of  strain 
rate. 


providing  different  strain  rates,  for  stoichiometric  DME/02/N2  mix¬ 
ture  with  15.8%  N2  under  1  atm.  It  is  seen  that,  for  the  low  strain 
rate  case,  the  flame  temperature  is  higher  such  that  thermal 
expansion  is  larger,  resulting  in  a  non-monotonic  velocity  profile 
in  the  flame  zone.  The  local  strain  rate  k  characterizing  the  flame 
residence  time  can  be  defined  as  the  local  maximum  strain  rate 
upstream  of  the  flame,  and  the  local  minimum  flow  velocity  as 
the  reference  flame  speed  of  the  cool  flame  SUiC00/  for  the  corre¬ 
sponding  strain  rate,  in  the  same  manner  as  that  of  the  hot  flame 
characterization  in  the  counterflow  configuration  [31].  It  is  seen 
that  in  case  A,  the  local  strain  rate  is  about  20.7/s  and  the  local 
minimum  velocity  is  6.6  cm/s,  and  as  such  demonstrates  that  the 
reference  cool  flame  speed  5U  C00/  is  much  slower  than  the  laminar 
flame  speed  under  the  same  condition.  In  case  B  for  which  the  exit 
velocity  doubles  to  30  cm/s,  the  flame  temperature  and  the  ther¬ 
mal  expansion  effect  on  the  flow  velocity  decrease.  Consequently 
the  nonmonotonic  behavior  in  the  velocity  profile  is  reduced  to 
an  inflection  point  for  the  higher  strain  rate,  leading  to  ambiguity 
in  evaluating  the  cool  flame  speed.  For  such  cases,  the  reference 
cool  flame  speed  corresponding  to  the  local  maximum  strain  rate 
is  defined  as  the  flow  velocity  at  the  inflection  point,  as  shown  in 
Fig.  11(a).  For  the  cases  with  intermediate  local  strain  rate  between 


those  corresponding  to  cases  A  and  B,  thermal  expansion  is  also 
intermediate,  which  causes  a  plateau  behavior  in  the  correspond¬ 
ing  velocity  profile. 

We  now  systematically  show  how  a  premixed  cool  flame 
responds  to  the  strain  rate  variation.  The  reference  cool  flame 
speed  as  a  function  of  local  strain  rate  is  given  in  Fig.  11(b),  for 
the  same  mixture  in  Fig  11(a).  Four  different  regimes  could  be 
identified  by  varying  the  strain  rate.  For  a  low  enough  strain  rate 
k ,  below  20/s,  ignition  induced  by  the  low-temperature  chemistry 
has  long  enough  residence  time  to  transition  to  an  intense  burning 
hot  flame,  therefore  a  steady-state  cool  flame  does  not  exist,  as 
demonstrated  in  Fig.  8.  With  k  increased  from  20  to  30/s,  a 
steady-state  cool  flame  can  be  initiated,  with  large  enough  thermal 
expansion  to  cause  a  nonmonotonic  velocity  profile  in  the  flame 
zone  as  case  A  in  Fig.  1 1(a).  In  this  regime,  the  reference  cool  flame 
speed  at  the  minimum  point  varies  linearly  with  the  strain  rate, 
and  hence  the  flame  stretch  rate,  although  in  an  insensitive  man¬ 
ner.  With  k  further  increased  to  the  range  of  30-40/s,  thermal 
expansion  decreases  as  a  result  of  the  reduced  residence  time 
and  the  nonmonotonic  velocity  profile  reduces  to  a  plateau  behav¬ 
ior,  where  the  reference  cool  flame  speed  is  even  more  insensitive 
to  strain  rate  variations.  With  k  further  increased,  the  plateau 
behavior  is  reduced  to  the  inflection  point  behavior  as  case  B  in 
Fig.  11(a),  for  which  the  reference  cool  flame  speed  again  shows 
linear  dependence  on  k  with  the  same  slope  as  in  the  regime  with 
local  minimum  velocity.  These  calculations  hence  shed  light  on  the 
dynamic  response  of  laminar  premixed  cool  flames  to  the  variation 
of  aerodynamic  strain  and  flame  stretch. 

Fig.  12  shows  the  chemical  structure  of  the  calculated  DME/02/ 
N2  cool  flame  with  15.8%  N2  and  18  cm/s  boundary  velocity.  It  is 
seen  that  the  major  species  are  those  involved  in  the  first-stage 
autoignition  process,  and  that  both  fuel  and  oxygen  are  only 
slightly  consumed,  producing  H20,  CO,  formaldehyde,  etc.,  in  the 
downstream.  Important  radicals  in  high  temperature  flames  such 
as  H,  OH  have  rather  small  concentrations,  with  peak  values  being 
around  0.1  ppm.  To  identify  the  critical  reactions  for  the  cool  flame 
propagation,  dimensionless  sensitivity  coefficients  of  each  reaction 
to  the  reference  cool  flame  speed,  defined  as  0lnSUiCOO//SlnAj  and 
normalized  by  the  maximum  absolute  value,  are  calculated  for 
the  case  of  Fig.  9  and  shown  in  Fig.  13.  It  is  seen  that  the  cool  flame 
speed  is  mainly  controlled  by  the  low-temperature  chemical  kinet¬ 
ics.  Reactions  promoting  cool  flame  propagation  are  mainly  the 
low-temperature  chain  branching  and  propagation  reactions, 
including  H  abstraction  of  DME,  R02  isomerization  into  QOOH, 


Fig.  12.  Mole  fractions  of  major  species  for  a  stoichiometric  DME/02/N2  cool  flames 
with  15.8%  N2  in  a  symmetric  counterflow  configuration  with  inlet  velocity  18  cm/s, 
under  1  atm. 
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Fig.  13.  Normalized  sensitivity  coefficients  demonstrating  important  reactions  for 
cool  flame  propagation  corresponding  to  the  case  shown  in  Fig.  12. 

and  02  addition  to  Q00H02,  etc.  On  the  contrary,  a  major  reaction 
retarding  cool  flame  propagation  is  the  //-scission  of  the  QOOH  rad¬ 
icals  forming  the  relatively  stable  intermediate  HCHO,  which  com¬ 
petes  with  the  low-temperature  chain  branching  reactions. 

It  is  then  of  interest  to  show  the  dependence  of  the  cool  flame 
temperature  and  speed  on  the  equivalence  ratio.  For  a  fixed  local 
strain  rate  of  25/s,  the  reference  flame  speeds  and  temperatures 
of  DME/02/N2  cool  flames,  with  02  and  N2  volume  ratio  4:1  and 
varying  equivalence  ratio,  are  shown  in  Fig.  14(a).  It  is  seen  that, 
contrary  to  the  hot  flame  cases,  the  cool  flame  temperature  gradu¬ 
ally  increases  with  equivalence  ratio  and  reaches  a  limit  under  rich 
conditions.  This  result  is  consistent  with  previous  homogeneous 
autoignition  studies  [10,32],  in  that  in  a  wide  range  of  equivalence 
ratios,  the  richer  the  mixture,  the  larger  the  temperature  increase 
at  the  end  of  the  first-stage  autoignition.  The  cool  flame  tempera¬ 
ture  should  be  bounded  by  the  “ceiling  temperature”,  which  is  the 
equilibrium  temperature  of  the  R02  radical  formation  reaction  R 
+  02  =  R02.  Under  rich  conditions,  the  dilution  effect  from  nitrogen 
is  weaker  and  can  lead  to  a  higher  cool  flame  temperature.  On  the 
other  hand,  instead  of  the  strong  dependence  on  flame  tempera¬ 
ture  for  the  high-temperature  laminar  flames  [27],  the  cool  flame 
speed  does  not  follow  the  same  trend  as  the  cool  flame  tempera¬ 
ture.  Also,  it  shows  a  modest  sensitivity  in  terms  of  the  equivalence 
ratio,  and  the  range  of  the  flame  speed  is  around  6-7  cm/s.  The 
maximum  value  is  around  stoichiometric,  similar  to  that  of  the 
laminar  flame  speed.  The  maximum  HCHO  mole  fraction  is  also 
calculated  and  shown  in  Fig.  14(b),  indicating  the  increase  of 
chemical  reactivity  of  the  cool  flames  with  increasing  equivalence 
ratios. 

4.  Stretched  cool  flame  in  counterflow  -  experimental 
observations 

Experiments  were  carried  out  in  a  counterflow  facility  used  in 
[12,16,33]  to  explore  if  the  predicted  premixed  cool  flame  indeed 
exists.  Briefly,  the  counterflow  setup  consists  of  two  vertically 
opposing  quartz  tubes  with  20  mm  diameters  and  separated  by 
20  mm.  A  room-temperature  DME/02  mixture  is  issued  from  the 
lower  tube,  flowing  against  an  N2  jet,  which  is  directed  downward 
by  the  upper  tube  and  electrically  heated  by  both  inner  and  outer 
heaters.  Both  upper  and  lower  jets  are  surrounded  by  N2  guard 
flows.  The  attempt  to  observe  the  cool  flame  was  conducted  by 
either  gradually  increasing  the  N2  temperature  to  about  700  K 


Equivalence  ratio 

Fig.  14.  (a)  Flame  speeds  and  temperatures,  (b)  Maximum  HCHO  mole  fraction  for 
the  cool  flame  calculations  for  atmospheric  DME/02/N2  mixtures  with  15.8%  N2  in  a 
symmetric  counterflow  configuration  with  various  equivalence  ratios  at  300  K,  and 
the  local  strain  rate  is  fixed  around  25/s. 

under  proper  flow  rate,  or  gradually  varying  the  flow  rate  with 
fixed  N2  boundary  temperature.  Since  a  cool  flame  is  hardly 
observable  by  eye  or  camera  in  the  visible  spectrum  range,  a 
high-sensitivity  monochrome  CCD  camera  with  high  relative 
response  in  the  UV  spectrum  (CCE-B013-U,  Mightex)  was  setup 
to  capture  the  global  chemiluminescence  from  the  low  tempera¬ 
ture  chemistry.  It  is  noted  that  extensive  work  has  been  performed 
on  the  chemiluminescence  from  the  excited  species  in  flames,  such 
as  excited  HCHO  molecule  [34],  excited  C02  molecule  [34-36], 
excited  CH  radical  [36-38]  and  OH  [36,37,39]  radical.  Since  the 
main  goal  of  the  current  experiment  is  to  prove  the  existence  of 
a  steady  laminar  premixed  cool  flame  in  the  strained  flow,  the  glo¬ 
bal  chemiluminescence  from  the  low  temperature  chemistry  is 
sufficient  to  be  used  as  an  indicator  of  the  location  and  burning 
intensity  of  the  cool  flame.  Consequently,  a  quantitative  validation 
of  the  chemistry  by  separate  measurements  of  the  emission  from 
each  of  these  intermediates  was  not  conducted  and  certainly  mer¬ 
its  further  efforts. 

To  quantify  the  flow  and  temperature  boundary  conditions,  an 
uncoated  K-type  (Chromel-Alumel)  thermocouple  was  used  to 
measure  the  maximum  exit  temperature  of  the  heated  N2  and 
the  flow  rate  was  measured  by  a  flow  calibrator  (Gilibrator-2  Flow 
Calibrator,  Standard  Flow  Cell,  Sensidyne).  From  discussions  in  Sec¬ 
tion  2,  it  is  noted  that  the  dilution  effect  from  N2  in  the  mixture  can 
be  essential  for  the  existence  of  a  premixed  cool  flame  even  at  a 
fixed  strain  rate.  Consequently,  mixtures  of  DME  and  02  without 
N2  addition  were  used  to  facilitate  the  experimental  observation. 
To  prevent  flashback  of  the  possible  hot  flames  generated,  two 
porous  flame  quenchers  (with  1  pm  average  pore  size)  were 
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connected  in  between  the  flow  tubes  and  downstream  of  the  mix¬ 
ing  chamber. 

To  validate  the  observation,  we  first  performed  the  experiment  of 
a  reactive  case  with  heated  N2  (686  K)  flowing  against  a  DME/02 
mixture  of  equivalence  ratio  1.5;  shown  as  the  left  image  of 
Fig.  15.  Clearly  a  thin  layer  with  chemiluminescence  is  observed  in 
the  flow  field.  When  02  is  replaced  by  the  same  amount  of  N2  to  keep 
the  same  flow  rate  and  boundary  temperatures,  the  reactive  layer 
disappeared  as  shown  in  the  right  image  of  Fig.  15.  The  cool  flame 
would  recur  if  the  amount  of  N2  in  the  mixture  was  replaced  by  02. 
Furthermore,  if  the  amount  of  DME  was  replaced  by  the  same 
amount  of  N2,  the  observed  cool  flame  again  disappeared  and 
recurred  when  DME  was  changed  back.  From  this  behavior,  we  con¬ 
clude  that  the  observed  reactive  layer  is  induced  by  the  chemistry 
between  DME  and  02  at  low  temperatures.  In  addition,  a  pass- 
band  filter  corresponding  to  the  chemiluminescence  from  the 
excited  HCHO  molecules,  with  characteristic  spectrum  peaking 
around  400  nm  [12],  was  utilized  to  verify  the  observation.  It  is  seen 
that  the  chemiluminescence  appears  at  the  same  location,  although 
the  intensity  is  much  weaker  with  the  filter  and  therefore  is  not 
shown. 

The  flow  rates  were  also  controlled  to  ensure  that  the  experi¬ 
mental  condition  matches  the  strain  rate  window  for  the  occur¬ 
rence  of  cool  flames,  as  discussed  in  Section  2.  Specifically,  the 
flow  rate  of  the  hot  N2  (686  K)  from  the  upper  nozzle  is  about 
1826cm3/min  while  that  of  the  DME/02  mixture  with  various 
equivalence  ratios  is  fixed  at  960  cm3/min,  from  the  bottom  nozzle 
with  ambient  temperature,  300  K.  Fig.  16  shows  the  images  taken 
for  the  steady  state  counterflow  cool  flames  observed  with  DME/02 
mixtures  with  various  equivalence  ratios  from  0.8  to  1.8. 

The  above  observation  indicates  that  there  are  two  qualitative 
trends  for  the  kinetic  and  dynamic  behavior  of  the  premixed  cool 
flames.  First,  with  increasing  equivalence  ratio  from  lean  to  rich, 
the  chemiluminescence  from  the  cool  flame  becomes  progressively 
more  intense.  This  is  consistent  with  the  calculated  concentration 
of  HCHO  (Fig.  14(b)),  which  is  a  characteristic  species  indicating 
the  low-temperature  chemistry,  showing  that  it  increases  with 
equivalence  ratio  for  the  same  flow  condition,  and  hence  stronger 
chemiluminescence.  However,  it  should  be  noted  that  the  chemilu¬ 
minescence  observed  is  partially  from  the  excited  HCHO  mole¬ 
cules,  instead  of  the  ground  state  HCHO  molecules.  Since  HCHO 


is  a  major  characteristic  species  showing  the  low-temperature 
chemistry,  the  qualitative  agreement  implies  stronger  low- 
temperature  reactivity  and  chemiluminescence  with  increased 
equivalence  ratio,  even  for  rich  mixtures.  Second,  the  position  of 
the  reaction  layer  almost  does  not  change  for  the  same  boundary 
flow  rate.  This  observation  then  demonstrates  the  insensitive 
dependence  of  the  cool  flame  speed  on  the  mixture  equivalence 
ratio.  As  shown  in  Fig.  14(a),  for  a  wide  range  of  equivalence  ratios 
(0.4-2.0),  the  variation  of  the  cool  flame  speed  is  within  about  15%, 
which  is  much  less  than  the  variation  of  normal  laminar  flame 
speeds  in  terms  of  the  equivalence  ratio. 

Compared  with  the  conventional  measurements  of  laminar 
flame  speeds,  a  quantitative  measurement  of  the  cool  flame  speed 
embodies  considerable  difficulty.  Specifically,  the  local  flow  veloc¬ 
ity  is  low  and  the  reaction  layer  is  quite  weak,  and  as  such  it 
responds  sensitively  to  any  flow  disturbance  from  intrusive  meth¬ 
ods,  such  as  thermocouple  measurement.  Furthermore,  seeding 
particle  loading  in  LDV  (Laser  Doppler  Velocimetry)  or  PIV  (Particle 
Image  Velocimetry)  measurements  could  modify  the  boundary 
conditions  or  the  local  flow  field  to  cause  flame  extinction  or 
autoignition.  It  is  also  noted  that  the  observed  reaction  layer  is 
from  the  chemiluminescence  of  the  cool  flame  chemistry,  corre¬ 
sponding  to  the  position  of  the  peak  concentration  of  excited 
HCHO  molecules  in  the  downstream  of  the  cool  flame.  Therefore, 
the  observed  cool  flame  is  still  at  the  downstream  of  the  location 
for  the  reference  velocity  as  shown  in  Fig.  11. 

Fig.  17  compares  the  experimental  image  with  detailed  calcula¬ 
tion  for  the  counterflow  configuration  with  hot  N2  (upper,  686  K) 
flowing  against  DME/02  mixture  with  equivalence  ratio  1.6.  The 
position  of  the  captured  chemiluminescence  well  corresponds  to 
the  peak  of  the  simulated  temperature  profile  and  also  the  HCHO 
profile.  It  is  seen  that  the  observed  chemiluminescence  is  at  a  loca¬ 
tion  very  close  to  the  stagnation  plane,  where  the  corresponding 
flow  velocity  is  quite  low.  However,  from  the  calculation,  it  is  seen 
that  the  location  for  the  reference  cool  flame  speed  is  about  0.4  cm 
away  from  the  chemiluminescence  zone,  leading  to  large  ambigu¬ 
ity  in  the  estimation  of  the  location  of  the  reference  cool  flame 
speed  and  therefore  the  reference  cool  flame  speed.  With  thermal 
expansion  from  the  cool  flame  as  an  additional  unknown,  it  is 
therefore  not  feasible  to  meaningfully  determine  the  cool  flame 
speed  based  on  the  experimental  chemiluminescence  observation. 


Fig.  15.  Observations  for  experiments  of  hot  N2  vs.  DME/02  mixture  (left)  and  hot  N2  vs.  DME/N2  mixture  (right)  with  the  same  flow  rate  and  same  boundary  temperatures. 
The  flow  rate  of  the  hot  N2  (686  K)  from  the  upper  nozzle  is  about  1826  cm3/min  while  that  of  the  DME/02  mixture  with  various  equivalence  ratios  is  fixed  at  960  cm3/min, 
from  the  bottom  nozzle  with  ambient  temperature,  300  K. 
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Fig.  16.  Images  taken  with  a  CCD  camera  for  atmospheric  cool  flame  observations  (shown  as  the  arrow)  in  counterflow  configuration  with  hot  N2  (upper,  686  K)  flowing 
against  DME/02  mixture  with  various  equivalence  ratios  of  0.8,  1.0,  1.2,  1.4,  1.6,  1.8  (from  left  to  right  and  from  top  to  bottom). 


■=> 

dme/o2/n2 


Heated  N2 


Fig.  17.  Comparison  for  the  premixed  cool  flame  image  in  counterflow  configuration  with  hot  N2  (upper,  686  K)  flowing  against  DME/02  mixture  with  equivalence  ratios  of 
1.6  with  detailed  numerical  simulation.  Dash  line  in  the  simulation  indicates  the  location  for  the  stagnation  plane. 


In  summary,  we  have  successfully  established  and  observed  a  computations.  Furthermore,  while  it  is  beyond  the  scope  and 
steady  cool  flame  in  the  aerodynamically  strained  counterflow,  experimental  capability  of  the  present  investigation  to  probe  the 

whose  characteristics  correspond  well  to  those  obtained  from  underlying  chemistry,  we  have  nevertheless  demonstrated  the  fea- 
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sibility  of  using  the  counterflow  to  further  study  the  cool  flame 
characteristics  including  its  chemistry  with  advanced  laser  instru¬ 
mentation,  such  as  HCHO-PLIF  [40]  and  OH-PLIF  [41]. 

5.  Conclusions 

In  the  present  investigation,  we  have  first  computationally 
studied  the  initiation  and  propagation  of  premixed  cool  flames 
for  DME/02/N2  mixtures  in  the  1-D  planar  configuration  as  well 
as  the  counterflow  configuration.  The  results  demonstrate  that  a 
residence  time  window  exists  for  the  occurrence  of  cool  flames; 
external  to  this  window  a  cool  flame  either  cannot  be  initiated  or 
will  transition  to  a  hot  flame.  Parametric  study  shows  that  the  cool 
flame  regime  is  sensitive  to  the  amount  of  dilution  in  the  flow,  and 
that  the  cool  flame  becomes  more  intense  at  richer  conditions.  Glo¬ 
bal  parameters  such  as  the  reference  cool  flame  speed  and  cool 
flame  temperature  have  been  studied  in  a  symmetric  counterflow 
for  various  strain  rates  and  equivalence  ratios,  with  the  controlling 
reactions  identified.  The  results  show  that  the  premixed  cool  flame 
speeds  are  relatively  low,  around  6-7  cm/s,  and  are  much  less  sen¬ 
sitive  in  terms  of  equivalence  ratio  variations.  Experimental  inves¬ 
tigations  have  also  been  carried  out  in  a  counterflow,  with  a  steady 
cool  flame  observed  with  a  CCD  camera  with  high  sensitivity  of  the 
UV  band.  Such  observations  qualitatively  demonstrate  stronger 
cool  flame  chemiluminescence  and  insensitive  change  in  the  flame 
position  with  increasing  equivalence  ratios,  which  correspond  well 
with  the  detailed  numerical  simulation. 
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Cool  flames,  governed  by  low-temperature  chemistry,  are  closely  related  to  engine  knock.  Since  the  low- 
temperature  chemical  kinetics  is  promoted  at  elevated  pressures,  the  ignition  and  extinction  of  non¬ 
premixed  cool  flame  at  elevated  pressures  were  experimentally  and  computationally  investigated  herein 
in  the  counterflow.  Specifically,  the  hysteretic  ignition  and  extinction  behavior  of  the  nonpremixed  cool 
flame  was  for  the  first  time  observed  and  quantified.  S-curve  analysis  was  conducted  to  demonstrate  the 
thermal  and  chemical  structure  of  the  cool  flame  and  to  elucidate  the  dominant  chemical  pathways  dur¬ 
ing  the  ignition  and  extinction  processes.  The  dominant  low-temperature  chemical  reactions  shift  from 
those  responsible  for  radical  runaway  to  exothermic  reactions  that  sustain  the  cool  flame.  Increasing  the 
ambient  pressure  and/or  the  oxygen  concentration  in  the  oxidizer  stream  promote  the  heat  release  from 
the  cool  flame,  and  hence,  result  in  an  extended  hysteresis  temperature  window  between  ignition  and 
extinction.  It  is  further  noted  that  while  the  observed  cool  flame  ignition  temperatures  were  well  pre¬ 
dicted  by  computation,  significant  discrepancies  existed  for  the  extinction  temperatures  based  on  the 
well-adopted  reaction  mechanism  used.  Possible  reasons  were  discussed  to  facilitate  further  cool  flame 
studies  and  the  development  of  the  low-temperature  chemistry. 
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1.  Introduction 

Cool  flames,  first  reported  in  1817  [1],  are  controlled  by  low- 
temperature  chemical  kinetics  that  have  been  extensively  studied 
ever  since  [2].  Such  low-temperature  chemistry  is  ubiquitous  for 
most  large  hydrocarbon  fuels  and  has  been  shown  to  be  related  to 
the  negative  temperature  coefficient  (NTC)  phenomenon  observed 
in  autoignition  processes  [3,4]  and  engine  knocks  [5].  The  fun¬ 
damental  understanding  of  low-temperature  chemistry  and  cool 
flame  dynamics  can  also  be  important  for  combustion  phasing  con¬ 
trol  in  the  recent  development  of  homogeneous  charge  compres¬ 
sion  ignition  (HCCI)  engines  [6,7]  and  reactivity  controlled  com¬ 
pression  ignition  (RCCI)  [8]  engines. 

Besides  the  potential  utilization  in  engines,  the  dynamics  of  the 
cool  flame  are  of  fundamental  interest.  Cool  flames  and  the  NTC 
phenomena  have  been  mostly  observed  in  homogeneous  systems 
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such  as  rapid  compression  machines  [9-11],  shock  tubes  [4,12,13], 
flow  reactors  [14,15]  and  stirred  reactors  [16,17].  Without  the  com¬ 
plexity  of  inhomogeneity  and  transport,  low-temperature  chemical 
kinetics  can  be  studied  in  these  systems  for  chemical  model  devel¬ 
opment.  However,  in  practical  combustors  such  as  the  diesel  and 
gasoline  direct  injection  (GDI)  engines,  combustion  processes  are 
governed  by  both  transport  and  chemical  kinetics.  For  such  con¬ 
ditions,  the  characteristic  mixing  time  scales  can  be  comparable 
to  the  chemical  reaction  time  scales,  and  the  convective-diffusive 
processes  may  affect  the  initiation  and  sustenance  of  cool  flames. 

Recognizing  that  the  inevitable  presence  of  nonuniformities  in 
practical  combustion  systems  requires  consideration  of  the  cou¬ 
pled  effects  of  chemistry  and  transport  on  cool  flames,  cool  flames 
in  nonpremixed  systems  have  been  recently  studied  in  the  coun¬ 
terflow  [18-21]  and  microgravity  droplet  combustion  [22-24].  In 
particular,  Law  and  Zhao  [18]  and  Zhao  and  Law  [19]  numerically 
investigated  the  nonpremixed  counterflow  cool  flame  of  n-heptane, 
and  identified  the  existence  of  a  secondary  S-curve  dominated 
by  low-temperature  chemistry  at  low  strain  rates  and/or  high 
pressures.  This  secondary  S-curve  is  grafted  onto  the  lower  branch 
of  the  primary  S-curve  and  has  its  own  ignition  and  extinction 
turning  points.  Experimentally,  nonpremixed  cool  flames  were 
observed  in  the  counterflow  at  atmospheric  pressure  [20],  includ¬ 
ing  those  employing  a  reactivity  promoter  [21,25].  The  ignition 
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temperature  of  dimethyl  ether  (DME)  was  quantified  with  infrared 
imaging  [20],  and  was  found  to  increase  with  increasing  strain 
rate  but  was  insensitive  to  the  DME  concentration  in  the  fuel 
stream.  In  microgravity  droplet  combustion,  it  was  found  that  the 
visible  flame  of  large  n-heptane,  n- octane,  and  n-decane  droplets 
could  transit  to  a  quasi-steady  low-temperature  burning  mode 
after  extinction  of  the  hot  flame  due  to  radiation,  suggesting  the 
existence  of  steady  droplet  burning  sustained  by  a  cool  flame  [24]. 

Since  low-temperature  chemistry  is  more  pronounced  at  ele¬ 
vated  pressures,  which  for  example  could  affect  flame  stabiliza¬ 
tion  [26],  we  have  conducted  a  systematic  experimental  and  com¬ 
putational  study  on  nonpremixed  DME  cool  flames  at  elevated 
pressures.  DME  was  chosen  because  it  is  an  alternative  diesel  fuel 
that  has  low  soot  emissions  due  to  its  high  cetane  number  and 
hydrogen-carbon  ratio  [27].  Moreover,  there  has  been  much  de¬ 
velopment  of  DME  chemical  models  with  low-temperature  chem¬ 
istry  that  are  also  computationally  affordable  [28,29].  Furthermore, 
DME  is  a  gaseous  fuel,  which  does  not  require  pre-vaporization 
and  insulation  in  the  experimentation.  The  first  objective  of  the 
present  study  is  to  experimentally  observe  and  quantify  the  igni¬ 
tion  and  extinction  states  of  strained  nonpremixed  cool  flames  and 
hopefully  to  substantiate  the  computationally  predicted  ignition- 
extinction  S-curve  hysteresis.  Upon  such  a  validation,  computation 
with  detailed  chemistry  and  transport  are  therefore  able  to  elu¬ 
cidate  the  thermal  and  chemical  structures  of  the  cool  flame  and 
to  demonstrate  their  evolution  during  the  ignition  and  extinction 
processes.  Finally,  effects  of  the  ambient  pressure  and  oxygen  con¬ 
centration  are  investigated  and  compared  between  experiment  and 
computation.  Possible  reasons  for  discrepancies  in  such  compar¬ 
isons  are  discussed  to  facilitate  further  studies. 

2.  Methodology 

2.1.  Experimental  methodology 

The  ignition  and  extinction  temperatures  of  DME  nonpremixed 
cool  flames  were  measured  in  the  same  counterflow  system  as  that 
of  Deng  et  al.  [20],  but  at  elevated  pressures  and  different  bound¬ 
ary  conditions.  Detailed  descriptions  of  the  experimental  system 
are  provided  in  [30,31].  Briefly,  the  counterflow  facility  consists  of 
two  symmetrical,  vertically  oriented,  opposing  quartz  tubes  with 
inner  diameter  of  2  cm  and  separated  by  2  cm.  A  fuel  stream  con¬ 
sisting  of  50%  DME  and  50%  nitrogen  at  room  temperature  is  is¬ 
sued  from  the  lower  nozzle  and  impinges  onto  the  heated  oxidizer 
stream  from  the  upper  nozzle.  Both  streams  are  shielded  from  the 
ambience  by  coflowing  nitrogen  streams.  The  exhaust  valve  of  the 
chamber  is  adjusted  to  balance  the  inflow  and  outflow  so  as  to 
maintain  the  desired  chamber  pressure. 

A  high-sensitivity  monochrome  CCD  camera  with  high  relative 
response  in  the  UV  spectrum  (CCE-B013-U)  was  used  to  capture 
the  global  chemiluminescence  from  the  low-temperature  chemi¬ 
cal  reactions,  which  indicated  the  intensity  of  the  cool  flame  [32]. 
The  exposure  time  of  the  UV  camera  is  the  same  6  s  for  all  cases. 
The  ignition/extinction  states  of  the  cool  flame  were  determined  by 
gradually  changing  the  oxidizer  boundary  temperature  for  a  given 
fuel/oxidizer  flow  rate,  until  the  chemiluminescence  signal  of  the 
cool  flame  respectively  emerged/disappeared  as  detected  by  the 
camera.  The  ignition  and  extinction  temperatures  were  quantified 
with  the  corresponding  oxidizer  boundary  temperatures  measured 
with  an  uncoated  K-type  (Chromel-Alumel)  thermocouple  after  ra¬ 
diation  correction  [30].  To  measure  the  oxidizer  temperature  with¬ 
out  the  presence  of  cool  flame,  oxygen  was  replaced  with  nitrogen 
while  the  flow  rates  and  heating  power  were  kept  the  same. 

The  ignition  and  extinction  states  of  the  cool  flames  were  stud¬ 
ied  at  2  and  3  atm  ambient  pressures,  with  the  oxygen  volume 
fraction  in  the  oxidizer  stream  varied  from  21%  to  25%.  The  global 


strain  rate  used  in  this  study  is  defined  as  the  pressure-weighted 
gradient  of  the  axial  flow  velocity  [33]. 

2.2.  Computational  methodology 

The  S-curve  analysis  [34],  in  which  a  system  response  such  as 
the  maximum  temperature  or  maximum  radical  concentration  is 
studied  versus  the  variation  of  an  imposed  parameter,  was  adopted 
here  to  investigate  the  steady-state  response  of  the  counterflow 
reactive  system  subjected  to  the  physical  effects  of  flow  strain  or 
heat  loss.  The  upper  and  lower  turning  points  of  the  S-shaped  re¬ 
sponse  curve  denote  the  system  extinction  and  ignition  states,  re¬ 
spectively. 

The  S-curve  analysis  was  conducted  computationally.  The  gov¬ 
erning  equations  for  the  counterflow  flame  computation  are  pro¬ 
vided  in  [35].  A  posteriori  analysis  based  on  the  optically  thin  radi¬ 
ation  model,  with  radiative  properties  based  on  the  RADCAL  model 
by  Grosshandler  of  NIST  [36]  shows  that  the  radiative  heat  loss 
is  minimal,  and  therefore  radiation  is  not  included  in  the  current 
computation.  The  damped  Newton  method  and  time  integration 
solution  scheme  were  adopted  in  a  modified  numerical  code  of 
Smooke  et  al.  [37],  solving  the  differential  equations  with  bound¬ 
ary  conditions  specified  on  both  sides  of  the  potential  flow.  S-curve 
marching  is  conducted  using  a  flame  controlling  method  [38].  A 
skeletal  DME  mechanism  [29],  reduced  from  the  detailed  mech¬ 
anism  of  Zhao  et  al.  [28]  with  detailed  transport  properties,  was 
used.  Details  about  the  reduction  method  and  validation  results 
against  the  full  mechanism  under  a  wide  range  of  conditions  were 
provided  as  supplementary  material  of  that  work. 

In  the  computation,  the  nozzle  separation  distance  was  held 
at  2  cm  to  match  the  experimental  system.  The  temperature  and 
composition  boundary  conditions  were  consistent  with  the  exper¬ 
iments.  To  generate  the  S-curve,  the  boundary  temperature  of  the 
heated  oxidizer  stream  was  changed  at  a  given  strain  rate.  The  crit¬ 
ical  ignition  and  extinction  states  of  the  cool  flame  were  predicted 
based  on  the  turning  points  of  the  secondary  S-curve  [18,20]. 

3.  Results  and  discussion 

In  the  following  discussion,  the  experimental  affirmation  of  the 
distinctive  ignition  and  extinction  behavior  of  cool  flame  is  first 
demonstrated,  and  thereby  validate  the  computational  prediction. 
Furthermore,  the  thermal  and  chemical  structures  of  the  reacting 
layer  in  the  counterflow  at  cool  flame  ignition,  steady  burning,  and 
extinction  conditions  are  demonstrated,  elucidating  the  dominant 
chemical  pathways  for  each  condition.  Finally,  the  effects  of  ambi¬ 
ent  pressure  and  oxygen  concentration  on  ignition  and  extinction 
are  discussed. 

3.1.  Hysteretic  ignition  and  extinction  behavior 

What  differentiates  a  strained  cool  flame  from  a  reacting  mix¬ 
ture  going  through  slow  oxidation  in  a  heated  flow  is  the  hys¬ 
teretic  ignition  and  extinction  behavior.  Such  hysteretic  behavior 
is  illustrated  with  an  S-curve  analysis  shown  in  Fig.  1.  Since  the 
hot  ignition  of  DME  in  the  counterflow  configuration  has  already 
been  reported  in  the  literature  [39],  only  the  cool  flame  portion 
of  the  S-curve  is  included.  At  a  fixed  strain  rate,  by  gradually  in¬ 
creasing  the  air  boundary  temperature  and  hence  the  reactivity 
and  heat  generation  rate  within  the  flow,  at  some  point  within  the 
flow  the  temperature  will  exceed  the  boundary  temperature  and 
eventually  lead  to  self-sustained  burning.  Such  a  transition  state 
I  is  defined  as  the  ignition  point,  and  the  condition  after  ignition 
that  is  located  on  the  upper  steady  branch  of  the  S-curve  is  desig¬ 
nated  as  point  S,  representing  a  steady  cool  flame.  The  difference 
in  the  maximum  temperatures  between  points  S  and  I  primarily 
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Fig.  1.  S  -Curve  analysis  of  the  response  of  nonpremixed  counterflow  of  50%  DME 
and  50%  nitrogen  versus  heated  air  at  the  2  atm  and  pressure-weighted  strain  rate 
of  80  s-1. 


7’air=  626  K 


7^  =  629  K 


7^  =  635  K 
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Fig.  2.  Flame  images  demonstrating  the  hysteretic  nature  of  ignition  and  extinc¬ 
tion  of  cool  flames  with  air  temperature.  DME  cool  flame  at  2  atm  and  pressure- 
weighted  strain  rate  of  84  s-1 ;  DME  volume  fraction  in  the  fuel  stream  is  50%,  and 
the  oxidizer  stream  is  air. 


results  from  the  heat  release  from  the  cool  flame.  As  this  tempera¬ 
ture  decreases  from  point  S,  the  maximum  temperature  in  the  flow 
field  also  decreases,  following  the  upper  branch  trajectory,  until  the 
flame  extinguishes  at  point  E.  The  difference  in  the  air  boundary 
temperature  between  points  I  (or  S)  and  E  demonstrates  the  hys¬ 
teresis  between  the  ignition  and  extinction  of  the  flame. 

To  capture  this  computationally  predicted  hysteresis,  the  igni¬ 
tion  and  extinction  temperatures  were  experimentally  measured 
in  the  counterflow  based  on  the  chemiluminescence  of  the  cool 
flame.  Figure  2  shows  representative  images  of  the  ignition  and 
extinction  detection  for  the  DME/air  cool  flame  at  2  atm  and 
pressure-weighted  strain  rate  of  84  s-1.  Ignition  of  the  cool  flame 
is  achieved  by  gradually  increasing  the  oxidizer  boundary  temper¬ 
ature.  As  shown  in  the  left  column  of  Fig.  2,  when  the  oxidizer 
boundary  temperature  is  slowly  increased  from  626  K  to  637  K 
(A-D),  the  chemiluminescence  from  the  low-temperature  chem¬ 
istry  suddenly  becomes  detectable  by  the  UV  camera,  indicating 
onset  of  the  cool  flame.  The  oxidizer  boundary  temperature  is 


Fig.  3.  Experimental  and  computed  ignition  and  extinction  temperatures  at  2  atm 
and  various  strain  rates.  The  computed  extinction  temperature  is  shifted  up  by 
250  K  for  better  illustration.  DME  volume  fraction  in  the  fuel  stream  is  50%,  and 
the  oxidizer  stream  is  air. 


then  gradually  reduced  with  a  temperature  interval  of  1  K,  while 
maintaining  quasi-steadiness  at  each  step.  The  right  column  of 
Fig.  2  then  shows  that  although  the  cool  flame  chemiluminescence 
remains  as  the  boundary  temperature  is  decreased  from  637  K  to 
626  K  (E  to  H),  it  suddenly  vanishes  at  626  K,  which  was  thus 
defined  as  the  extinction  temperature. 

The  same  procedure  was  followed  to  obtain  the  ignition  and  ex¬ 
tinction  temperatures  at  various  strain  rates  for  comparisons  with 
the  computationally  predicted  values,  as  shown  in  Fig.  3.  Results 
from  experiment  and  computation  then  both  show  that  the  igni¬ 
tion  and  extinction  temperatures  increase  with  increasing  strain 
rate,  due  to  reduced  residence  time.  Furthermore,  the  ignition  tem¬ 
perature  not  only  is  higher  than  the  extinction  temperature  at  a 
given  strain  rate,  it  is  also  less  sensitive  to  the  strain  rate,  resulting 
in  a  less  pronounced  ignition-extinction  hysteresis.  Noting  that  the 
repeatability  of  the  measurement  is  within  2  K,  which  is  within  the 
marker  size  in  the  figure,  and  the  error  bar  represents  the  uncer¬ 
tainty  of  the  radiation  correction  using  different  models,  compari¬ 
son  between  the  experimental  data  should  be  made  based  on  the 
upper  or  lower  bound  of  the  uncertainty  bar  across  all  the  mea¬ 
surements. 

It  is  also  apparent  from  Fig.  3  that  while  the  experimental  and 
computational  results  separately  exhibit  the  anticipated  physics, 
namely  the  extinction  temperature  is  lower  than  the  ignition  tem¬ 
perature,  and  they  both  increase  with  increasing  strain  rate,  the 
quantitative  comparison  between  them  is  overall  poor,  both  in 
magnitude  as  well  as  the  strain-rate  sensitivity.  We  shall  defer  the 
detailed  analysis  of  these  discrepancies  to  Section  3.3. 

3.2.  Analysis  of  thermal  and  chemical  structures 

In  order  to  elucidate  the  dominant  chemical  pathways  and  evo¬ 
lution  of  the  ignition  and  extinction  processes,  the  three  character¬ 
istic  points  on  the  S-curve  of  Fig.  1,  which  respectively  represent 
the  states  prior  to  ignition  (point  I),  steady  burning  (point  S),  and 
prior  to  extinction  (point  E),  were  chosen  for  structural  analysis. 

As  shown  in  Fig.  4,  the  heat  release  just  prior  to  the  initia¬ 
tion  of  the  cool  flame  is  negligible,  and  the  maximum  tempera¬ 
ture  in  the  flow  field  is  set  by  the  oxidizer  boundary.  For  both 
steady  flames  at  points  S  and  E,  there  are  reaction  kernels  delin¬ 
eated  by  the  heat  release  profiles.  Based  on  the  full  width  at  half 
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Fig.  4.  Thermal  structures  in  the  computed  flow  field  at  three  representative  points 
on  the  S-curve  shown  in  Fig.  1.  The  solid  lines  refer  to  heat  release  rate,  while  the 
dashed  lines  refer  to  temperature.  Fuel  stream  boundary  conditions  are  specified  at 
x  =  0  in  the  computation. 

maximum  of  the  heat  release  profile,  a  typical  cool  flame  thickness 
is  about  0.05  cm  for  the  present  conditions.  At  point  E,  the  heat  re¬ 
lease  rate  of  the  cool  flame  is  higher  than  that  at  point  S,  for  low- 
temperature  chemistry  is  favored  at  relatively  lower  temperatures. 
However,  the  steep  gradient  of  the  heat  release  profile  indicates 
that  heat  loss  due  to  strain  at  point  E  is  more  pronounced  com¬ 
pared  to  point  S.  As  the  boundary  temperature  further  decreases, 
such  heat  loss  increases,  and  chemical  reactions  cannot  keep  up 
with  heat  loss  from  the  reaction  zone,  leading  to  extinction  due  to 
the  limited  residence  time  in  the  strained  flow. 

The  profiles  of  three  representative  species,  namely 
methoxymethylperoxy  (CH30CH202),  formaldehyde  (CH20),  and 
hydroxyl  (OH),  were  then  investigated  to  elucidate  the  chemical 
structures  at  these  three  points,  I,  S,  and  E,  as  shown  in  Fig.  5. 
The  methoxymethylperoxy  radical  was  chosen  because  it  is  a 
representative  species  for  the  low-temperature  chemistry  [40,41]. 
Comparing  the  species  profiles  of  the  three  states,  it  is  seen  that 
the  CH30CH202  mass  fraction  increases  upon  flame  initiation, 
and  similar  to  the  heat  release  rate  profile,  its  peak  increases 
as  extinction  is  approached.  Moreover,  the  profile  broadens  for 
low-temperature  chemistry  is  favored  at  reduced  oxidizer  bound¬ 
ary  temperature.  Formaldehyde  was  selected  because  it  is  a 
major  product  of  the  cool  flame,  with  its  intensity  manifested 
through  its  chemiluminescence  [32].  It  is  then  seen  that:  neg¬ 
ligible  CH20  is  formed  prior  to  ignition;  significant  amount  is 
formed  in  the  steady  flame;  and  the  concentration  and  hence 
the  intensity  of  chemiluminescence  is  reduced  as  extinction  is 
approached.  Finally,  the  hydroxyl  radical  was  chosen  because  it 
represents  high-temperature  flame  chemistry  and  is  also  an  im¬ 
portant  radical  formed  during  the  chain  branching  reactions  of  the 
low-temperature  chemistry.  It  is  then  noted  that  the  peak  mass 
fraction  of  OH  is  several  orders  of  magnitude  smaller  than  those 
in  a  typical  hot  flame  [41].  The  ignition  of  cool  flame  is  initiated 
with  a  peak  formation  of  OH,  which  subsequently  decreases  upon 
ignition.  With  subsequent  heat  release  from  the  cool  flame,  a 
second  peak  of  OH  mass  fraction  emerges  on  the  oxidizer  side  of 
the  cool  flame  peak,  which  is  responsible  for  the  initiation  of  the 
hot  flame  at  higher  boundary  temperatures  [18]. 

Finally,  sensitivity  analysis  was  conducted  to  elucidate  the 
evolution  of  the  dominant  chemical  pathways  during  the  ignition 
and  extinction  processes.  The  maximum  temperature  was  chosen 
as  the  target  for  sensitivity  analysis,  and  the  ratio  of  the  relative 
change  of  the  maximum  temperature  to  that  of  the  Arrhenius 


Fig.  5.  (a)  Methoxymethylperoxy  radical,  (b)  formaldehyde,  and  (c)  hydroxyl  radical 
profiles  in  the  computed  flow  field  at  three  representative  points  on  the  S-curve 
shown  in  Fig.  1 .  Fuel  stream  boundary  conditions  are  specified  at  x  =  0  in  the  com¬ 
putation. 
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ch3och3  +  ho2  <=>  ch3och2  +  h2o2 
ch3och2  +  o2  <=>  ch3och2o2 

CH3  +  H02  <=>  CH30  +  OH 

ch3och2  <=>  ch2o  +  ch3 

CH20  +  OH  <=>  HCO  +  H20 

H  +  02  (+M)  <=>  H02  (+M) 

2CH3OCH202  <=>  o2  +  ch3ocho  +  CH3OCH2OH 
H02  +  OH  <=>  HzO  +  02 
CH20CH202H  <=>  OH  +  2CH20  - 
2CH3OCH202  <=>  02  +  2CH3OCH20 
CH3OCH3  +  OH  <=>  CH3OCH2  +  h2o 

ho2ch2ocho  <=>  och2ocho  +  OH 
ch2och2o2h  +  o2  <=>  o2ch2och2o2h 
ch3och2o2  <=>  ch2och2o2h  ( 

-L0 

Normalized  Sensitivity  Coefficient 

Fig.  6.  Reaction  sensitivity  analysis  at  three  representative  points  on  the  S-curve 
shown  in  Fig.  1. 

factor  of  each  reaction  being  perturbed  was  defined  as  the  sen¬ 
sitivity  coefficient.  The  sensitivity  coefficients  of  all  the  reactions 
considered  in  the  chemical  mechanism  were  normalized  by  the 
maximum  value  such  that  the  normalized  sensitivity  coefficient  is 
bounded  by  -1  and  1  to  allow  for  comparisons  between  the  cases. 

The  normalized  sensitivity  coefficients  ranked  by  their  ab¬ 
solute  values  for  states  I,  S,  and  E  in  Fig.  1  are  shown  in 
Fig.  6.  At  all  three  states,  low-temperature  chemistry  is  dom¬ 
inant,  with  the  absence  of  the  typical  high-temperature  chain 
branching  reactions  such  as  H  +  02  <^=>-  OH  +  0.  Similar  to  the 
ignition  analysis  in  Deng  et  al.  [20],  the  isomerization  reac¬ 
tion  CH30CH202  <=>  CH20CH202H  and  the  oxygen  addition  re¬ 
action  CH20CH202H  +  02  02CH20CH202H  promote  ignition, 

while  the  reaction  CH20CH202H  <^=>-  OH  +  2CH20  retards  it.  The 
species  profiles  demonstrated  in  Fig.  5  also  support  such  findings, 
for  the  mass  fraction  of  CH30CH202  for  the  extinction  state  E  is 
higher  than  that  for  the  steady  state  S,  while  the  opposite  trend 
is  observed  for  CH20.  Near  the  extinction  state  E,  low-temperature 
chemistry  is  favored  as  the  boundary  temperature  decreases,  and 
therefore,  more  CH30CH202  and  less  CH20  is  formed.  However,  for 
the  steady  cool  flames  corresponding  to  points  S  and  E,  the  relative 
importance  of  these  reactions  shifts.  First,  CH20CH202H  +  02  <^=>- 
02CH20CH202H  becomes  even  more  important  in  steady  cool 
flames.  Since  the  maximum  temperature  was  chosen  as  the  target 
to  evaluate  the  sensitivity,  the  exothermic  reaction  CH20CH202H  + 
02  <<=>>  02CH20CH202H  demonstrates  a  large  positive  sensitivity 
coefficient.  Conversely,  at  ignition,  reactions  that  finally  lead  to 
radical  runaway  are  more  important,  for  heat  release  from  the  mix¬ 
ing  layer  is  negligible  at  this  state,  as  shown  in  Fig.  4.  Second, 
the  isomerization  reaction  CH30CH202  <<=>■  CH20CH202H  becomes 
important  again  near  extinction,  for,  at  extinction,  the  radical  pro¬ 
duction  rates  barely  keep  up  with  the  transport  losses  and  the 
chain  carrying  limiting  step  becomes  crucial  again. 

3.3.  Uncertainty  analysis 

As  noted  in  Fig.  3,  while  the  computation  is  able  to  capture 
the  experimental  observation  of  the  hysteretic  feature  of  ignition 
and  extinction,  as  well  as  the  increasing  trend  of  ignition  and  ex¬ 
tinction  temperatures  with  increasing  strain  rate,  the  quantitative 
agreement  is  rather  poor.  Specifically,  the  ignition  temperature  is 
slightly  overpredicted;  its  sensitivity  to  strain  rate,  indicated  by  the 
slope  of  the  ignition  temperature  profile,  is  substantially  underpre¬ 
dicted;  and,  most  importantly,  the  extinction  temperature  is  signif¬ 
icantly  underpredicted.  Upon  extensive  exploration  of  the  various 
experimental  and  modeling  factors  that  could  contribute  to  such 
substantial  disagreements,  the  uncertainty  of  the  low-temperature 
chemistry  used  in  the  computation  has  surfaced  to  be  the  domi¬ 
nant  factor. 


Fig.  7.  Revisit  of  Fig.  3  with  modified  reaction  rates  in  the  chemical  mechanism. 
Note  that  the  line  plots  with  the  modified  reaction  rates  are  presented  without 
shifting. 


To  demonstrate  the  effects  of  such  uncertainty  and  since  the 
ignition  temperatures  of  the  cool  flames  are  fairly  well  captured, 
it  is  reasonable  to  inspect  reactions  with  large  sensitivity  coef¬ 
ficients  in  Fig.  6  for  the  cool  flames.  This  leads  to  the  identi¬ 
fication  of  the  reactions:  CH20CH202H  +  02  <<=>•  02CH20CH202H 
and  CH20CH202H  OH  +  2CH20.  Since  cool  flames  appear  to 
be  more  robust  to  extinguish,  we  have  repeated  the  calculation  by 
modifying  the  preexponential  factors  of  these  two  reactions  by  re¬ 
spectively  reducing  and  increasing  to  70%  and  200%  of  their  origi¬ 
nal  values. 

Figure  7  then  shows  that  both  qualitative  and  quantitative 
agreements  between  experiments  and  computations  are  achieved 
with  the  modified  reaction  parameters.  Specifically,  although  the 
magnitude  of  the  ignition  temperature  is  not  very  sensitive  to  the 
modifications,  as  expected,  its  sensitivity  to  increasing  strain  rate 
is  enhanced,  for  the  modifications  essentially  slow  down  the  low- 
temperature  chemistry  such  that  the  sensitivity  to  finite  residence 
time  is  more  pronounced.  More  importantly,  it  is  seen  that  good 
agreement  is  also  achieved  for  the  highly  sensitive  extinction  tem¬ 
peratures,  significantly  boosting  their  values  but  without  changing 
the  sensitivity  to  the  strain  rate  variation. 

While  the  above  results  appear  to  be  encouraging,  we  hasten 
to  clarify  that  we  are  not  suggesting  modified  kinetic  parameters 
for  certain  reactions.  What  we  have  demonstrated  is  the  sensitive 
nature  and  potential  uncertainty  of  the  low-temperature  chem¬ 
istry,  in  that  substantial  change  in  the  global  response  can  re¬ 
sult  from  even  small  changes  in  these  preexponential  factors.  It 
is  noted  that  while  the  Zhao  et  al.  model,  based  on  which  the 
adopted  skeletal  model  was  developed,  was  still  the  most  well 
accepted,  additional  experimental  data  have  become  available  re¬ 
cently  to  guide  the  modification  of  the  low-temperature  chem¬ 
ical  pathways  and  kinetic  parameters  [42,43].  It  is  nevertheless 
also  noted  that,  during  the  development  of  the  low-temperature 
chemical  model,  the  validation  data  is  limited  to  those  from  ho¬ 
mogeneous  systems.  The  current  study  aims  at  providing  valida¬ 
tion  targets  in  inhomogeneous  systems  by  substantiating  the  exis¬ 
tence  of  the  nonpremixed  cool  flame  and  revealing  the  evolution  of 
low-temperature  chemistry  during  the  ignition  and  extinction  pro¬ 
cesses.  Clearly,  additional  validation  data  on  cool  flames  covering 
a  wider  range  of  conditions,  which  are  inherently  present  in  flame 
systems,  is  needed. 
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Fig.  8.  S  -curve  analysis  for  various  ambient  pressures  and  pressure-weighted  strain 
rate  of  100  s_1.  DME  volume  fraction  in  the  fuel  stream  is  50%,  and  the  oxidizer 
stream  is  air. 


Another  factor  that  could  potentially  affect  the  accuracy  of  the 
comparison  between  the  experimental  and  computational  results  is 
the  experimental  detection  limit  of  the  UV  camera,  which  could  be 
too  high  compared  to  the  chemiluminescence  emission  of  the  cool 
flame  near  extinction.  Since  the  hysteresis  temperature  window 
was  observed  in  the  experiment  at  various  conditions,  the  detec¬ 
tion  threshold  of  the  UV  camera  should  be  lower  than  the  chemilu¬ 
minescence  intensity  of  the  steady  cool  flame  upon  ignition.  Con¬ 
sequently,  the  ignition  temperatures  for  these  cases  should  be  well 
captured.  However,  the  lower  bound  of  the  hysteresis  temperature 
window  could  be  limited  by  the  detection  threshold,  and  there¬ 
fore,  the  actual  extinction  temperatures  might  not  be  captured  ac¬ 
curately. 

It  is  further  noted  that  Shan  and  Lu  conducted  a  computational 
study  on  the  ignition  and  extinction  of  DME/air  mixtures  at  high 
pressures  in  a  perfectly  stirred  reactor  (PSR)  [44].  In  their  work, 
eigen-analysis  of  the  chemical  Jacobian  matrix  was  performed,  and 
physically  unstable  segments  were  observed  on  the  branches  of  an 
S-curve  for  DME/air  corresponding  to  both  strong  and  cool  flames. 
Oscillations  were  observed  in  the  simulation  for  the  unstable  cases 
under  perturbation.  In  the  present  experiments,  such  oscillations 
were  not  observed  when  extinction  occurred  even  in  the  presence 
of  unavoidable  system  perturbations  such  as  fluctuations  from  the 
heat  supply  and  flow  meters,  with  the  repeatability  of  ignition  and 
extinction  temperature  measurements  within  2  K,  as  indicated  in 
Section  3.1.  Consequently,  oscillation  was  not  observed  when  ex¬ 
tinction  occurred  in  the  current  study.  Therefore,  we  do  not  expect 
that  the  large  discrepancy  between  the  current  experimental  and 
computational  cool  flame  extinction  results  is  due  to  the  instabil¬ 
ity  of  the  cool  flame. 

3.4.  Effects  of  pressure  and  oxygen  concentration 

The  effects  of  ambient  pressure  and  oxygen  concentration  in 
the  oxidizer  stream  on  the  ignition  and  extinction  of  the  cool 
flames  were  also  investigated.  Pressure  effects  were  first  compu¬ 
tationally  studied  by  fixing  the  oxygen  mole  fraction  in  the  oxi¬ 
dizer  stream  at  21%  and  fixing  the  pressure-weighted  strain  rate 
at  100  s-1,  as  shown  in  Fig.  8.  It  is  seen  that,  as  the  pressure  in¬ 
creases  from  2  to  5  atm,  the  ignition  temperature  decreases  and 
the  heat  release  from  the  cool  flame  becomes  more  pronounced, 
as  indicated  by  the  temperature  differences  between  the  igni¬ 
tion  turning  point  and  the  point  on  the  cool  flame  branch  with 


Fig.  9.  Ignition  and  extinction  temperatures  at  various  strain  rates  and  pressures 
in  experiments  and  computations  with  (a)  the  original  and  (b)  modified  chemical 
models.  Some  of  the  computed  extinction  temperatures  are  shifted  up  for  better 
illustration.  DME  volume  fraction  in  the  fuel  stream  is  50%,  and  the  oxidizer  stream 
is  air. 


the  same  boundary  temperature.  Moreover,  the  extinction  turning 
point  shifts  to  a  lower  boundary  temperature,  resulting  in  an  ex¬ 
tended  hysteresis  temperature  window.  Conversely,  the  extent  of 
the  low-temperature  chemistry  governed  S-curve  hysteresis  dimin¬ 
ishes  with  decreasing  pressure,  leading  to  its  absence  at  0.5  atm. 

The  conclusion  that  elevated  pressure  promotes  low- 
temperature  chemistry  is  consistent  with  previous  studies 
with  n-heptane  [4,18]  and  more  generally  discussed  in  Pilling’s 
book  [45].  Such  promotion  effect  can  be  explained  with  the 
sensitivity  analysis  in  Section  3.2.  Qualitatively,  sensitivity  anal¬ 
ysis  conducted  for  elevated  pressures  shows  similar  dominant 
chemical  pathways  as  Fig.  6.  At  elevated  pressures,  the  balance 
of  the  reaction  CH20CH202H  +  02  <=>  O2CH2OCH2O2H  shifts 
forward  and  promotes  the  formation  of  the  important  inter¬ 
mediate  radicals  for  low-temperature  chemistry.  Moreover,  the 
CH20CH202H  OH  +  2CH20  reaction  is  retarded  at  elevated 
pressures. 

Figure  9  further  shows  that  the  experimental  ignition  and  ex¬ 
tinction  temperatures  decrease  at  elevated  pressures.  It  is  seen 
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Fig.  10.  S-curve  analysis  for  various  oxygen  concentrations  and  the  pressure- 
weighted  strain  rate  of  100  s_1.  DME  volume  fraction  in  the  fuel  stream  is  50%, 
and  the  ambient  pressure  is  2  atm. 


that,  while  the  effects  of  elevated  pressure  on  ignition  tempera¬ 
tures  are  predicted  by  computations  qualitatively,  the  effects  on 
extinction  temperatures  are  significantly  overpredicted  when  us¬ 
ing  the  original  mechanism  (Fig.  9a).  The  comparison  is  improved 
by  using  the  two  modified  reactions  (Fig.  9b),  although  the  degree 
of  improvement  is  less  satisfactory  as  for  the  case  of  2  atm  pres¬ 
sure,  shown  in  Fig.  7.  We  emphasize  again  that  we  prefer  leaving 
the  comparison  as  is,  without  further  tuning  the  reactions  as  we 
do  not  believe  it  is  justified  within  the  scope  of  the  present  study. 

Figure  10  shows  that  increasing  the  oxygen  concentration  ex¬ 
tends  the  hysteresis  temperature  window  of  the  cool  flame  at  the 
same  ambient  pressure  and  strain  rate,  while  the  ignition  tempera¬ 
ture  is  almost  unaffected  except  at  very  low  concentrations.  This  is 
because,  with  increased  oxygen  concentration  and  hence  decreased 
inert  concentration,  the  heat  release  from  the  cool  flame  becomes 
more  pronounced,  which  results  in  decreased  extinction  tempera¬ 
ture.  Dominant  chemical  pathways  for  these  conditions  are  similar 
to  those  in  Fig.  6. 

The  insensitivity  of  the  cool  flame  ignition  temperature  is  fur¬ 
ther  confirmed  with  the  experimental  measurements,  as  shown  in 
Fig.  11a.  However,  the  reduction  effect  of  increased  oxygen  con¬ 
centration  on  the  cool  flame  extinction  temperature  is  again  over¬ 
predicted  by  the  computation,  while  improved  agreements  are 
achieved  with  the  two  modified  preexponential  factors  (Fig.  lib). 

4.  Conclusions 

The  ignition  and  extinction  of  nonpremixed  DME/air  cool  flames 
at  elevated  pressures  were  experimentally  and  computationally  in¬ 
vestigated  in  the  counterflow.  For  the  first  time,  the  hysteretic  ig¬ 
nition  and  extinction  behavior  of  the  nonpremixed  cool  flame  was 
experimentally  observed  and  quantified.  Results  further  show  that 
although  low-temperature  chemistry  is  crucial  for  the  initiation 
and  sustain  of  the  cool  flame,  the  dominant  chemical  pathways 
shift  from  reactions  responsible  for  low-temperature  radical  run¬ 
away  to  cool  flame  heat  release  reactions  upon  ignition.  The  heat 
release  from  the  cool  flame  is  able  to  sustain  itself  at  lower  oxi¬ 
dizer  boundary  temperature,  and,  therefore,  results  in  the  hystere¬ 
sis  temperature  window  between  ignition  and  extinction. 

Increasing  ambient  pressure  and/or  oxygen  concentration  in  the 
oxidizer  stream  promotes  the  heat  release  from  low-temperature 
chemistry  and  extends  the  hysteresis  between  ignition  and 


Pressure- Weighted  Global  Strain  Rate  [s’1] 


Fig.  11.  Ignition  and  extinction  temperatures  at  various  strain  rates  and  oxygen 
concentrations  in  experiments  and  computations  with  (a)  the  original  and  (b)  mod¬ 
ified  chemical  models.  Some  of  the  computed  extinction  temperatures  are  shifted 
up  for  better  illustration.  DME  volume  fraction  in  the  fuel  stream  is  50%,  and  the 
ambient  pressure  is  2  atm. 


extinction.  Although  the  influences  on  the  cool  flame  ignition 
temperature  were  well  predicted  by  computation,  the  influences 
on  extinction  were  significantly  overpredicted.  Possible  reasons 
for  such  discrepancies  were  discussed,  including  the  uncertainties 
from  experiments  and  chemical  models.  The  need  for  improved 
comprehensiveness  of  the  chemical  kinetics  model  is  emphasized. 
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